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NOMENCLATURE 


A  -  Arrhenius  frequency  factor  or  prefactor 

AP  -  Ammonium  perchlorate 

a  -  Distance  from  center  of  flame  to  edge  of  oxidizer  crystal 

BDP  -  Beckstead,  Derr,  Price  combustion  model 

b  -  Distance  from  center  of  the  flame  to  outer  edge  of  burner 

C  -  Specific  heat  of  the  gas 

Cg  -  Specific  heat  of  the  solid 

C  -  Constant  (also  c) 

D  -  Oxidizer  crystal  diameter 

DQ  -  Actual  oxidizer  crystal  diameter 

D'  -  Statistically  average  intersection  diameter 

D  -  Diffusivity 

E  *  Activation  Energy 

F  -  Distribution  function 

fra  -  Fraction  of  binder  reacting 

GDF  -  Granular  diffusion  flame  model 

HMX  -  Cyclotetramethylenetetranitramine ,  C^HgNgOg 

HTPB  -  Polybutadiene 

AH  -  Heat  of  reaction 

h  -  See  Figure  6,  measure  of  non-planarity  of  oxidizer  crystal 

-  Bessel  function  of  the  ith  kind 
K  -  Kelvin 

K  -  Reaction  rate  constant  (also  k) 

L  -  See  Figure  6,  Oxidizer  projection  distance 

m  -  Mass  flow  rate 

n  -  Mass  interaction  parameter 

P  -  Pressure  (also  p) 

PEM  -  Petite  Ensemble  Model 

Q  -  Energy  released  and/or  absorbed 

R  -  Gas  constant 

R  -  See  Figure  8,  oxidizer  burning  surface  radius 

r  -  Burn  rate 

r  -  Mean  burn  rate 

r  -  Radial  coordinate 

S  -  Surface  area 

T  -  Temperature 

t  -  Time 

t^a  -  Constant  in  binder  burr.through  time  equation 

u  -  Gas  velocity 

V  -  Volume 

v  -  Gas  velocity 

W  -  Molecular  Weight 

w  -  Reaction  rate 
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x  -  Vertical  coordinate 

x  -  Flame  standoff  distance 

x  -  Mean  flame  standoff  distance 

XN  -  Diameter  exponent  in  volume  of  fuel  equation 
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CHAPTER  I 
INTRODUCTION 


A  composite  solid  propellant  typically  consists  of  oxidizer  particles, 
solid  metal  additives,  burning  rate  modifiers,  catalysts  and  acoustic 
suppressants  all  well  mixed  within  a  hydrocarbon  liquid  fuel  binder  which, 
after  mixing,  is  polymerized  to  a  solid.  Composite  solid  propellants  are  used 
in  rocket  motors  ranging  in  size  from  shoulder  launched  tactical  weapons  to 
large  space  vehicle  solid  boosters.  Efforts  to  increase  solid  motor 
performance  has  lead  to  the  use  of  nitramine  high  energy  ingredients.  The 
introduction  of  these  ingredients,  often  explosive  in  pure  form,  into 
composite  propellants  offers  improved  specific  impulse,  ISp,  which  translates 
into  greater  range  and/or  payload  for  a  solid  rocket  system. 

Combustion  modeling  can  be  used  to  tailor  propellant  ingredients  in  such  a 
fashion  as  to  achieve  desired,  or  optimum,  steady  and  non-steady  burning 
characteristics.  In  addition,  propellant  temperature  sensitivity  and  erosive 
burning  effects  can  be  studied.  Without  a  theoretical  combustion  model; 
experimental  correlations,  costly  experimental  results,  past  experience,  and 
trial  and  error  methods  are  all  the  propellant  designer  has  to  accomplish  his 
task. 

Early  composite  solid  propellant  models  focused  on  predicting  burning 
characteristics  of  propellants  consisting  of  just  binder  and  a  single 
oxidizer.  Later,  effects  of  aluminum  in  composite  propellants  were  accounted 
for  in  combustion  models.  More  recently,  models  have  been  developed  which 
describe  high  energy  mono-propellants  and/or  explosive  burning 
characteristics.  The  goal  of  this  research  is  to  develop  a  composite  solid 
propellant  combustion  model  which  includes  ammonium  perchlorate  (AP)  crystals 
and  high  energy  cyclotetramethylene-tetranitramine  (HMX)  crystals  in  a 
polybutadiene  (HTPB)  binder.  The  model  will  be  used  to  predict  the  burning 
rate,  pressure  exponent,  and  temperature  sensitivity  for  a  variety  of 
propellants  consisting  of  various  combinations  of  HMX  and/or  AP  crystals. 
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CHAPTER  II 
BACKGROUND 


Theoretical  modeling  of  solid  propellant  combustion  serves  a  number  of 
useful  purposes  in  propellant  development.  It  represents  the  physics  and 
chemistry  of  the  combustion  processes  in  mathematical  terms  so  that  the 
combustion  characteristics  of  formulations  can  be  predicted.  It  strives  to 
explain  the  combustion  properties  which  are  observed,  providing  a  diagnostic 
tool  to  help  address  anomalous  or  undesirable  combustion  behavior.  Composite 
propellants,  in  their  simplest  form,  consist  of  a  dispersion  of  particulate 
oxidizers  within  a  matrix  of  fuel.  Even  if  all  the  particles  were  the  same 
size  and  shape  (ie.  spherical)  and  even  if  the  dispersion  of  the  oxidizer  in 
the  fuel  was  uniform,  the  combustion  of  the  propellant  would  involve  a 
multitude  of  subprocesses.  Some  of  these  processes  are  heating  of  the 
condensed  phase,  decomposition  of  the  oxidizer  and  the  fuel,  melting, 
vaporization,  mixing  in  the  vapor  phase,  and  gas-phase  combustion.  Thus,  any 
attempt  at  understanding  the  combustion  behavior  must  identify  the  key 
processes  that  control  the  burning.  Combustion  modeling  seeks  to  distinguish 
governing  mechanisms  from  those  of  secondary  of  little  importance,  so  that 
chemists  and  engineers  may  be  guided  more  productively  in  the  means  to 
overcome  combustion  limitations  and  achieve  propellant  tailoring  goals.  A 
comprehensive  model  which  describes  the  combustion  of  a  solid  propellant 
should  be  based  on  the  key  physical  and  chemical  processes  involved  in  the 
combustion  zone.  This  review  of  composite  propellant  modeling  will  review 
past  modeling  efforts.  In  addition,  the  observed  physics  and  chemistry  of 
nitramine  combustion  will  be  discussed. 


2 . 1  Non -Nitramine  Combustion  Models 


2.1.1  Granular  Diffusion  Model 

The  granular  diffusion  flame  (GDF)  model  was  formulated  to  describe  the 
burning  characteristics  of  ammonium  perchlorate  (AP)  based  composite 
propellants.  The  granular  diffusion  flame  model  was  formulated  in  1958  and 
revised  in  1960. [1,2]  Based  on  enlarged  pictures  of  the  burning  zone  of 
simple  composite  propellants,  the  authors  of  the  model  attempted  to 
theoretically  predict  the  effects  of  pressure,  oxidizer  particle  size,  and 
fuel-to-oxidizer  ratio  upon  the  mean  burning  rate.  The  pictures  clearly 
indicated  heterogeneity  both  near  the  surface  and  in  the  vapor  phase.  Instead 
of  the  popular  procedures  of  attempting  to  "average"  the  heterogeneity  in  the 
vapor  phase,  the  model  considered  the  heterogeneity  as  a  principal  feature. 

The  model  assumed  that  vapors  of  fuel  or  oxidizer  or  both  are  released  in  the 
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form  of  pockets  of  a  certain  mass  content.  The  average  mass  content  of  a 
pocket  was  assumed  to  be  much  smaller  than  that  of  an  average  oxidizer  crystal 
and  to  be  independent  of  the  temperature  and  pressure.  However,  the  size  of 
the  pockets  was  assumed  to  be  related  to  the  average  oxidizer  particle  size. 
The  rate  at  which  the  pockets  are  consumed  in  a  chemical  reaction  as  they  pass 
through  the  flame  zone  was  assumed  to  be  controlled  by  diffusional  mixing  and 
the  chemical  kinetics.  The  model  assumes  mixing  of  the  oxidizer  and  the  fuel 
occurs  only  in  the  gas  phase,  no  sub-surface  reactions,  the  heat  of  combustion 
is  released  in  a  thin  reaction  zone  in  the  gas  phase  and  a  linear  temperature 
profile  from  the  propellant  surface  temperature  to  the  temperature  of  the 
product  combustion  gases.  The  GDF  model  is  the  only  model  to  assume  a  gas 
phase  temperature  profile. 


2. 1.1.1  Equation  Development 

The  following  energy  equation  for  the  GDF  model  is  based  on  the  quasi-one- 
dimensional  model  depicted  by  Figure  1. 


Figure  1 

Theoretical  Model  of  Steady  State  Flame  for 
Granular  Diffusion  Flame  Model 
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“CS(TS-T0)  -  Ag(§i)g  +  i(Qr-Qs)  -  mQr 


(1) 


where : 

m  -  Mass  burning  rate,  ra  -  r ps,  ps  is  propellant  density 

Cs  -  Specific  heat  of  solid  (average  of  fuel  and  oxidizer) 

Ts  -  Temperature  at  surface  of  solid  (average) 

TQ  -  Initial  temperature  of  the  propellant 

X  -  Average  thermal  conductivity  of  flame  gas  at  surface 

O 

(g^)  -  Temperature  gradient  in  flame  gas  at  surface 

Qr  -  Heat  of  combustion  of  the  propellant 
Qs  -  Net  heat  release  (positive)  for  gasification 
of  the  propellant 


Using  Figure  1  and  approximating  the  temperature  gradient,  equation  (1) 
becomes : 


®[Cs(Ts-T0)-Qs] 


TfTs 


(2) 


where  L  is  the  granular  diffusion  flame  height.  The  model  assumes  that  L  can 
be  written  as  follows : 


L  -  Lj,  +  I*  (3) 

where  Lp  is  the  diffusional  mixing  distance  and  is  the  kinetic  reaction 
distance.  Since  the  mass  diffusion  coefficient  of  the  reaction  species  is 
inversely  proportional  to  pressure  while  the  chemical  reaction  rate  is 
directly  proportional  to  pressure,  the  molecular  diffusion  rate  is  much  faster 
than  the  chemical  reaction  rate  at  low  pressure.  Thus,  the  kinetic  reaction 
distance  will  be  greater  than  the  diffusional  mixing  distance  at  low  pressure, 
Lj^  »  Ljj.  At  high  pressures  the  diffusional  processes  are  slower  than  the 
kinetic  reactions  and  the  combustion  process  is  controlled  by  the  diffusional 
mixing,  Lp  »  L^. 

In  the  determination  of  Ljr  for  the  low  pressure  premixed  flame,  the  model 
treated  the  flame  as  a  stream  of  velocity  u  in  which  a  second  order  reaction 
was  taking  place.  The  mean  gas  velocity  in  the  flame  zone  is  given  by: 

u  -  m/Pg  (M 

where  pg  is  the  density  of  the  gas  in  the  flame  zone.  The  flame  zone 
thickness  for  the  premixed  flame  is  given  by: 

L  -  ur  (5) 

where  r  is  the  time  it  takes  the  chemical  reaction  to  occur.  Assuming  a 
global  second  order  reaction  for  the  premixed  gases,  the  reaction  time  is 
given  by: 

r  -  [(1  -  £)2PgAexp(-E/RTg)]  (6) 

where  an  Arrhenius  expression  is  assumed  for  the  rate  of  reaction,  Tg  is  the 
average  gas  phase  reaction  zone  temperature  and  e  is  a  function  of  the 
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products  of  reaction,  i.e.,  «  is  equal  to  zero  at  the  propellant  surface  and  e 
is  equal  to  unity  at  the  completion  of  the  reaction.  For  simplicity,  the 
model  assumed  that  the  average  value  for  (1-e)^  over  the  flame  zone  is  unity. 
Physically  this  means  that  most  of  the  overall  reaction  takes  place  near  the 
propellant  surface. 


Combining  equation  (4),  (5)  and  (6),  yields  the  following  for  the  flame 
zone  thickness  for  the  low  pressure  premixed  flame: 


Lr  “  ™/ ( PgAexp ( - E/RTg) ] 


(7) 


when  equation  (2)  is  solved  for  the  total  mass  flux  and  the  result  substituted 
into  equation  (7),  the  following  thickness  is  obtained: 


% 


VTf  ‘  TS> 


LtCs(Ts  -  T0)  -  Qs ] J 


1/2 


PgtAexpC-E/RTg)]1/2 


(8) 


For  the  high  pressure  diffusion  flame,  the  model  assumed  that  the  mass  of 
a  pocket  could  be  expressed  as : 


(9) 


where  D  is  the  oxidizer  particle  diameter.  The  lifetime  of  a  pocket  is 
determined  by  the  rate  of  gas  diffusion  to  the  surrounding  flame.  The 
molecular  diffusivities  of  the  fuel  and  the  oxidizer  are  assumed  to  be  the 
same  and  are  averaged  over  the  gaseous  reaction  zone.  Under  these  conditions 
the  lifetime  of  a  fuel  pocket  is: 


D 


2 


T 


(10) 


where  Dg  is  the  average  diffusivity  of  the  oxidizer  and  fuel  vapors.  The 
thickness  of  the  flame  zone  is  the  product  of  the  velocity  of  the  gas  stream 
and  the  lifetime  of  the  fuel  pocket.  Thus: 


lD  “  ^T 


m  D 


pg  Dg 


(ID 


Combining  equation  (9)  and  (11)  yields  the  diffusion  flame  zone  thickness: 


k>  " 


,2/3m 


'i/3>6 


(12) 


Solving  equation  (3)  for  the  total  mass  flux  and  substituting  this  into 
equation  (12)  yields: 


h)  “ 


VTf  •  V 


L(Cs(Ts  -  T0)  -  QS]J 


1/2  ^1/3 


'T'Y2 


(13) 
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Expressions  for  the  flame  zone  thickness  in  the  two  extreme  cases  can  now 
be  written.  For  the  general  intermediate  case,  the  model  assumed  that  the 
flame  zone  thickness  varies  with  pressure  partly  as  if  it  were  reaction  rate 
controlled,  Lj^,  and  partly  as  if  it  were  diffusion  controlled,  LD,  and  that 
the  flame  zone  thickness  could  be  expressed  as  the  sum  of  the  above  two. 

Thus,  the  general  expression  for  the  flame  zone  thickness,  from  equation  (3), 
becomes : 


Ag(Tf  -  Ts).- 

1/2  (  1  PV3 

.[CS(TS  t0)-qs3. 

{pg[Aexp(-E/RTg)]1/2  p5/6Dl/2 

(14) 


Using  the  ideal  gas  law,  the  expression  D  -  K^Tg/2/p  for  the  diffusivity,  an 
average  gas  phase  temperature  and  letting  m  -  r ps  where  ps  is  the  density  of 
the  solid  propellant,  an  expression  for  the  burning  rate  can  be  obtained  by 
substituting  equation  (14)  for  L  into  the  energy  equation,  equation  (3): 


1  a 


r  P 


(15) 


where  the  constants  a  and  b  are  given  by: 

Ag(Tf  -  Ts)  U/2 


a  — 


b  - 


[CS(TS-T0)-QS] 

VTf  *  v 


's“g 


Lfcs<Ts-To>-QslJ 


[Aexp(-E/RTg)]1/2 

1/2  R5/6X  1/12  1/3 

a _ s _ 

kJ/2 


(16) 


(17) 


The  constants  a  and  b  are  independent  of  pressure,  p,  and  are  related  to 
and  Ljj.  If  »  Lq  (reaction  rate  controlled)  a  »  b  and  the  burning  rate 
exponent,  n,  in  the  burning  rate  expression  that  follows: 


r  -  cp 


n 


(18) 


approaches  1.  If  Ljj  »  Lj^  (diffusional  mixing  controlled)  b  »  a  and  the 
burning  rate  exponent,  n,  approaches  1/3.  If  Lp  is  of  the  same  order  as  LK, 
then  n  is  between  1  and  1/3.  Experimentally,  with  AP/inert  binder 
propellants,  these  results  are  obtained. 


2. 1.1. 2  Discussion 

The  granular  diffusion  flame  model  used  equation  (15)  to  correlate 
experimental  burning  rate  data  for  a  range  of  ammonium  perchlorate  propellant 
formulations.  The  equations  provided  good  agreement  between  theory  and 
experiment  in  predicting  the  effect  of  oxidizer  crystal  size,  fuel  to  oxidizer 
ratio,  and  pressure  on  burning  rate.  The  GDF  model  developed  a  burning  rate 
expression  that  provided,  in  a  closed  form,  a  means  of  estimating  the  effects 
of  varying  the  important  propellant  combustion  parameters  on  burning  rate. 
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The  model  was  later  modified  to  take  into  account  a  distended  AP 
monopropellant  flame.  This  was  done  to  make  the  model  applicable  to  low 
pressure  combustion. [3]  The  model  was  also  modified  to  model  polydisperse 
oxidizer  distributions  by  applying  a  statistical  technique  to  the  differing 
particle  sizes. [4] 


2.1.2  Hermance  Heterogeneous  Reaction  Model 

The  Hermance  heterogeneous  reaction  (HR)  model  was  based  on  a  detailed 
combination  of  the  steady  state  decomposition  processes  of  the  fuel  and 
oxidizer.  The  model  was  first  formulated  in  1966  and  later  revised  in 
1967. [2,5,6]  The  model  postulates  an  exothermic,  heterogeneous  reaction 
between  the  solid  fuel  and  oxidizer  decomposition  products  and  assumes  that 
the  gas  phase  flame  position  is  the  sum  of  the  lengths  associated  with 
diffusional  mixing  and  chemical  reaction.  The  heterogeneous  reaction  model 
was  a  significant  advance  in  the  state  of  the  art  of  combustion  modeling  since 
it  was  the  first  model  to  deal  directly  with  the  heterogeneity  of  the  solid 
phase  by  taking  into  consideration  arbitrary  particle  sizes. 

The  model  assumes  that  the  steady  state  combustion  of  a  composite 
propellant  occurs  in  three  regions;  the  propellant  surface,  the  gas  phase 
flame  zone  and  the  region  between  the  surface  and  the  flame.  The  surface 
chemical  processes  are  endothermic  fuel  pyrolysis,  exothermic  oxidizer 
decomposition  and  an  exothermic  heterogeneous  chemical  reaction  between  the 
fuel  binder  and  decomposed  oxidizer  in  small  regions  surrounding  individual 
oxidizer  crystals.  Each  of  these  reactions  produces  a  mass  flux  from  the 
propellant  surface. 

In  the  gas  phase  flame  zone,  an  exothermic  chemical  reaction  occurs 
between  the  oxygen  rich  products  from  the  decomposed  oxidizer  and  the 
paralyzed  fuel.  Heat  from  the  flame  zone  is  fed  back  to  the  propellant 
surface  by  conduction.  In  the  region  between  the  surface  and  the  flame  zone, 
the  surface  decomposition  products  mix  by  diffusion  and  undergo  a  chemical 
kinetic  delay  before  ignition  at  the  flame.  In  the  original  model  only  the 
ignition  delay  was  considered. [5]  The  processes  depend  on  the  pressure  level 
and/or  the  temperature  at  the  location  of  the  process  and  are  linked  together 
by  the  temperature  distribution  in  the  gas  and  solid  phases. 


2 . 1 . 2 . 1  Equation  Development 

One  of  the  more  important  aspects  of  the  model  is  the  calculation  of  the 
intersection  diameter  of  an  oxidizer  crystal  with  the  burning  surface.  If  the 
propellant  burning  surface  is  visualized  as  a  plane  passing  through  a  randomly 
packed  bed  of  spherical  oxidizer  crystals,  it  can  be  shown  that  the 
statistical  average  intersection  diameter  is: 

D'  -  fT/?  Dc  (19) 

where  DQ  is  the  actual  crystal  diameter  and  D'  is  the  statistical  intersection 
diameter. [7] 
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The  propellant  burning  rate  can  be  determined  from  the  total  mass  flux 
issuing  from  the  propellant  surface.  It  is  the  sum  of  the  mass  fluxes  of 
gaseous  species  produced  by  each  of  the  surface  decomposition  processes 
multiplied  by  the  fraction  of  the  total  propellant  surface  area  associated 
with  each  of  the  processes.  Thus: 

r  -  ^[mf(Sf/S0)  +  m0X(S0X/S0)  +  msr(Ssr/SQ) ]  (20) 

where  p is  the  density  of  the  propellant,  m  is  the  mass  flux  and  S  is  surface 
area.  rThe  subscripts  f,  ox,  o  and  sr  designate  fuel,  oxidizer,  total 
propellant  and  surface  reaction  respectively.  Assuming  a  planar  surface,  the 
ratio  of  fuel  and  oxidizer  surface  area  to  the  total  propellant  surface  area 
can  be  expressed  as  a  function  of  the  volume  fraction  of  oxidizer,  f ,  in  the 
propellant: 


s  S 

“  (1  -  f )  and  — -  f 
-’o  °o 


(21)  and  (22) 


The  area  on  which  the  surface  reaction  occurs  is  calculated  by  postulating 
that  an  oxidizer  crystal  decreases  in  size  during  decomposition  as  shown  in 
the  top  of  Figure  2.  This  size  reduction  produces  a  fissure  of  depth  e 
between  the  oxidizer  crystal  and  the  fuel  binder  which  is  the  region  where  the 
heterogeneous  reaction  takes  place.  Arrhenius  type  functions  are  used  to 
describe  the  surface  regression  of  both  fuel  and  oxidizer. 


The  model  divides  the  steady  state  combustion  process  of  composite 
propellants  into  three  regions  as  shown  in  the  bottom  of  Figure  2 .  The 
differential  form  of  the  energy  equation  is  solved  for  each  region 
simultaneously.  The  required  boundary  conditions  result  by  matching 
temperature  and  heat  flux  at  the  interface  of  each  section.  The  analysis 
results  in  three  equations  with  three  unknowns.  They  are  the  burning  rate 
equation,  the  surface  temperature  equation  and  the  flame  temperature  equation. 
These  equations  in  condensed  form  are: 


bi  Pc 


-  aiexp(-Ef/E*s)  +  —— 

ft1/*- 


c3 

dl 

exp 

'-<Eox  +  Esr>‘ 

pin 

r 

E*s  - 

(23) 


bi  Pfl 


b2  + 


>V2 


[c3  dl 

exp 

'-<Eox  +  Esr>‘ 

i  a 

I  Cl, 

* 

u 

- J 

E*s  - 

a 

^exp(-Ef/E0)  +  a3exp(-£*) 


(24) 


-  6S  +  a3(l  -  exp(-£*)  (25) 

where  a^ ,  a3«  &!»  t>2  ’  ^3  ancl  dl  are  constants  containing  the  physical  and 

chemical  properties  of  the  propellants;  E  is  the  activation  energy,  £  is  the 
non-dimensional  flame  standoff  distance  and  P  is  the  ambient  pressure.  The 
burning  rate,  surface  temperature  and  flame  temperature  are  calculated  by 
numerical  iteration  using  the  above  three  equations. 


21 


NWC  TP  6992 


Figure  2 

Physical  Description  of  the  Hermance  Heterogeneous 
Reaction  Model 
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2. 1.2. 2  Discussion 

The  Hermance  heterogeneous  reaction  model  solved  the  above  three 
simultaneous  equations  with  a  set  of  input  data  which  represented  the 
combustion  properties  of  a  polysulfide  AP  composite  propellant.  The  model 
obtained  a  reasonable  fit  to  experimental  burning  rate  data  taken  over  the 
pressure  range  from  1  to  400  atmospheres.  The  model  qualitatively  predicts 
the  effect  of  oxidizer  particle  size  on  the  propellant  burning  rate  at 
pressures  above  200  atmospheres. 

At  low  pressure  the  model  predicts  very  little  heterogeneous  reaction 
while  at  high  pressure  the  crevice  becomes  huge  and  the  heterogeneous  reaction 
dominates,  as  depicted  in  Figure  3.  Moreover,  the  calculated  surface 
temperature  is  almost  constant  with  increasing  pressure  indicating  that  the 
temperature  dependent  Arrhenius  type  function  will  be  essentially  constant. 
Thus,  the  burning  rate  characteristics  predicted  by  the  HR  model  are 
determined  almost  completely  by  the  crevice  and  its  formation.  Yet  there  is 
no  physical  evidence  to  support  the  formation  of  a  crevice  as  postulated  by 
the  model . [ 5 ] 


2.1.3  Beckstead,  Derr,  Price  (BDP)  Model 

The  goal  of  the  Beckstead,  Derr,  Price  model  was  to  develop  a  new 
analytical  model  of  composite  propellant  combustion  based  on  a  realistic 
physical  model,  utilizing  information  derived  both  from  experimental  studies 
concerning  the  surface  structure  of  extinguish  .d  propellants  and  previous 
analytical  modeling  ef forts . [ 1-6 , 8 - 12 ]  The  goal  was  not  to  predict  exact 
burning  rates,  but  of  predicting  changes  in  burning  rate  for  a  given  change  in 
propellant  formulation. 

After  extensive  experimental  observation  of  the  surface  structure  of  a 
burning  composite  solid  propellant,  Beckstead,  Derr  and  Price  formulated  a 
multiple  flame  composite  propellant  combustion  model  based  on  the  postulated 
existence  of  three  flames  in  the  region  of  each  exposed  oxidizer 
particle .[ 13 , 14]  The  three  flames  are  the  AP  flame,  the  primary  flame  and  the 
final  diffusion  flame. 

The  AP  flame  occurs  between  the  decomposition .products  of  the  AP  forming 
C>2  and  inert  products.  Thus: 

NH-j  +  HCLO^  - *  Inert  Products  +  C>2  (26) 

Since  the  AP  decomposition  products  are  premixed,  the  AP  flame  standoff 
distance  is  dependent  only  upon  chemical  kinetics. 

The  primary  flame  occurs  between  fuel  pyrolysis  products  and  the  AP 
decomposition  products.  Thus: 

Fuel  Products  +  HClO^  (decomposition  products)  - * 

Primary  Flame  Combustion  Products  (27) 

The  primary  flame  height  is  dependent  on  diffusional  mixing  as  well  as 
kinetics . 
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LOW  PRESSURE -VERY  SMALL  CREVICE  AREA 


HIGH  PRESSURE-VERY  LARGE  CREVICE  AREA 


Figure  3 

The  Pressure  Dependence  of  the  Hermance  Crevice 
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The  final  diffusion  flame  occurs  between  the  binder  decomposition  products 
and  the  oxidizing  products  from  the  AP  flame.  Since  the  reactants  of  the 
final  diffusion  flame  are  preheated,  the  kinetics  are  very  fast,  thus,  the 
final  flame  height  is  dependent  only  on  diffusional  mixing.  The  geometric 
relationship  of  the  assumed  flame  structure  is  shown  in  Figure  4. 

The  flame  standoff  distances  are  a  function  of  the  combustion  pressure. 

At  low  pressures  (<100  atm)  the  propellant  is  considered  to  burn  as  a  premixed 
flame  with  the  oxidizer  and  binder  decomposition  products  mixing  completely 
before  a  reaction  occurs.  As  the  pressure  increases,  the  mixing  path  length 
increases,  and  the  reaction  path  length  decreases  so  that  the  two  reaction 
paths  (ammonia  +  oxidizing  products  and  binder  products  +  oxidizing  products) 
become  competitive.  At  higher  pressures  (>100  atm)  the  ammonia  reacts  with 
the  oxidizing  products  before  the  binder  products  can  diffuse  into  the 
oxidizer  stream  and  react.  Due  to  the  fuel  rich  nature  of  composite  solid 
propellants,  a  final  diffusion  flame  always  occurs  above  the  AP  flame.  The 
flame  structure  variation  with  pressure  is  depicted  in  Figure  5. 


Figure  4 

BDP  Multiple  Flame  Structure 
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LOW  PRESSURE  (flF>|0) 


Figure  5 

Multiple  Flame  Structure  Variation  with  Pressure 
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At  the  propellant  surface  the  initial  decomposition  step  of  both  the 
binder  and  the  oxidizer  is  endothermic.  However,  while  still  adsorbed  on  the 
surface,  either  set  of  products  can  undergo  a  condensed  phase  reaction  before 
passing  into  the  gas  phase.  Products  from  the  surface  decompositions  then 
pass  into  the  gas  phase  and  begin  mixing  and  reacting. 

The  burning  surface  of  the  exposed  oxidizer  particle  is  assumed  to  be 
spherical  in  shape.  Due  to  an  ignition  delay,  the  burning  oxidizer  surface 
protrudes  above  the  planar  fuel  surface  at  low  pressure  but  is  recessed  below 
the  planar  fuel  surface  at  high  pressure.  The  oxidizer  regression  is  assumed 
to  be  the  overall  rate  controlling  mechanism,  and  an  overall  average 
temperature  is  defined  for  the  entire  burning  surface. 

In  addition  to  the  above  assumptions,  radiation  heat  transfer  is 
neglected,  an  average  value  for  the  solid  and  the  gas  phase  specific  heat  is 
assumed,  and  the  mass  diffusion  coefficient  and  thermal  conductivity  of  the 
reacting  gases  is  averaged  over  the  gas  phase  reaction  zone. 


2. 1.3.1  Conservation  of  Mass 


In  the  model's  original  form,  the  oxidizer  was  assumed  to  be  monomodal, 
monodisperse  and  spherical.  However,  the  model  has  since  been  extended  to 
bimodal  and  trimodal  monodisperse  propellants . [15-18]  In  the  original  BDP 
model  the  propellant  surface  temperatures  of  the  AP  and  binder  were  assumed 
equal.  Later  versions  of  the  model  accounted  for  differing  surface 
temperatures. [19-20]  The  case  of  monomodal  monodisperse  particle  distribution 
and  an  uniform  surface  temperature  will  be  discussed  here . 

It  is  assumed  that  global  kinetics  adequately  describe  the  decomposition 
of  the  oxidizer  and  binder.  Once  a  surface  temperature  is  assumed,  the  mass 
flux  rates  of  both  the  binder  and  AP  or  oxidizer  are  calculated  from  the 
following  Arrhenius  expressions: 

®ox  “  rox^ox  “  ^oxexP(*^ox/^s)  (28) 

mf  -  rfpf  -  Afexp(-Ef/RTS)  (29) 

where  A  is  the  Arrhenius  frequency  factor,  E  is  the  activation  energy,  R  is 
the  gas  constant  and  Ts  is  the  average  surface  temperature.  The  subscripts  ox 
and  f  represent  the  oxidizer  and  fuel,  respectively. 


The  conservation  of  mass  equation  can  be  written  as  follows: 


ff 

fs 

box 

-So- 

+  ™ox 

So  - 

(30) 


This  represents  the  total  mass  flux  of  the  burning  propellant.  SQ  is  the 
total  propellant  area  and  r  is  the  total  burn  rate.  The  propellant  density  is 
Pp.  The  oxidizer  mass  fraction  is  defined  as: 
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(31) 


where  a  is  the  oxidizer  mass  fraction.  Rearranging  Equations  (30)  and  (31)  a 
new  expression  for  the  total  mass  flux  can  be  derived: 


(32) 


In  these  expressions  everything  is  known  except  the  area  ratios.  The 
propellant  surface  geometry  will  be  discussed  next. 


2. 1.3. 2  Surface  Geometry 


The  goal  of  this  section  is  to  find  an  expression  for  the  area  ratio: 


S 


ox 


S 


o 


This  ratio  is  the  oxidizer  particle  surface  area  ratio  to  the  area  of  the 
total  propellant.  Because  of  the  random  nature  of  the  oxidizer  particles  a 
statistical  description  of  the  propellant  surface  must  be  used.  This  was 
first  used  in  the  fore  mentioned  heterogeneous  reaction  model. [5]  If  the 
propellant  burning  surface  is  visualized  as  a  plane  passing  through  a  randomly 
packed  bed  of  spherical  oxidizer  crystals,  it  can  be  shown  that  the 
statistical  average  intersection  diameter  is: [7] 

D'  -  JTJl  D0  (19) 


where  DQ  is  the  actual  oxidizer  crystal  diameter. 


Referring  to  Figure  6,  the  exposed  oxidizer  surface  area  to  the  total 
propellant  area  ratio  can  be  determined  solely  on  geometrical  principles. 
This  ratio  can  be  expressed  as: 
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ox 
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o 


(33) 


where  the  term  h/DQ  is  given  by: 
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Figure  6 

Geometrical  Relationship  of  the  Oxidizer  Crystal 
to  the  Burning  Area 
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-  -  £(1  ±  7173) 

Do 

In  the  above  expressions  f  is  the  volume  fraction  of  oxidizer,  h  is  the 
distance  a  oxidizer  crystal  either  protrudes  above  or  is  depressed  below  the 
surface  of  the  propellant,  t^gn  is  the  ignition  delay  time  of  the  oxidizer 
crystal  and  DQ  is  the  oxidizer  crystal  diameter.  rQX  and  r£  are  the  burning 
rates  of  the  oxidizer  and  fuel  respectively,  given  by  equations  (28)  and  (29). 
The  +  or  -  in  equation  (34)  refers  to  the  positive  or  negative  particles,  ie. 
whether  the  lower  hemisphere  or  the  upper  hemisphere  is  cut  by  the  planar 
surface.  It  is  assumed  that  there  are  an  equal  number  of  positive  and 
negative  crystals. 


1  - 


ox 


r<r 


+  r 


-ign 


ox 


(34) 


2 . 1 . 3 . 3  Energy  Balance 

Figure  7  shows  the  heat  transfer  from  the  three  flames  above  the 
propellant's  surface,  the  surface  decomposition  energies  and  the  state  values 
of  enthalpy  at  the  propellant  surface  and  deep  into  the  propellant.  From 
Figure  7,  the  one -dimensional  energy  balance  at  the  burning  surface  is: 

mrCp(Ts-To)  -  -<5>0X(S0X/So)Ql  -  mT(Sf/S0)Qf  + 

^fQpF™TexP^'^PF  +  (l'/?f  )™ox^ox/^o^^APex^^ "£Ap)  +  QFFex£^  "^FF^  (35) 

where  Cp  is  the  average  specific  heat  of  the  gases  and  the  solid  and  TQ  is  the 
initial  propellant  temperature,  QL  and  Qf  represent  the  energy  required  to 
vaporize  the  oxidizer  and  fuel  respectively,  QpF,  QAp,  and  QFF  represent  the 
energy  released  in  the  final  flame,  and  AP  flame  and  the  primary  flame 
respectively,  /?£  is  the  fraction  of  reactants  that  react  in  the  primary  flame, 
and  CpF>  ^AP’  aTU*  ^FF  are  t*ie  non-dimensional  flame  standoff  distances  for  the 
final  flame,  the  AP  flame  and  the  primary  flame,  respectively. 

After  rearranging  equation  (35)  and  combining  with  equation  (31),  the 
equation  for  the  average  surface  temperature  can  be  written: 


T  -  T 

4s  lo 

i 

1 

/-S 

-a)Qp/Cp  +  Pf (Qpp/^p)exP('^PF^ 

+  (l-/9f)((QAP)/Cp)exp(-e*p) 

+  (Qff/Cp)6Xp(*^FF^ 

(36) 

The  heat  release  terms  can 

be  expressed  by: 

Qpp  “ 

Cp(Tf  -  Tq) 

+  qQl  +  (l-a)Qf 

(37) 

QAP  ” 

Cp(TAP  *  T0) 

+  Ql 

(38) 

Qff  - 

(Cp/a) [ (Tf  - 

Tq)  -  o(TAp  -  Tc)  +  (1-cOQf/Cp] 

(39) 

where  Tf  and  TAp  are  the  final  flame  temperature  and  the  AP  monopropellant 
flame  temperature  respectively. 
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Heat  Transfer  from  Flames  to  Propellant  Surface 


-  Heat  Transfer  from  Surface  Decomposition  Reactions 


-  The  Enthalpy  Terms  at  the  Surface  and  Deep  into  the  Propellant 


Figure  7 

The  BDP  Surface  Energy  Balance 
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2. 1.3. 4  Flame  Standoff  Distances 

The  non-dimensional  primary  flame  standoff  distance  is  the  sum  of  the 
diffusional  mixing  distance  plus  the  kinetic  reaction  distance  and  is  given 
by: 


+  x*  ) 
PF 


c_mT 

pj--* 


-x  + 
PD 


cpmT2 


(40) 


In  this  expression  Xpp  is  the  effective  diffusional  standoff  distance,  Xpp  is 
the  kinetic  reaction  distance,  A  is  the  thermal  conductivity,  kpF  is  the 
reaction  rate  constant  for  the  primary  flame  and  6  is  the  reaction  order.  The 
non-dimensional  AP  monopropellant  flame  standoff  distance  consists  of  a 
kinetic  reaction  distance  only  and  is: 
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(41) 


The  final  flame  standoff  distance  is  assumed  purely  a  diffusional  flame  but  is 
assumed  to  occur  above  the  AP  monopropellant  flame.  Hence,  is  the  sum  of  the 
AP  kinetic  flame  standoff  distance  plus  the  diffusional  mixing  distance  of  the 
final  flame  and  is  given  by: 


cpmox  *  _* 


^FF  “  X  <XAP  +  XFD>  “ 


cpmox 
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Vox  _* 

X  XFD 


(42) 


Using  global  kinetics  the  generic  form  of  the  kinetic  reaction  distance 
can  be  expressed  as : 


*  Cpm  * 

s - v— x 


(43) 


To  determine  the  diffusional  flame  standoff  distance  the  authors  assumed 
that  the  burning  oxidizer  crystal  surrounded  by  fuel  could  be  represented  by 
Bunsen  burner  type  of  flame.  An  analysis  reported  by  Burke  and  Schumann  on 
this  type  of  flame  and  later  modified  by  Williams,  was  used  to  arrive  at  the 
following  equation  for  the  average  diffusional  flame  standoff 
distance . [13,21,22] 
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(44) 


In  this  equation  C2  is  a  constant  related  to  the  propellant  properties,  7  is  a 
constant  relating  to  the  stoichiometry  of  the  flame,  is  the  first  non  zero 
root  of  the  first  order  Bessel  function  and  b  is  a  characteristic  dimension  at 
the  surface.  The  term  A^  is  a  constant  that  determines  the  effective  planar 
flame  height.  The  Burke -Schumann  flame  height  analysis  determines  the  maximum 
height  of  the  diffusion  flame.  The  diffusional  heights  used  in  equation  (40) 
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and  (42)  are  effective  planar  flame  heights.  In  the  equations  (43)  and  (44) 
the  appropriate  mass  flux,  reaction  order,  rate  constant,  and  physical 
constants  are  used  to  arrive  at  correct  flame  standoff  distances  used  in 
equations  (40) ,  (41)  and  (42) . 


2. 1.3. 5  The  Competing  Flames 


There  is  one  final  term  needed  in  the  energy  equation  which  has  not  been 
defined.  This  term,  /?£,  is  the  fraction  of  the  total  reactants  issuing  from 
the  propellant  surface  that  react  in  the  primary  flame.  This  term  is  required 
since  there  is  competition  for  the  oxidizer  decomposition  products  between  the 
primary  flame  and  the  AP  monopropellant  flame.  Based  solely  on  geometry,  and 
assuming  that  the  primary  diffusion  flame  is  parabolic  in  shape,  this  term  can 
be  calculated  with  the  following  expression: 


(45) 


In  this  expression  the  numerator  is  the  difference  between  the  primary  flame 
and  AP  monopropellant  flame  kinetic  reaction  distance.  The  denominator  is  the 
actual  primary  diffusion  flame  height  given  by  equation  (44)  with  A^  equal  to 
one . 


2. 1.3. 6  Discussion 

Comparisons  were  made  with  a  series  of  unimodal  AP  polysulfide 
propellants . [13]  The  predicted  dependence  of  burning  rate  on  oxidizer 
concentration  was  nearly  the  same  as  observed  experimentally,  where  the 
dependence  of  the  burning  rate  on  the  oxidizer  particle  size  was  greater  than 
observed  experimentally.  The  dependence  of  the  burning  rate  on  the  initial 
temperature  was  found  to  be  In  qualitative  agreement  with  the  available 
experimental  data. 

The  Beckstead,  Derr,  Price  multiple  flame  composite  propellant  combustion 
model  laid  the  foundation  for  the  subsequent  Petite  Ensemble  Model  which  will 
be  described  next. 


2.1.4  Petite  Ensemble  Model  (PEM) 

In  both  the  Beckstead,  Derr,  Price  multiple  flame  model  and  the  Hermance 
heterogeneous  reaction  models,  it  was  assumed  that  the  macroscopic  behavior  of 
an  ensemble  of  different  flames  surrounding  several  sized  oxidizer  particles 
can  be  replaced  by  a  multiple  flame  structure  centered  about  a  single, 
characteristically  sized  oxidizer  particle.  In  actuality,  the  burning 
propellant  surface  is  much  more  complicated  with  many  different  sized  oxidizer 
crystals  all  burning  simultaneously  and  reacting  with  a  binder  medium.  A 
statistical  approach  to  deal  with  a  more  complicated  burning  surface  structure 
was  first  presented  in  1974.(23]  The  combination  of  a  statistical  formalism 
within  a  Beckstead,  Derr,  Price  multiple  flame  model  was  called  the  Petite 
Ensemble  Model  (PEM) . [24-28] 
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In  the  PEM,  the  burning  surface  of  a  composite  propellant  is  assumed  to 
consist  of  a  random  arrangement  of  polydisperse ,  or  multi-sized,  oxidizer 
particle/fuel  binder  surface  pairs.  It  is  further  assumed  that  each  oxidizer 
particle  at  the  propellant  surface  has  associated  with  it  a  portion  of  the 
total  available  fuel  binder  that  is  assumed  to  surround  every  oxidizer 
crystal.  Therefore,  each  oxidizer  particle/fuel  binder  pair  will  produce  a 
unit  flame  structure  centered  about  the  oxidizer  particle.  If  it  is  assumed 
that  all  of  the  oxidizer/fuel  pairs  burn  independently  of  each  other,  then  the 
burning  propellant  surface  can  be  mathematically  rearranged  into  imaginary 
families  of  monodisperse ,  or  all  one  size  oxidizer  particles,  propellants. 
These  monodisperse  propellants  are  subsequently  referred  to  as  pseudo 
propellants . 

The  statistically-based  PEM  calculates  the  overall  burning  rate  of  a 
polydisperse  propellant  with  a  summation  scheme  incorporating  the  calculated 
pseudo  propellant  burning  rates.  That  is,  if  the  available  combustion  model 
can  determine  the  burning  rate  of  each  of  the  monodisperse  pseudo  propellants, 
or  more  precisely,  of  each  of  the  oxidizer/fuel  pairs,  then  with  a  knowledge 
of  the  oxidizer  particle  size  distribution,  the  overall  propellant  burning 
rate  can  be  determined. 


2. 1.4.1  Statistical  Formulation 

Before  presenting  the  expression  describing  the  overall  propellant  burning 
rate  as  a  function  of  the  pseudo  propellant  burning  rates,  a  description  of 
the  oxidizer  size  and  size  distribution  will  be  presented.  A  composite 
propellant  is  a  heterogeneous  mixture  of  finely  ground  oxidizer  particles  all 
mixed  in  with  a  liquid  fuel  binder  which  is  polymerized  into  the  solid 
propellant.  The  oxidizer  particles  used  in  composite  propellants  range  in 
size  from  less  than  a  micron  to  several  hundred  microns  in  diameter. 

The  number  density  of  oxidizer  particles  in  a  composite  propellant  usually 
varies  in  a  log  normal  distribution  fashion  with  the  size  of  the  particle. 

This  variation  of  particle  number  density  with  particle  diameter  is  termed  a 
particle  size  distribution.  Figure  8  depicts  a  typical  particle  size 
distribution  for  a  nominal  diameter  of  10  microns.  Composite  propellants 
usually  contain  several  distribution  modes  in  order  to  achieve  high  solids 
loading  to  obtain  high  oxidizer  to  fuel  ratios.  A  log  normal  distribution 
function  can  be  expressed  by  the  following  equation: 


1 

-1 

■^n®ox*^n^ox 
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The  mean  diameter,  Dox  is  the  fifty  percent  weight  mean  diameter,  since  half 
of  the  oxidizer's  mass  is  composed  of  diameters  less  than  DQX.  The  term  aQX 
represents  the  actual  width  of  the  distribution  about  this  mean  oxidizer 
diameter.  A  mode  width  parameter,  oQX,  equal  to  unity  corresponds  to  a 
monodisperse  oxidizer  mode;  that  is  all  oxidizer  particles  in  that  mode  are  of 
one  diameter,  BQX.  A  mode  width  parameter  exceeding  unity  represents  a 
polydisperse  oxidizer  mode  distribution.  In  a  polydisperse  mode,  the  oxidizer 
mass  is  distributed  about  the  mean  particle  size;  the  larger  the  value  of  oox, 
the  wider  the  distribution  of  oxidizer  mass  about  Dox. 
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By  utilizing  the  log  normal  distribution  function  given  in  equation  (46), 
each  oxidizer  mode  within  a  composite  propellant  can  be  represented  by  the  two 
independent  parameters,  EQX  and  <7QX,  and  the  oxidizer  particle  size 
distribution  can  easily  be  characterized  by  a  simple  log  normal  expression  of 
the  form  given  by  equation  (46) .  An  expression  for  the  overall  polydisperse 
propellant  burning  rate,  r,  can  now  be  presented: 


r 


-  t(rdFd>dInDox 


D, 


ox 


(47) 


In  this  expression,  r^  and  are  the  burning  rate  and  distribution  function, 
respectively ,  of  each  of  the  monodisperse  pseudo  propellants  of  a  given 
oxidizer  diameter  (as  designated  by  the  subscript  d) . 


2. 1.4. 2  Separate  Surface  Temperatures 

In  the  development  of  the  basic  equations  comprising  the  Petite  Ensemble 
Model  (PEM) ,  many  of  the  assumptions  evoked  in  the  Beckstead,  Derr,  and  Price 
(BDP)  model  discussed  previously  are  retained.  This  is  especially  true  in  the 
description  of  the  surface  geometry  of  the  oxidizer  crystal  embedded  in  the 
fuel  binder.  In  the  subsequent  analysis,  it  will  be  assumed  that  the 
propellant  is  composed  of  spherical  oxidizer  crystals  and  a  fuel  binder. 

In  describing  the  propellant  surface,  it  is  assumed  that  the  oxidizer 
crystals  are  randomly  mixed  within  a  fuel  binder.  Like  the  Beckstead, 
Derr,  Price  multiple  flame  model  the  statistical 
average  intersection  diameter,  D' ,  is: 

D'  -  7275  D0  (19) 

As  in  the  BDP  model  the  oxidizer  mass  burning  rate  and  the  fuel  binder  burning 
rate  are  assumed  to  be  represented  adequately  by  an  Arrhenius  expression  based 
on  the  surface  temperature.  Unlike  the  original  BDP  model,  the  surface 
temperatures  of  the  oxidizer  and  binder  are  different.  The  Arrhenius 
expressions  are  as  follows: 

“ox  ~  roxpox  "  AoxexP(*Eox/^Ts ,ox) 

mf  -  r fpf  -  Afexp(-Ef/RTsf)  (49) 

In  these  expressions  Tg  Qx  and  Tg  £  are  the  surface  temperature  of  the 
oxidizer  and  fuel,  respectively. 

2. 1.4. 3  Surface  Energy  Balance 

Like  the  BDP  multiple  flame  model  the  entire  PEM  analysis  of  the  burning 
behavior  of  a  solid  propellant  rests  upon  the  solution  of  the  energy  balance 
equation  written  at  the  propellant  surface.  The  energy  transmitted  back  to 
the  surf*«,e  from  the  individual  flames  situated  above  the  surface  must  be 
equal  to  the  heat  required  to  bring  the  propellant  from  an  initial  propellant 
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temperature,  TQ,  up  to  a  given  temperature  at  the  surface  and  to  supply  heat, 
if  required,  for  constituent  gasification  at  that  surface.  It  has  already 
been  shown  that  the  rate  of  gasification  of  each  of  the  reactants  is  tied 
directly  to  its  respective  surface  temperature.  The  surface  energy  equation 
couples  the  solid  phase  and  gas  phase  processes  occurring  simultaneously 
during  the  burning  of  solid  propellant. 

Assuming  the  oxidizer  has  a  surface  temperature,  Tg  ,  and  the  fuel 
binder  has  a  possible  different  surface  temperature,  Ts’b,  the 
surface  energy  balance  takes  the  following  form: 

aCP.°X<TS,0X-T0>  +  a-«>S,f(TS,b-T0>  “  QQL  +  (l'a)Qf 

+  ^fQppexp(-CpF)  +  (l-y3f)QQApexp(-^p)  +  (1-0  f  )QFFexp( -£pF)  (50) 

The  terms  on  the  left-hand  side  of  equation  (50)  are  simply  the  sensible  heat 
terms  associated  with  bringing  the  oxidizer  and  fuel  binder,  in  their 
respective  relative  proportions,  up  to  their  surface  temperatures.  It  is 
assumed  throughout  this  analysis  that  the  oxidizer  and  fuel  binder  have 
different  values  of  solid  phase  specific  heat,  C  ox  and  Cp  p,  respectively. 
The  next  two  terms  represent  net  energy  either  absorbed  or  released  at  the 
surface,  Qb  is  the  latent  heat  of  oxidizer  surface  decomposition  and  Qp  is  the 
heat  of  pyrolysis  for  the  fuel  binder.  The  convention  adopted  for  this 
analyses  is  that  the  terms  QL  and  Qp  are  positive  if  the  phase  change  reaction 
is  endothermic  and  negative  if  the  phase  change  reaction  is  exothermic. 

Finally,  one  has  three  terms  associated  with  the  heat  released  by  the 
three  flames.  The  energy  actually  released  by  each  of  the  individual  flames 
per  unit  mass  of  reactants  is  given  the  designation  Qpp  for  the  primary  flame, 
QAp  for  the  oxidizer  monopropellant  flame  and  Qpp  for  the  final  flame.  The 
actual  quantity  of  energy  transferred  back  to  the  propellant  surface  can  be 
shown  to  be  equal  to  the  quantity  of  energy  released  in  the  flame  times  the 
exp(-(  ),  where  £  is  a  non-dimensional  flame  standoff  distance.  Separate 
energy  balances  can  be  utilized  to  evaluate  each  of  the  three  flame  heat 
release  terms: 

^AP  "  Cp,ox^s  .ox'^o^  +  ^p,g(TAP'^s,ox)  +  (51) 

Qpp  ”  QCp,ox^s  .ox'^o^  +  (l-o)Cp,f(Ts,b"T0) 

+  <1-“>Qf  +  QCp.g<TrTs,ox>  +  a-aKp  gOTp-Ts  b)  (52) 

QPP  -  (l-a)Cpf(Tsb-T0)  +  (l-a)Cpg(Tp-Ts>b) 

+  aCpg(Tp-TAp)  +  (l-a)Qf  (53) 

In  these  expressions,  Tp  and  TAp  are  adiabatic  flame  temperatures  for  the 
overall  propellant  mixture  and  for  pure  AP,  respectively. 


2. 1.4.4  Flame  Standoff  Distances 

The  flame  standoff  distances  are  computed  performing  a  modified  Burke - 
Schumann  flame  analysis  like  in  the  Beckstead,  Derr,  Price  model  described 
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earlier.  However,  there  are  several  important  improvements.  The  modified 
Burke -Schumann  solution  for  the  shape  of  a  Bunsen  burner  type  flame  surface 
is: [21,22] 
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In  the  equation  v 


is  the  stoichiometric  ratio  given  by: 


(54) 


^actual 

u  —  - 

^stoich 


where  the  actual  stoichiometric  ratio  is: 
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(56) 


and  ^sto^ch  is  related  to  the  type  of  fuel  and  oxidizer  considered  and  can  be 
readily  found  within  the  literature.  The  variable  c  is  the  ratio  of  the  inner 
oxidizer  radius,  a,  to  the  outer  fuel  binder  radius,  b.  <)>n  represents  roots 
of  the  first  order  Bessel  function,  J^.  JQ  is  a  Bessel  function  of  the 
zeroth  kind.  ^  is  a  function  of  diffusivity,  gas  stream  velocity,  and  the 
outer  radius  b;  ip  -  D/(vb).  £  is  the  non-dimensional  radial  distance  and  rj  is 

the  non-dimensional  axial  distance. 


In  the  BDP  model  and  early  versions  of  the  PEM  only  the  first  term  of  the 
Bessel  function  series  is  computed  in  equation  (54). [13,24]  This  is  called 
the  short  flame  approximation.  In  later  versions  of  the  Petite  Ensemble 
model,  terms  are  computed  in  the  series  until  the  terms  become  very  small 
yielding  a  more  exact  solution  to  the  diffusional  height,  Xq. [25-28]  Equation 
(54)  is  solved  by  setting  the  non-dimensional  radial  distance,  f,  to  zero  and 
solving  for  the  non-dimension  axial  distance,  tj  . 

The  second  difference  is  that  the  averaged  effective  flame  height  is 
exactly  computed  instead  of  using  the  assumed  constant,  Afh>  of  equation  (44) . 
If  is  assumed  that  the  flame  shape  over  a  oxidizer  surrounded  by  a  fuel  binder 

is  parabolic,  it  can  be  shown  that  a  distributed  heat  release  all  along  the 

parabolic  flame  front  surface  can  be  represented  by  a  planar  flame  sheet  at 
some  fraction  of  the  diffusion  flame  height.  This  parameter,  e,  is  calculated 

from  the  non-dimensional  flame  height,  £  ,  by  the  relation: 

[ln(Z*exp(Z*))-ln(exp(Z*)- 1) ] 

«  -  - T - (57) 

i 

so  the  effective  planar  flame  heights  are  given  by: 

*PD  "  ePDxPD  and  XFD  ”  eFDxFD  (58) 

In  this  expression,  Xpp  and  Xp^  are  the  primary  and  final  diffusional 
flame  effective  planar  flame  standoff  distances  and  Xpp  and  Xpp  are  the  actual 
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diffusion  flame  heights  determined  at  the  center  of  the  two  flames. 

The  final  expressions  for  the  net  flame  standoff  distances  are  the  same  is 
seen  previously  and  are  given  by: 
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The  final  term  needed  for  equation  (50)  is  which  represents  the 
fraction  of  reactants  that  react  in  the  primary  flame.  This  term  is  based 
solely  on  geometry  and  is  computed  in  the  same  manner  as  in  the  BDP  model: 


(45) 


2. 1.4. 5  Discussion 

The  Petite  Ensemble  Model  was  compared  against  AP  polybutadiene  (HTPB) 
propellants  which  contained  up  to  4  distributions  of  AP  particles  ranging  in 
size  from  0.7  to  400  microns .[ 29 , 30]  The  model  successfully  predicted  the 
burning  rate  and  exponent  of  these  formulations  with  over  seventy- five  percent 
of  the  data  being  within  ten  percent  of  the  predicted  values. [25]  The  model 
successfully  predicted  the  effects  of  particle  size  on  propellant  burning  rate 
as  well  as  the  dependence  of  particle  size  distribution  on  burning  rate. [26] 
The  model  has  also  been  used  to  predict  temperature  sensitivity,  and  more 
importantly,  predict  the  effects  of  particle  size  and  propellant  parameters  on 
temperature  sensitivity  of  composite  propellants .[ 27 ] 
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2 . 2  The  Physics  and  Chemistry  of  HMX  Combustion 

2.2.1  Physical  Description  of  HMX 

The  cyclic  nitramine  cyclotetramethylenetetranitramine ,  commonly  called 
HMX,  is  an  important  ingredient  in  propellants  used  in  gun  and  solid  rocket 
propulsion  systems.  HMX  has  also  been  called  1 , 3 , 5 , 7- tetranitro-octahydro 
1,3,5, 7  - tetrazocine ;  1 , 3 , 5 , 7- tetra-nitro  1,3,5, 7- tetrazacyclo- octane ; 
octahydro-1 , 3 , 5 , 7  tetranitro-1 , 3 , 5 , 7- tetrazocine ;  and,  in  the  Soviet 
literature,  octogen.[31] 

HMX  was  first  discovered  in  the  1920' s  simultaneously  in  the  United  States 
and  Germany  where  it  was  found  in  small  amounts  during  the  production  of  the 
explosive  RDX,  a  compound  similar  to  HMX. [32]  Because  HMX  melted  64°C  higher 
than  RDX,  it  was  called  the  "high  melting  explosive"  and,  hence,  the  name  HMX. 
It  wasn't  until  the  early  1940 's  when  large  quantities  of  explosives  were 
required  by  the  war  effort  that  HMX  was  produced  in  any  appreciable 
quantities .[ 32 ]  The  first  important  use  of  KMX  was  the  detonating  explosive 
for  the  first  atomic  weapons  due  to  its  thermal  stability  and  high  energy. 

HMX  is  a  cyclic  molecule  consisting  of  four  methylene  groups  and  four 
nitramine  groups  as  shown  in  Figure  9.  The  chemical  formula  of  HMX  is 
C^HgNgOg .  HMX  in  its  pure  form  is  an  explosive  and  can  be  very  hazardous  due 
its  sensitivity  to  shock.  There  are  several  advantages  in  using  HMX  as  a 
propellant  ingredient.  Due  to  the  high  energy  content  of  HMX,  its  high 
density  and  the  low  molecular  weight  of  gas  produced  during  its  combustion, 
high  values  of  specific  impulse  for  rocket  propellants  and  impetus  for  gun 
propellants  can  be  achieved.  In  contrast  to  propellants  based  on  ammonium 
perchlorate  (AP) ,  nitramine  propellants  do  not  produce  hydrochloric  acid  (HC1) 
unless,  AP  is  incorporated  into  the  propellant  to  serve  as  a  ballistic 
modifier.  Besides  being  corrosive,  HC1  in  the  exhaust  provides  nucleation 
sites  for  moisture  droplets  to  condense  upon,  thereby  producing  a  visible 
contrail  or  secondary  smoke.  The  decomposition,  ignition  and  deflagration  of 
HMX  will  be  described  in  the  following  sections. 

2. 2. 1.1  HMX  Decomposition 

Before  the  combustion  and  deflagration  behavior  of  HMX  can  be  modeled,  it 
is  necessary  to  have  an  understanding  of  the  decomposition  chemistry.  An 
understanding  of  the  detailed  chemical  decomposition  might  well  make  it 
possible  to  suggest  new  types  of  additives  for  combustion  modification. 
Combustion  models  require  decomposition  kinetic  parameters  as  inputs  and  the 
product  distributions  from  the  decomposition  reactions  are  needed  for  detailed 
flame  models. 

The  decomposition  of  solid  HMX  has  been  studied  by  heating  the  HMX  and 
observing  what  evolves  from  the  surface.  Techniques  used  to  measure  the 
decomposition  products  include  mass  spectrometry,  gas  chromatography,  and 
various  laser  diagnostic  techniques, [33-51]  The  major  nitrogen  species  found 
are  N2,  N2O,  NO2,  NO,  and  NH-j  and  the  major  carbon  species  found  are  CH2O 
(formaldehyde),  CO,  CO2  and  HCN.  The  relative  yield  of  the  products  is  a 
function  of  the  decomposition  temperature  and  heating  rate. [52-60]  Figure  10 
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Chemical  Structure  of  the  HMX  Molecule 
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shows  the  change  in  HMX  decomposition  products  produced  as  the  surface 
temperature  is  varied. [33]  In  studies  in  which  the  heating  rates  are  low,  ie. 
less  than  100°C  per  second  and,  the  primary  species  are  CH2O,  NO,  N2O,  CO2  and 
CHO. [33-43]  In  studies  where  the  heating  rates  are  greater  than  1000°C  per 
second  the  primary  species  are  NO,  CO  and  HCN. [44-51]  In  an  actual  burning 
propellant  the  heating  rates  are  extremely  high,  on  the  order  of  10^°C  per 
second. 


Several  reaction  mechanisms  have  been  proposed  on  how  HMX  decomposes. 
Referring  to  Figure  9,  one  scenario  is  that  C-N  bond  fission  initiates  the 
decomposition,  and  that  in  the  gas  phase  the  N-NO2  bond  rupture  is 
important. [61]  A  two  step  flame  is  proposed:  a  near  field,  only  microns  off 
the  surface,  of  exothermic  nitramine  decomposition  to  produce  H2CO,  NO2  and 
N2O  at  a  flame  temperature  of  about  1100 °K;  and  a  far  field  driven  partly  by 
the  oxidation  of  H2CO  by  NO2.  This  scenario  has  been  extended  to  include  a 
fizz  zone  near  the  surface  of  the  H2CO  +  NO2  reaction  and  a  preparation  and 
final  flame  zone  involving  reduction  of  NO  and  oxidation  of  CO. [58]  These 
analyses  do  not  include  HCN  as  an  active  specie. 

A  second  scenario  that  includes  HCN  as  an  important  specie,  is  that  the 
HMX  molecule  initially  eliminates  an  NO2  molecule  which  greatly  destabilizes 
the  ring  (refer  to  Figure  9)  to  produce  rapid  ring  cleavage  forming  3 
molecules  of  H2CNNO2  and  ^CN  +  NO2,  with  the  H2CNNO2  further  decomposing  to 
H2CO  and  N2O  or  ^CN  and  N02-[60]  The  H2CN  rapidly  decomposes  to  H  +  HCN  and 
the  H  attacks  NO2  to  from  NO  and  OH. 

At  low  heating  rates,  l’C/sec,  HMX  deflagrates  at  temperatures  a  few 
degrees  above  the  liquefaction  point  of  280°C. [31,62]  Since  actual  heating 
rates  during  HMX  combustion  are  considerably  higher,  HMX  surface  temperatures 
during  combustion  are  considerably  higher.  Little  information  is  available  on 
the  actual  surface  temperature  during  actual  HMX  combustion,  however,  some 
work  has  been  done  with  embedded  thermocouples  in  which  the  surface 
temperature  was  measured  to  be  around  500°C. [63-65] 

As  HMX  is  slowly  heated  to  between  185 °C  -  1908C,  it  undergoes  an 
endothermic  crystallographic  phase  change  process  from  beta  to  delta  HMX. [62] 
Further  heating  leads  to  slightly  exothermic  decomposition,  beginning  at  about 
271#C.[31]  This  process  is  interrupted  by  a  sharp  endotherm  at  the 
liquefaction  point,  which  is  then  quickly  followed  by  exothermic  deflagration 
at  about  280°C.(33]  The  nature  and  proximity  of  the  liquefaction  and 
deflagration  renders  it  difficult  to  define  heats  of  fusion  and  decomposition 
for  use  in  combustion  models.  If  the  heating  rate  is  slow  enough,  the  sample 
is  observed  to  melt  in  bulk  and  commence  vigorous  boiling  indicative  of  gas 
evolution  in  the  melt  phase.  [31 , 33]  The  liquid  is  brown  in  color,  and  quickly 
disappears  in  a  puff  of  brown  smoke.  The  reaction  from  the  melt  phase  is  an 
exothermic  process. [66] 


The  surface  decomposition  activation  energies  and  frequency  factors  are 
important  parameters  in  combustion  modeling  which  are  needed  in  the  Arrhenius 
expressions  which  determine  the  surface  regression  rate  as  a  function  of 
surface  temperature  (see  equations  28,  29,  48  and  49).  The  measured 
activation  energies  associated  with  decomposition  of  HMX  range  from  10  to  256 
kcal/mole  and  frequency  factors  range  from  1.0x10^  to  unrealistically  high 
values  of  2.7xl02^  1/sec. [67-73]  This  wide  range  of  values  may  be  attributed 
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to  various  causes.  Sources  of  error  include  an  insufficient  temperature  range 
over  which  the  measurement  is  made,  uneven  heating  of  the  sample  of  HMX, 
effects  of  pressure,  sublimation  and  gas  phase  reactions .[ 66 , 70 , 73 ]  In  view 
of  the  relatively  narrow  temperature  range  over  which  the  principal  events 
occur;  ie.  crystallographic  phase  change  at  around  190°C,  solid  phase 
decomposition  at  271°C,  liquefaction  at  284®C  and  finally  deflagrating  between 
285  and  500°C;  it  is  possible  that  a  reported  activation  energy  reflects  an 
average  of  several  processes  which  change  with  temperature  and  physical  state. 
A  statistical  analysis  of  available  decomposition  data  suggests  a  " 
activation  energy  of  53  Kcal/mole  and  a  frequency  factor  of  8.9x10^ 

1/sec. [67]  Where  the  surface  of  deflagrating  HMX  is  a  mixture  of  liquid  and 
gas,  it  may  be  necessary  to  treat  two  processes  of  decomposition  in  combustion 
models . 

The  solid  phase  specific  heat  of  pure  HMX  has  been  determined  using  a 
differential  scanning  calorimeter.  A  value  of  0.25  cal/gr-°K  was  determined 
at  298°K.[74]  The  results  showed  that  the  differences  in  specific  heat  values 
for  single  crystals,  powdered  and  partially  decomposed  samples  were  small. [75] 
In  addition,  the  solid  phase  specific  heat  was  found  to  be  a  linear  function 
of  temperature  around  298  to  450°K.  The  solid  phase  thermal  conductivity  is 
also  a  linear  function  of  temperature  and  has  a  value  of  0.022  cal/cm-sec- °K 
at  298°K. [76,77] 


2. 2. 1.2  HMX  Ignition 

The  ignition  of  HMX  has  been  studied  using  a  xenon  arc -image  furnace  and 
CO2  laser  apparatus  to  apply  known  heat  fluxes  to  the  crystal  surface. [78-80] 
When  a  sample  of  HMX  is  subjected  to  an  energy  flux,  several  processes  occur. 
At  first,  no  reactions  occur.  This  is  a  period  of  inert  heating  of  the  sample 
and  is  shown  as  Region  I  of  Figure  11.(79]  At  some  time,  depending  on  the 
level  of  the  flux  and  the  sample  itself,  the  sample  will  start  to  gasify, 
producing  pyrolysis  products.  The  occurrence  of  this  gasification  is  shown  by 
the  "first  light,  first  gasification"  line  of  Figure  11.  Unlike  the  behavior 
of  AP  propellants  which  almost  immediately  ignite  upon  gasification  at  the 
fluxes  of  most  tests,  HMX  displays  significant  pre-ignition  behavior.  The 
pyrolysis  gases  flow  from  the  surface  but  do  not  immediately  ignite.  This  is 
the  pre-ignition  region,  shown  as  Region  II  of  Figure  11.  If  conditions  are 
favorable,  the  flame  "snaps  back"  to  an  equilibrium  position,  and  the  sample 
is  said  to  have  ignited.  That  is,  the  energy  release  from  the  flame  is 
sufficient  to  further  pyrolyze  the  surface  providing  more  reactants  and  the 
process  can  become  self-sustaining.  The  point  at  which  the  combustion  can  be 
sustained,  upon  interruption  of  the  external  flux,  is  said  to  be  full  or 
go/no-go  Ignition.  The  "go/no-go"  ignition  locus  is  shown  in  Figure  11 
dividing  Regions  II  and  III. [79]  As  is  shown  in  Figure  11,  there  is  a 
significant  delay  between  gasification  of  the  solid  and  ignition  of  HMX. 

The  experimentally  observed  behavior  shown  in  Figure  11  has  also  been 
analytically  modeled  using  two  different  approaches.  The  first  approach  was 
to  use  Bradley's  Unified  Ignition  Theory  which  includes  both  mass  and  thermal 
diffusion  effects. [81]  This  approach,  because  of  the  completeness  of  the 
analytical  description,  is  difficult  to  run  so  a  simpler,  modified  flame  sheet 
analysis  was  also  used. [82]  The  results  of  calculations  using  both  methods 
are  also  indicated  in  Figure  11.  These  analyses  were  then  used  to  predict 
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ignition  behavior  of  HMX  at  other  fluxes  and  pressures .[ 79 ]  It  was  concluded 
that  gas -phase  reactions  and  diffusion  processes,  specifically  the  reaction  of 
pyrolysis  products  to  final  products  and  the  rate  at  which  these  reactions 
occur,  control  the  ignition  event.  Other  HMX  ignition  studies  found  that  the 
time  to  ignite  depended  on  the  atmosphere  and  upon  the  pressure,  further 
evidence  that  ignition  takes  place  in  the  gas  phase. [83,84] 

A  different  study  which  examined  the  ignition  of  catalyzed  RDX  (a  compound 
similar  to  HMX)  concluded  that  the  gas  phase  reactions  were  also  responsible 
for  controlling  the  ignition  behavior.  The  ignition  was  studied  using  a 
reflected  shock  technique  by  shock-heating  N2  at  reflected  shock  temperatures 
of  600-1200°K  and  pressures  of  0.1-1. 5  MPa. [85]  For  temperatures  greater  than 
750°K,  ignition  delays  were  independent  of  pressure  and  almost  independent  of 
temperature.  For  temperatures  less  than  750°K,  temperature  dependence 
increased  markedly  and  the  HMX  ignition  delays  showed  an  inverse  pressure 
dependence.  It  was  concluded  that  for  temperatures  less  than  750°K  the  gas 
phase  reactions  and  diffusional  processes  control  ignition. [85] 
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Figure  11 

Ignition  of  HMX  Pressed  Pellets  Using  Xenon  Arc- image  Furnace 
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2. 2. 1.3  HMX  Deflagration 


Motion  pictures  and  quench  techniques  have  been  applied  to  observe  the 
surface  characteristics  of  deflagrating  HMX.  At  one  atmosphere,  there  is  a 
thick  bubbling  melt  layer.  The  melt  layer  thins  at  higher  pressure,  becoming 
about  4/m  thick  at  1000  psia  and  not  discernible  above  1700  psi.[38]  The 
presence  of  gas  bubbles  in  the  melt  layer  of  deflagrating  HMX  is  observed  in 
electron  microscope  photography  of  extinguished  sample  surfaces. [86] 


At  1700  psi,  the  characteristic  thermal  wave  thickness  in  the  HMX  is  about 
6/im  and  the  heating  rate  becomes  about  10^ °C  per  second,  so  depending  upon  the 
surface  temperature,  a  melt  layer  of  less  than  1/im  would  is  reasonable .[ 62 ] 
Another  interesting  surface  feature  of  deflagrating  HMX  is  the  presence  of 
dark  spots.  These  dark  spots  are  probably  a  polymeric  residue  of 
decomposition. [62]  HMX  is  slightly  fuel-rich  and  similar  spots  have  been 
observed  on  the  surfaces  of  deflagrating  nitrocellulose  propellants  which  are 
slightly  fuel -rich. 


The  gas  phase  at  pressures  below  1000  psia  consists  of  a  transparent  blue 
flame  and  white  smoke. [62 ,80]  The  white  smoke  was  originally  believed  to  be 
condensed  H2CO,  formaldehyde.  Recent  investigations  have  shown  this  to  be 
incorrect.  Using  Coherent  Anti-Stokes  Raman  (CARS)  temperature  measurements 
and  planar  laser  induced  fluorescence,  NO2 ,  NO,  CN,  NH  and  OH  were  seen  in  the 
product  gases  but  there  was  no  evidence  of  H2CO. [48-51] 

Flame  standoff  distances  are  difficult  to  determine  since  the  reaction 
zone  is  not  a  simple  flame  sheet.  At  pressures  above  2000  psia,  the  gas  phase 
consists  of  a  bright  orange  flame. [62,79]  This  orange  color  may  be  due  to 
additional  soot  in  the  combustion  products . [51]  This  change  would  seem  to  be 
significant,  yet  the  burning  rate-pressure  curve  is  not  significantly  altered 
at  ambient  temperature. [31,62] 


The  burning  rates  of  HMX  have  been  measured  using  single  crystals  and 
pressed  pellets  of  HMX  from  atmospheric  pressure  to  50000  psi  and  over  various 
initial  propellant  temperatures . [86-91]  Data  near  one  atmosphere  are 
considered  to  be  uncertain  due  to  the  nature  of  the  bubbling  decomposition  at 
low  pressure. [62]  These  burning  rate,  pressure  and  temperature  dependence 
data  are  needed  by  combustion  modelers  to  check  limiting  conditions  of  the 
model  and  to  deduce  values  of  unknown  constants  by  fitting  the  model  to  the 
data. 


One  interesting  result  that  has  been  observed  in  the  combustion  of  HMX 
are  the  discontinuities  in  the  burning  rate  dependence  on  pressure. [86,88-91] 
Figure  12  shows  the  burning  rate  of  HMX  versus  pressure  from  100  to  50000 
psia. [89]  As  can  be  seen  in  Figure  12,  there  is  an  abrupt  change  in  the 
burning  rate  at  around  700  psi  where  the  pressure  exponent  becomes  very  large 
over  a  narrow  pressure  range  and  elevates  the  burning  rate  by  several  orders 
of  magnitude.  Very  small  HMX  particles  do  not  show  this  behavior  and  have  a 
much  smoother  burning  rate  versus  pressure  curve. [90] 

One  explanation  for  this  burning  rate  behavior  is  attributed  to  the 
formation  of  non-uniform  melt  layer  over  the  surface  of  the  HMX.  Smaller 
particles  are  more  easily  heated  and,  hence,  are  covered  by  a  melt  layer  of 
molten  HMX.  The  melt  layer  over  the  small  particles  prevents  the  permeation 
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of  hot  gases  into  the  HMX  solid  particle  bed.  In  larger  crystals  the  melt 
layer  is  uneven,  thinner  and  more  readily  broken  by  protruding  HMX  crystals. 

As  a  consequence,  depending  on  particle  size  and  beyond  a  given  pressure,  the 
liquid  layer  in  larger  crystals  is  non-uniform  and  allows  hot  gases  to 
diffuses  into  the  HMX  and  increase  the  mass  burning  rate  by  increasing  the 
heat  transfer  to  the  crystals. [92,93] 

A  second  theory,  which  is  an  extension  of  the  first,  states  that  the  added 
heat  transfer  to  the  larger  crystals  causes  the  crystals  to  crack  due  the 
added  thermal  stresses. [94-97]  The  thermal  cracking  increases  the  exposed 
crystal  surface  area.  The  greater  surface  area  increases  the  burning  rate  by 
allowing  even  greater  heat  transfer  into  the  solid. 


PRESSURE.  Kpsi 


Figure  12 

Apparent  Burning  Rates  of  HMX 
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The  above  explanations  indicate  that  the  slope  breaks  shown  in  Figure  12 
result  from  structural  changes  in  the  HMX  crystals,  and  not  changes  in  the 
actual  burning  rate  pressure  dependence  of  pure  HMX.  In  an  effort  to  prove 
that  structural  phenomenon  like  crystal  cracking  is  an  important  aspect  in  HMX 
combustion,  an  attempt  was  made  to  make  tougher,  less  fragile  crystals  by  an 
acetone  pressing  technique . [96 , 98 ]  The  burning  rate  of  these  samples  was 
measured  in  a  closed  combustion  bomb  at  pressures  from  1  atmosphere  to  50,000 
psi.  The  burning  rate  versus  pressure  curve  is  shown  in  Figure  13.  In  Figure 
13,  there  are  no  slope  breaks  and  the  pressure  exponent  is  approximately  one. 


Figure  13 

Pure  HMX  Burning  Rate  Data  with  no  Crystal  Cracking 
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2.2.2  HMX  Composite  Propellants 

Propellants  containing  cyclotetramethylenetetranitramine  (HMX)  are  of 
considerable  interest  for  solid  rocket  propulsion.  The  average  molecular 
weight  of  the  combustion  products  of  composite  nitramine  propellants  is  lower 
than  that  of  more  conventional  AP  based  composite  propellants.  Consequently, 
the  propellant  performance  is  greater  than  AP  propellants  at  equivalent  flame 
temperatures  since  propellant  performance  is  inversely  proportional  to  the 
molecular  weight  of  the  exhaust  products.  Furthermore,  propellants  with  HMX 
crystals  tend  to  have  higher  densities  than  AP  propellants  and,  thus,  their 
density  impulse  is  greater.  Other  advantages  gained  by  using  nitramines  in 
propellant  formulations  include  excellent  thermal  stability,  low  propellant 
flame  temperature  and  non-toxic,  non-corrosive  combustion  products.  However, 
numerous  problems  have  been  encountered  by  the  substitution  of  HMX  for  AP  in 
propellant  formulations .  These  include  high  bum  rate  exponents ,  exponent 
shifts,  low  burning  rate,  and  difficulty  in  tailoring  these  low  burning  rates 
and  high  pressure  exponents.  In  addition,  nitramine  composite  propellants  are 
much  more  difficult  to  ignite  when  compared  to  AP-based  composite 
propellants . [ 99 , 100 ] 

The  pressure  exponent  of  nitramine  composite  propellants  exhibits  complex 
behavior  as  pressure  changes.  The  pressure  exponent  for  HMX  inert  binder 
propellants  ranges  from  0.6  to  0.7  at  pressures  below  3000  psi.  The  value  can 
exceed  1.5  between  3000  and  10000  psi  and  decreases  to  a  value  of 
approximately  1.0  as  pressure  is  increased  beyond  10000  psi. [101-111]  For 
comparison,  AP  composite  propellants  exhibit  exponents  ranging  from  0.3  t">  0.9 
to  3000  psi  and  are  less  sensitive  to  pressure  changes. 

Figure  14  illustrates  the  slope  break  phenomena  typically  observed  for  HMX 
mixed  with  an  inert  binder. [103]  In  Figure  14,  the  burning  rate  is  a  function 
of  the  HMX  particle  size.  The  two  lines  represent  propellants  containing 
coarse  195pm  and  fine  5pm  HMX  particles,  respectively.  At  low  pressures  the 
larger  particle  size  propellant  exhibits  burning  rates  significantly  lower 
than  that  of  the  pure  nitramine . [98 ]  At  about  3000  psi,  a  sudden  increase  in 
the  pressure  exponent  occurs,  and  the  burning  rate  rapidly  approaches  that  of 
the  pure  nitramine.  A  second  discontinuity  then  occurs,  and  thereafter  the 
burning  rate  appears  to  follow  that  expected  for  the  pure  HMX  (see  Figure  13). 
Fine  particle  propellant  exhibits  higher  burning  rate  at  low  pressure.  A 
milder  slope  break  occurs  at  a  pressures  higher  than  the  coarse  particle  slope 
break.  At  very  high  pressures  the  burning  rate  approaches  that  of  the  pure 
nitramine . 

One  explanation  advanced  for  the  slope  break  phenomenon  of  HMX  propellants 
is  that  at  high  pressures  the  HMX  crystals  are  burning  independently  of  the 
propellant  binder . [103 , 104]  At  low  pressures  the  propellant  strand  is  able  to 
burn  as  a  complete  propellant  because  the  inert  binder  melts  and  fills  the 
crystal  interstices,  wetting  and  encapsulating  the  crystals.  A  new  composite 
propellant  is  formed  that  did  not  fully  exist  in  the  original  cured 
propellant.  At  some  higher  pressure  the  burning  rate  of  the  HMX  is  fa.  t 
enough  to  lead  the  development  of  the  binder  melt  such  that  the  the  strand  is 
able  to  burn  as  an  HMX  strand  depending  upon  the  extent  of  crystal  contact  and 
void  interstices  in  the  original  propellant. 
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Figure  14 

Burning  Rate  Characteristics  of  Pure  HMX  and  HMX 
Inert  Binder  Propellant 


Other  explanations  for  the  slope  break  behavior  of  nitramine  propellants 
include  crystal  fracture,  change  from  condensed  phase  to  gas  phase  control 
when  the  thickness  of  the  thermal  wave  approaches  the  crystal  dimension, 
condensed  phase  reaction  at  the  binder  nitramine  interface  and  an  ignition 
delay  for  the  nitramine  crystals  which  becomes  unimportant  for  higher 
pressures. [110]  However,  none  of  these  mechanisms  explain  all  of  the  observed 
features  of  nitramine  propellant  combustion 

Increasing  the  binder  energy  and,  therefore  the  binder  burning  rate,  to  a 
value  close  to  that  of  the  HMX  increases  the  low  pressure  burning  rate.  The 
use  of  energetic  binders  minimizes  the  magnitude  of  the  large  particle  slope 
changes  but  still  results  in  high  pressure  exponents .[ 112 , 113 ]  Bimodal  HMX 
(ie.  a  mixture  of  coarse  and  fine)  propellants  behave  very  much  like 
propellants  containing  fine  HMX  only. [114] 

An  explanation  for  the  small  variation  in  burning  rate  at  pressures  less 
than  3000  psi  is  that  HMX  composite  propellants,  although  heterogeneous  in 
physical  structure,  burn  more  as  a  homogeneous  propellant. [115-117]  During 
burning,  the  crystalline  HMX  melts  together  with  the  polymeric  binder  on  the 
propellant  burning  surface  and  forms  a  chemically  energetic  liquid  mixture. 
Because  of  this  melting,  the  combustion  waves  in  the  gas  phase  are 
homogeneous . 
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The  burning  rates  of  HMX  propellants  are  usually  much  lower  than 
comparable  propellants  containing  AP.  HMX  inert  binder  propellants  have 
burning  rates  of  0.05  to  0.20  in/sec  at  1000  psia  versus  rates  of  0.3  to  1.5 
in/sec  for  similar  inert  binder  propellants  containing  AP. [30,118]  In  order 
to  achieve  0.3  to  1.5  in/sec  variation  in  burning  rate  the  AP  particle  size  is 
adjusted  in  the  propellant.  As  a  general  rule  fine  AP  particles  will  give 
higher  burning  rates.  However,  for  HMX  propellants,  variations  in  particle 
size  have  minimal  effect  on  changing  the  burning  rate  as  can  been  seen  by 
comparing  the  burning  rate  range  above  with  that  of  AP  propellants.  In  order 
to  increase  the  burning  rate,  propellant  formulators  have  turned  to  additives 
such  as  Ferrocene  and  metal  borohydrides . [ 119 , 120]  In  some  cases  the  burning 
rates  were  accelerated  by  a  factor  of  2.2  using  sodium  borohydride  but  the 
pressure  exponent  was  high  and  these  rates  are  still  below  that  of  AP  based 
propellants . [ 120 ] 
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2 . 3  HMX  Combustion  Models 

2.3.1  HMX  Monopropellant  Models 

One  approach  to  the  problem  of  understanding  the  combustion  of  HMX 
propellants  has  been  to  study  the  combustion  of  HMX  itself,  as  a 
monopropellant.  Two  kinetic  flame  models  and  a  Beckstead,  Derr,  Price  type 
model  will  be  discussed  below  which  attempt  to  model  the  combustion  of  HMX  as 
a  monopropellant. 


2. 3. 1.1  Kinetic  Flame  Model 

Detailed  analysis  of  solid  mono-propellant  deflagration  traditionally 
incorporate  a  single  overall  reaction,  with  an  imposed  reaction  order.  This 
type  of  modeling  ensures  that  the  calculated  burning  rate  will  exhibit  the 
correct  pressure  dependence,  but  may  often  be  at  odds  with  observations  of  the 
chemical  kinetics  of  the  process,  and  for  example,  cannot  explain  a  variable 
pressure  exponent.  A  complex  flame  structure,  such  as  that  indicated  for 
nitramines,  cannot  be  expected  to  be  adequately  represented  by  a  single 
overall  reaction  model.  An  HMX  monopropellant  model  based  on  nitramine  flame 
chemistry  was  developed  to  address  these  considerations . [61 , 121-124] 

The  model  postulates  a  high- gradient ,  relatively  thin  near  field  where 
primary  decomposition  dominates  and  an  extended  far  field  dominated  by 
secondary  reactions  where  gradients  are  small  in  comparison.  The  kinetic 
flame  model  demonstrated  that  the  presence  of  the  considerable  weaker 
secondary  reactions  can  enhance  the  primary  decomposition  reaction.  This 
points  to  a  mechanism  of  indirect  coupling  between  the  far  field  processes  and 
the  burning  rate,  through  which  the  far  field  effect  is  amplified.  Figure  15 
shows  the  relative  importance  of  the  near  field  and  far  field  reaction  at 
varying  distances  from  the  propellant  surface .[ 122] 

The  kinetic  flame  model  was  aimed  primarily  at  testing  the  combined 
effects  of  particular  chemical  mechanisms,  kinetics  data  and  available 
thermophysical  properties  rather  than  calculating  the  burning  rate.  In 
addition,  the  model  allows  for  distinct,  simultaneous  gas  phase  reactions. 

The  following  assumptions  were  made.  First,  initial  decomposition  may 
occur  in  both  the  vapor  phase  and  condensed  phase,  following  unimolecular- 
overall  kinetics.  Second,  the  mechanism  in  the  condensed  phase  is  probably 
initiated  by  C-N  bond  fusion,  whereas  in  the  gas  phase  N-NO2  bond  rupture  is 
the  initial  step  leading  consequently  to  different  decomposition  mechanisms 
for  the  vapor  and  condensed  states.  Finally,  a  large  number  of  secondary 
reactions  can  occur,  leading  to  the  final  combustion  product  and  flame 
temperatures . 

In  the  vapor  phase,  two  major  reaction  categories  wore  defined.  The  first 
is  primary  decomposition  reactions  forming  reactive  products  such  as  CH2O, 

NO2 ,  and  N2O.  Second,  these  products  may  undergo  further  reactions  that  are 
somewhat  slower  and  are  typically  of  second  order.  CH2  and  NO2  reactions  lead 
to  final  complete  combustion  products  such  as  H2O,  CO,  CO2,  N2  and  NO.  The 
two  reactions  are  as  follows: 
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NARROW  NEAR-FIELD  REGION: 
HIGH  GRADIENTS 

CQNVECTI VE -D I FFUS I VE- REACT I VE 


Figure  15 

Gas  Phase  Diagram,  Showing  Che  Near  Field  and  Far  Field 
Configuration  for  the  HMX  Deflagration  Wave  (the  near  field 
region  is  magnified  for  clarity) 
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C^HgNgOg  - *  jN02  +  ^20  +  4CH2O  +  2N2  (59) 

NO2  +  ^CH20  - *  NO  +  yCO  +  yC02  +  ^H20  (60) 

All  the  reactions  occur  simultaneously  and  are  thermally  and  chemically 
coupled.  By  introducing  suitable  approximations  for  the  specific  heat  and 
thermal  conductivity  as  a  function  of  temperature,  the  above  reactions 
containing  nine  separate  species  can  be  represented  by  a  non-linear  system  of 
differential  equations  which  are  solved  numerically.  The  consideration  of  two 
overall  chemical  reactions  has  led  to  the  physical  picture  of  the  gaseous 
deflagration  region.  Figure  16  shows  the  mass  fraction  of  species 
participating  in  the  primary  and  secondary  reactions,  equation  (59),  versus 
distance  from  the  surface  at  three  different  pressures .[ 122 ] 

Later  revisions  to  this  model  accounted  for  gas  bubbles  in  the  melt  layer 
and  the  addition  of  a  equilibrium  nitramine  evaporation  law  to  facilitate 
independent  calculation  of  the  mass  burning  rate. [124]  It  was  assumed  that 
the  bubbles  are  small  relative  to  the  overall  melt  layer  thickness;  the 
thermal  relaxation  times  within  the  bubbles  and  between  the  bubbles  and  the 
liquid  can  be  considered  instantaneous;  constant  temperature  exists  throughout 
the  bubbles;  and  the  mean  intrinsic  thermophysical  properties  of  specific 
heat,  thermal  conductivity  and  density  are  constant. 

In  the  model's  original  form  the  experimental  burning  rate  was  an  input 
into  the  model. [121]  To  facilitate  independent  calculation  of  the  burning 
rate,  an  additional  physical  constraint  was  requires.  This  physical  condition 
is  provided  by  the  equilibrium  evaporation  law  at  the  melt/gas  interface.  The 
process  of  evaporation,  at  any  instant,  involves  molecules  of  the  substance 
leaving  the  liquid  in  an  outward  flux,  aout,  and  an  influx  of  molecules 
effectively  returning  to  the  surface,  min.  When  the  net  effect  is  mout  -  min, 
the  vaporization  is  termed  at  equilibrium.  This  was  assumed  to  occur  at  the 
surface  of  deflagrating  nitraraines. 

A  primary  weakness  of  this  type  of  model  is  a  lack  of  complete  and 
accurate  kinetic  reaction  rate  data.  Without  such  data,  experimental  burning 
rates  are  needed  to  curve  fit  the  insufficient  kinetic  data.  Because  of  these 
problems,  the  model  is  largely  empirical  in  nature.  As  more  reaction  rate 
data  are  available,  this  type  of  model  will  become  more  important. 

A  second  flaw  in  the  model  is  the  reaction  scheme  employed  for  the  near 
field,  equation  (59).  The  reaction  includes  CH2O,  formaldehyde,  as  a 
important  specie.  CH2O  has  been  shown  to  not  be  a  dominate  specie  in  the 
decomposition  of  HMX  at  realistic  heating  rates,  and  hence,  burning  rates. 
Formaldehyde  is  present  at  very  low  heating  rates,  but  these  heating  rates  are 
not  representative  of  actual  HMX  combustion. 


53 


NWC  TP  6992 


DISTANCE  FROM  SURFACE,  ^ 


Figure  16 

Near  Field  Distribution  for  HMX  Deflagration 
at  Three  Pressures 
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2. 3. 1.2  Enhanced  Kinetic  Flame  Model 

The  enhanced  kinetic  flame  model  simulating  HMX  combustion  is  an  extension 
of  kinetic  flame  model  described  above. [125-126]  It  builds  upon  the  previous 
work  in  the  following  respects.  First,  it  updates  the  combustion  chemistry  to 
include  more  chemical  reactions.  Second,  it  allows  for  an  alternate  scheme 
based  upon  the  formation  of  ^CN  and  NO2  as  well  as  the  more  conventional 
scheme  where  CH2O  (formaldehyde)  and  ^0  are  formed.  Finally,  several 
simplifications  were  made  to  reduce  the  computer  cost  and  numerical 
convergence  problems . 

Although  the  precise  nature  of  the  chemistry  is  not  known,  the  combustion 
model  assumes  two  possible  reaction  schemes:  one  based  upon  the  formation  of 
CH2O,  and  one  based  upon  the  formation  of  ^CN  or  HCN.  In  the  current  model, 
the  two  schemes  are  treated  separately  for  liquid  phase  and  vapor  phase 
decomposition. 

The  first  scheme,  which  may  involve  C-N  bond  rupture  in  the  liquid  phase 
and  N-N  bond  rupture  in  the  vapor  phase,  is  presented  in  Table  1.  The 
reactions  are  shown  as  sequences  of  global  steps,  with  0  being  the  fraction  of 
HMX  that  vaporizes.  The  reactions  following  liquid  phase  decomposition  follow 
in  a  straightforward  manner,  with  the  assumption  that  the  water-gas  reaction 
is  the  mechanism  by  which  final  products  are  achieved  at  the  end  of  the 
combustion  zone.  The  final  product  distribution  was  established  by 
thermochemical  calculations.  The  sequence  following  vapor-phase  decomposition 
is  based  on  the  premise  that  the  reaction  with  NO2  is  relatively  rapid  and 
that  the  reaction  with  NO  (high  activation  energy)  occurs  at  higher 
temperatures  further  out  in  the  combustion  zone.  The  formation  of  CO2  and  H2 
present  at  thermochemical  equilibrium  is  derived  from  the  water-gas  reaction. 

The  alternative  scheme,  which  may  involve  a  molecular  rearrangement, 
depolymerization  or  N-N  bond  rupture,  is  presented  in  Table  2.  This  scheme  is 
more  difficult  to  establish,  and  is  shown  to  be  more  complicated.  The  first 
step  in  liquid  phase  decomposition  is  a  global  representation  of  the  initial 
formation  of  H2CN  and  NO2,  A  rapid  decomposition  to  HCN  and  reaction  of  the 
H2.  The  subsequent  steps  follow  the  same  reasoning  as  above.  The  first  step 
in  vapor  phase  decomposition  is  also  a  global  representation  of  a  process  that 
includes  H2CN  formation,  but  the  mechanism  produces  CH2O  as  well. 

The  analysis  assumes  constant  thermal  properties.  The  analysis  also 
neglects  the  bubbles  in  the  liquid  layer  that  were  considered  by  the  previous 
kinetic  flame  model. [124]  These  simplifications  help  to  achieve  a  closed  form 
solution  for  the  liquid  layer. 

Calculated  burning  rates  are  shown,  together  with  experimental  data,  in 
Figure  17. [126]  These  calculations  were  based  upon  Table  1  chemistry.  It  is 
observed  that  the  burning  rate,  pressure  exponent  and  temperature  sensitivity 
are  predicted  quite  well  by  this  model.  Values  of  surface  temperature 
calculated  ranged  from  569°K  at  0.2  MPa  to  698°K  at  20  MPa.  Values  of  0 
ranged  from  0.92  to  0.98,  0  decreasing  with  increasing  pressure  because  of  the 
relatively  high  activation  energy  (high  sensitivity  to  temperature)  of  liquid 
phase  decomposition.  Values  of  0  close  to  1  indicate  that  the  major  portion 
of  HMX  decomposition  occurs  in  the  vapor  phase,  and  that  the  combustion  is 
controlled  by  vapor  phase  decomposition.  The  CH2O-N2O  reaction  is  not  a 
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Table  1 

HMX  Combustion  Chemistry:  CH2O  Scheme 

A.  From  liquid  phase  decomposition: 

1.  (l-0)C4H8NgO8  ♦  4<1-/8)CH20  +  4(1-/9)N20 

2.  4(1-)9)CH20  +  4(1-0)NO2  - ♦  4(1-0)  [CO  +  H20  +  N2] 

3.  (1-0)CO  +  (1-/3)H20  ♦  (1-0)CO2  +  (1-0)H2 

B.  From  vapor  phase  decomposition: 

1.  /3C4HgN808  - »  3y9N02  +  3)9N20  +  4/9CH20  +  20N2 

2.  30CH2O  +  3/3N02  - *  ^9N0  +  ^0CO  +  tj/9H20 

3.  ^CH20  +  §/9N20  - ♦  §/3C0  +  §/9H20  +  ^0N2 

4.  ^CH20  +  §0NO  - ♦  rj0CO  +  §/3H20  +  §0N2 

5.  0CO  +  £H20  - ►  pC02  +  0H2 

Table  2 

HMX  Combustion  Chemistry:  HCN  Scheme 

A.  From  Liquid  phase  decomposition: 

1.  a-j0)C4H8N8O8  »  (1-0)  [4HCN  +  2N0  +  2H20  +  2N02] 

2.  2(1-0)HCN  +  2(1-/3)N02  - ♦  (1-0)  [2C0  +  H20  +  NO  +  Jn2] 

3.  2(1-0)HCN  +  3(1-0)NO  - *  (1-0) [2C0  +  H20  +  |n2] 

4.  (1-0)CO  +  (1-0)H2O  ♦  (1-0)CO2  +  (1-0)H2 

B.  From  vapor  phase  decomposition: 

1.  0C4H8N8Og  - ►  30CH2O  +  30N2O  +  0HCN  +  ^0NO 

+  ^0H2O  +  ^N02 

2.  0Ctt2O  +  \0i:O2  ♦  0CO  +  0ti2O  +  J/3N2 

3.  0HCN  +  |/9N20  - ♦  0CO  +  £/3H20  +  20N2 

4.  |0ch2o  +  |^n2o  - ►  |0CO  +  |/SH20  +  |^n2 

5.  ^CH20  +  ^0NO  - *  ^0CO  +  ^0H2O  +  ^0N2 

6 .  fiCO  +  0^2®  - *  PCO2  +  ^2 
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factor  because  too  little  is  formed  in  the  liquid  layer  and  the  reaction  rates 
are  too  slow  at  the  calculated  surface  temperatures.  The  good  agreement  with 
data  was  achieved  by  adjusting  the  kinetics  of  vapor  phase  decomposition. 

Calculations  incorporating  Table  2  chemistry  produced  the  same  results. 
Since  the  same  kinetics  were  used  for  liquid  phase  decomposition,  the 
essential  result  that  /3  is  close  to  1  was  obtained  and  the  fact  that  the 
decomposition  is  now  endothermic  rather  than  exothermic  made  a  negligible 
difference.  However,  the  heat  release  of  the  Table  2  vapor  phase 
decomposition  is  much  less  than  that  of  Table  1,  so  it  was  necessary  to 
increase  Table  2's  vapor  phase  rate  constants  in  order  to  reproduce  the 
burning  rates.  Reactions  between  liquid  phase  decomposition  products  at  the 
surface  continue  to  be  negligible.  Therefore,  meaningful  differences  in  the 
chemistries  will  not  show  up  until  farther  out  in  the  gas  phase. 

The  enhanced  kinetic  flame  model  was  significant  since  it  made  use  of  more 
accepted  chemical  reaction  schemes,  ie.  the  presence  of  HCN  in  the 
decomposition  products.  The  model  was  also  able  to  accurately  predict  HMX 
monopropellant  burning  rates.  However,  Like  the  kinetic  flame  model  described 
earlier,  this  model  relies  heavily  on  reaction  rate  data  which  is  difficult  to 
obtain  and  often  nonexistent. 


PRESSURE,  Mpa 


Figure  17 

Comparison  of  Calculated  and  Experimental 
HMX  Burning  Rates 
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2. 3. 1.3  Beckstead,  Derr,  Price  Monopropellant  Models 

Because  of  the  widespread  use  of  the  Beckstead,  Derr,  Price  model  with  AP 
composite  propellants  it  was  logical  to  attempt  to  apply  the  model  methodology 
to  HMX  monopropellant  combustion. [82 ,127-134]  For  the  monopropellant  case, 
three  basic  equations  are  sed.  The  surface  regression  rate  is  described  by  a 
one-step  Arrhenius  expression: 


r p  -  Aexp ( - E/RTS )  (61) 

where  r  is  the  burning  rate,  p  is  the  propellant  density,  A  is  the  prefactor, 

E  is  the  activation  energy  and  Tg  is  the  surface  temperature.  The  surface 
temperature  is  derived  from  an  energy  balance  at  the  surface  and  is  expressed 
by: 


Ts  “  To  *  <W  +  (Qf/Cp)exp(-Cpmx*/Ag)  (62) 

In  this  expression  TQ  is  the^initial  propellant  temperature,  Ag  is  the  thermal 
conductivity  of  the  gases,  x  is  the  kinetic  flame  standoff  distance,  is 
the  energy  of  decomposition  and  Qf  is  the  gas  phase  energy  release  given  by: 

Qf  -  Cp(Tf  -  T0)  -  Ql  (63) 

where  Tf  is  the  adiabatic  flame  temperature. 

The  HMX  burning  rates  determined  from  these  equations  were  limited  to  a 
narrow  range  of  pressures.  The  simplicity  of  the  equations  did  not  allow  for 
the  complex  behavior  of  HMX  decomposition  and  deflagration.  This  simplified 
Beckstead,  Derr,  Price  approach  did  provide  a  basic  starting  point  to  develop 
more  sophisticated  models  based  on  the  above  three  equations.  A  modified 
Beckstead,  Derr,  Price  model  was  developed  which  attempted  to  accurately 
define  the  self  deflagration  rates  of  HMX  over  the  range  of  initial 
temperatures  from  198 “K  to  423 *K  and  over  the  pressure  range  from  atmospheric 
to  8000  psia. [132] 

The  heat  balance  utilized  in  the  original  BDP  model  included  the  surface 
heat  term,  Q^.  This  was  initially  envisioned  as  a  latent  heat  at  the  surface 
as  in  an  ablating  material.  It  seen  became  apparent  that  condensed  phase 
reactions  were  of  some  importance  to  the  process  and  the  term  became  the 
net  surface  heat  release,  and  could  be  either  exothermic  or  endothermic . 
However,  this  term  was  considered  to  have  a  fixed  value  for  a  given  material, 
in  the  manner  of  a  state  property.  Because  of  this,  two  competing  condensed 
phase  reactions  were  considered,  one  exothermic,  the  other  exothermic.  The 
first  reaction  is: 


C4H808Ng  - -  4CH20  +  4N20  (64) 

This  reaction  is  exothermic  in  nature  with  a  heat  of  reaction  of  170  cal/gra. 
The  second  reaction  is: 


c4h8°8N8  - *  4CHN  +  4no2  +  2h2  (65) 

This  reaction  is  endothermic  and  has  a  heat  of  reaction  of  -482  cal/gm. 
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The  net  result  of  these  reactions  was  lumped  into  a  net  surface  heat 
release  term,  but  it  was  no  longer  a  pseudo-state  material  property.  It  was 
determined  by  the  racio  of  the  two  reaction  rates,  and  becomes  a  rate-variable 
term. 

The  mass  flows  deriving  from  the  endothermic  and  exothermic  reactions  in 
the  condensed  phase  were  cast,  as  done  previously,  in  an  Arrhenius  form: 

MoxF  -  A1  exp( -E1/RTg)  (66) 

Moxg  -  A2  exp(-E2/RTg)  (67) 

where  Mox  is  the  mass  flow  rate  and  the  subscripts  F  and  G  refer  to  the 
exothermic  and  endothermic  reaction  products,  respectively.  Furthermore, 
the  products  resulting  from  these  condensed  phase  reactions  were  treated  as 
reacting  independently  in  the  gas-phase,  leading  to  a  two-flame  model 
analogous  to  the  treatment  of  composite  AP/HTPB  propellants  in  the  original 
model.  The  ratio  of  the  mass  flows  was  determined  by  the  kinetics  of  the 
condensed  phase  reactions,  rather  than  by  the  composition  of  the  propellant. 
The  parameter  f)  was  used  to  represent  the  mass  flow  proceeding  from  the 
endothermic  process.  The  net  surface  heat  release  can  now  be  defined  as: 

Ql  -  /^end  +  ^-^exo  <68> 

In  order  to  retain  the  simplicity  of  the  original  BDP  analysis,  the 
condensed  phase  reactions  are  lumped  at  the  surface,  and  the  gas -phase 
reactions  are  concentrated  into  a  "flame  sheet"  at  a  rate-sensitive  standoff 
distance.  This  leads  to  an  implicit  relationship  involving  a  surface  energy 
balance  where  the  solution  is  reached  by  iterating  on  the  surface  temperature. 
The  energy  equation  is  as  follows: 

? 

I  CpdT  +  0AHend  -  (1-/J)  AHex(J  - 
T 

/3(Qf+AHend)exp(-$*nd)  +  (1-0) (Qf-AHexo)exp(-£*xo)  (69) 


In  this  expression  Cp  is  the  solid  phase  specific  heat,  AH  is  the  heat 
of  reaction,  and  £  Is  the  non-dimensional  flame  standoff  distance. 

Later  versions  of  this  approach  included  distributed  condensed  phase 
reactions. [133,134]  The  condensed-phase  reactions  are  modeled  as  first  order 
distributed  reactions  with  dependence  on  the  vertical  distance  above  the 
propellant  surface.  The  condensed  phase  heat  feed  back  to  the  surface  is  more 
complicated  since  heat  can  be  released  or  absorbed  in  a  distributed  region 
above  the  surface.  A  solution  was  obtained  when  the  heat  transfer  to  the 
surface  equaled  the  heat  required  to  raise  the  propellant  to  the  surface 
temperature.  Figure  18  shows  burning  rate  results  over  a  wide  range  of 
pressures  for  the  deflagration  rate  of  HMX.[132] 
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Figure  18 

HMX  Burning  Rate  Versus  Pressure  at  298 #K 
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2.3.2  HMX  Composite  Combustion  Models 
2. 3. 2.1  Multiple  Flame  Models 

The  Beckstead,  Derr,  Price  multiple  flame  model  has  also  used  to  model 
composite  propellants  containing  HMX. [103,135-138]  Model  modifications  take 
account  of  the  different  surface  and  flame  structures  of  nitramine  propellants 
and  the  different  physical  and  chemical  properties  of  nitramines  relative  to 
AP  propellants  and  AP.  The  assumed  flame  structure  consists  of  an  primary 
flame  and  the  HMX  monopropellant  flame,  as  shown  in  Figure  19.  The  primary 
flame  occurs  between  the  decomposing  HMX  and  the  binder.  There  is  no  assumed 
final  diffusion  flame  between  the  HMX  monopropellant  flame  products  and 
unburned  binder  since  HMX  burns  nearly  stoichiometric  as  a  monopropellant  so 
it  has  no  free  oxygen  in  its  combustion  products . [103 , 135-138] 

New  expressions  were  developed  for  the  area  of  the  exposed  HMX  crystal  to 
describe  both  the  planar  and  cratered  surface  states.  For  the  planer  state, 
the  particle  volume  fraction  was  corrected  by  a  floating  parameter  to  attempt 
to  account  for  the  peripheral  binder  interference.  For  the  planar  state  the 
exposed  oxidizer  surface  area  to  total  surface  area  with  out  this  correction 
is : 


Sox/S0  “  f  (70) 

where  f  is  the  volume  fraction  of  oxidizer. [103]  This  expression  is  much 
simpler  than  the  expression  for  AP  propellants,  see  equation  (33).  Based  on 
correlations  with  experimental  data  there  is  more  interference  with  coarse 
particles  than  with  fine  particles.  Additionally,  an  change  in  the  surface 


Figure  19 

HMX  Multiple  Flame  Structure 
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structure  occurs  if  the  burning  rate  is  such  that  the  thermal  wave  across  a 
spherical  particle  is  too  steep  to  heat  the  bulk  particle  to  its  melting 
point.  The  criterion  is  derived  by  a  simple  transient  heating  analysis  of  a 
spherical  particle.  The  shift  in  the  surface  structure  is  the  proposed 
mechanism  for  a  major  shift  in  pressure  exponent  observed  with  nitramine 
propellants.  The  equation  for  the  area  ratio  that  is  corrected  for  these 
effects  is: 


Sox/So  -  f<KDoa>  <71> 

where  DQ  is  the  particle  diameter  and  K  and  a  are  constants. 

Because  of  the  experimentally  observed  differences  in  surface  temperature 
between  HMX  and  inert  binders,  separate  surface  temperatures  were  used  in  the 
energy  energy  equation: [137,138] 

aCp.HMX(Ts,HMX‘To)  +  (1-Q)cp,b(Ts,b‘To)  ”  +  (1‘a)Qf 

+  ^fQpFexP^ *^PF^  +  (^‘^f )aQHMXexP( "^HMX^  (72) 

The  equation  is  similar  to  the  energy  equations  used  previously  by  the 
Beckstead,  Derr,  Price  multiple  flame  model  and  the  Petite  Ensemble  Model, 
equations  (35)  and  (50),  except  that  there  is  no  heat  feedback  term  from  a 
final  diffusion  flame. 

In  other  respects  the  model  is  identical  to  the  previously  mentioned  BDP 
model,  ie.  the  surface  temperature  is  iterated  on  until  the  heat  transfer  from 
the  HMX  monopropellant  flame  and  HMX  primary  flame  equals  the  heat  needed  to 
raise  the  propellant  from  its  initial  temperature  to  the  surface  temperature. 
Results  from  this  type  of  model  are  shown  in  the  previously  discussed  Figure 
14  by  the  solid  lines. [103]  The  exponent  breaks  are  successfully  modeled  and 
the  burning  rate  agreement  is  very  good. 


2. 3. 2. 2  Time  Averaged  Multiple  Flame  Model 

The  time  averaged  multiple  flame  combustion  model  is  another  derivative  of 
the  Beckstead,  Derr,  Price  multiple  flame  model  but  has  several  underlying 
differences  which  separate  this  model  from  previously  discussed  multiple  flame 
models. [139-142]  The  time  averaged  multiple  flame  model  treats  the  surface 
structure  of  a  HMX  composite  propellant  in  a  completely  different  fashion  than 
Beckstead,  Derr,  Price  combustion  model. 

The  authors  recognized  that  the  combustion  characteristics  of  HMX  when 
compared  with  those  of  AP  exhibit  numerous  differences.  First,  HMX 
monopropellant  burns  30  to  40  percent  faster  with  much  higher  adiabatic  flame 
temperatures  than  AP  monopropellant  at  similar  pressures  and  litial 
temperatures.  The  monopropellart  flame  temperature  for  HMX  is  approximately 
3200°K  compared  to  only  1400#K  for  AP.  Second,  the  adiabatic  flame 
temperature  of  AP  propellants  is  between  2500  and  3000#K  depending  on  the  AP 
concentration.  In  comparison,  the  adiabatic  flame  temperature  of  HMX 
propellants  is  between  1800  and  2200°K.  Finally,  AP  composite  propellants 
typically  burn  at  rates  higher  than  the  AP  monopropellant  rate  with  a  very 
strong  dependence  on  particle  size  and  concentration.  HMX  propellants  burn  at 
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rates  3-10  times  slower  than  their  respective  monopropellant  rates  with  little 
dependence  on  particle  size  or  concentration. 

Differences  in  burning  rates  and  flame  temperatures  are  easily  explained 
since  HMX  burns  nearly  stoichiometrically  as  a  monopropellant.  The 
introduction  of  any  other  substance  such  as  a  fuel  binder  will  lower  the  HMX 
composite  propellant  flame  temperature  and  burning  rate.  The  monopropellant 
AP,  on  the  other  hand,  burns  stoichiometrically  with  the  addition  of  inert 
binder  with  a  resulting  flame  temperature  greater  than  that  for  pure  AP. 

The  time  averaged  multiple  flame  model  authors  phenomenological 
interpretations  of  these  observations  are  outlined  in  Figure  20. [140]  With 
fine  HMX  as  depicted  at  the  bottom  of  Figure  20,  diffusion  mixing  distances 
can  be  expected  to  be  very  short  (i.e.,  less  than  the  monopropellant  flame 
height) .  The  relatively  cool  HMX/binder  primary  flame  dominates  the 
combustion  and  suppresses  the  rate  below  that  of  the  monopropellant. 

However,  with  coarse  HMX  seen  at  the  top  of  the  Figure  20,  the  primary 
flame  should  not  be  able  to  dominate  because  of  the  increased  diffusional 
mixing  distances.  The  authors  initially  assumed  that  larger  crystals  burn  at 
a  rate  approximating  that  of  the  monopropellant,  and  the  propellant  rate  would 
also  be  expected  to  be  close  to  that  of  the  monopropellant.  However,  from 
window  bomb  movies,  coarse  crystals  appear  to  ignite  very  slowly  but  burn  out 
very  rapidly  leaving  depressions  on  the  propellant  surface .[ 143 , 144]  The 


FINE  HMX 


Figure  20 

Phenomenological  Picture  of  HMX  Combustion 
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relatively  short  burn  time  combined  with  the  long  ignition  time  results  in  an 
unstable,  erratic  burning  at  the  propellant  surface.  Statistically,  very  few 
crystals  burn  at  the  same  time,  with  the  end  result  being  a  much  reduced 
burning  rate.  Large  crystals  reside  on  the  surface  as  the  binder  burns  by  and 
then  ignite  and  disappear  rapidly. 

In  summary,  the  rate  controlling  mechanism  for  fine  HMX  crystals  is  the 
relatively  cool  primary  flame  between  the  HMX  and  the  binder.  For  large 
crystals  the  controlling  mechanism  is  the  HMX  ignition  delay  time  coupled  with 
the  fuel  binder  pyrolysis  and/or  combustion  that  controls  the  overall  burning 
rate  of  the  propellant. [140] 

The  time  averaged  multiple  flame  model  uses  a  Beckstead,  Derr,  Price 
multiple  flame  approach  but  instead  of  space  averaging  the  burning  rate,  a 
time  averaging  approach  is  used. [139-142]  In  the  original  BDP  model  the 
burning  surface  area  of  the  fuel  and  various  oxidizer  fractions  were  used 
throughout  the  model.  The  basic  need  for  the  surface  areas  comes  from  the 
space -averaged  overall  continuity  equation: 

mTSo  “  moxSox  +  "bsb  <73> 

Where  m^,  mQX  and  m^  are  the  mass  fluxes  of  the  total  propellant,  the  oxidizer 
and  the  binder  respectively.  SQ,  SQX  and  are  their  respective  burning 
surface  areas.  Solving  for  the  average  linear  burning  rate  yields: 

?  -  ^("WWV  +  ”b<w)  <«> 

so  that  the  surface  areas  are  normalized  with  the  total  burn  surface  area. 

In  addition,  the  oxidizer  surface  area  is  assumed  to  be  a  spherical  sector 
such  as  shown  in  Figure  21.(140]  This  is  probably  a  reasonable  assumption  as 
long  as  the  oxidizer  and  binder  rates  are  similar  and  the  protrusion  or 
recession  relative  to  the  intersection  plane  is  small.  However,  when  the 
rates  differ  significantly  or  when  the  ignition  delay  time  is  too  long, 
unrealistic  surface  areas  result.  Furthermore,  for  such  limiting  conditions 
the  calculated  surface  area  becomes  large,  reaching  a  maximum  at  the 
conditions  shown  in  Figure  21.  This  results  in  an  abrupt  discontinuity  in  the 
calculation  of  the  oxidizer  burning  surface  area  at  the  limiting  conditions. 

AP  propellants  with  very  wide  particle  size  distributions  (ie.  a  propellant 
containing  both  coarse  and  fine  AP  crystals),  double-base  propellants  (double¬ 
base  propellants  contain  nitrocellulose  gelled  with  a  nitrate  ester  such  an 
nitroglycerin) ,  and  propellants  containing  HMX  will  all  result  in  conditions 
where  either  the  oxidizer  binder  rates  are  significantly  different  or  where 
widely  differing  oxidizer  distributions  will  lead  to  the  unrealistic  limiting 
conditions.  Also,  scanning  electron  microscope  (SEM)  examination  of  quenched 
samples  shows  crystals  with  relatively  flat  surfaces,  although  often  with  a 
rounded  protrusion  in  the  center. [140]  Therefore,  an  assumption  of  a  planar 
oxidizer  surface  is  more  realistic  than  a  spherical  sector,  considering  the 
limiting  conditions  shown  in  Figure  21  or  SEM  photographs  of  quenched  samples. 

Utilizing  a  planar  area  derivation,  calculations  were  performed  with  a 
conventional  space  averaged  multiple  flame  model  that  proved  to  be  numerically 
unstable.  The  reason  for  the  instability  can  best  be  visualized  by  examining 
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ORIGINAL  BOP  CALCULATION 


PLANAR  SURFACE 


Figure  21 

Comparison  of  Surface  Area  Calculations  for  a 
Spherical  Sector  Versus  Flat  Sector 


a  map  of  a  normalized  ignition  delay  time  versus  the  ratio  of  oxidizer- to- 
binder  burning  rate  as  shown  in  Figure  22. [140]  In  this  figure  it  is  assumed 
that  the  statistical  intersection  line  drawn  parallel  to  the  surface 
intersects  half  the  particles  in  the  upper  half  and  half  the  particles  in  the 
lower  half.  The  shaded  areas  are  the  regimes  where  finite  SQX  values  are 
calculated.  In  the  nonshaded  areas,  SQX  values  of  zero  are  calculated  for  the 
various  configurations  indicated.  The  upper  shaded  area  corresponds  to  one- 
half  of  the  oxidizer  crystals  burning  and  the  other  half  not  yet  ignited.  The 
lower  shaded  area  corresponds  to  the  other  half  of  the  crystal  ignited  with 
the  first  configuration  burned  out. 

From  this  map,  it  is  apparent  that,  for  large  ignition  delay  times  o:  for 
cxidizer- to-binder  rate  ratios  greater  than  2,  both  oxidizer  configurations 
will  not  be  burning  simultaneously.  Either  one  or  the  other  or  neither 
configuration  could  result.  This  is  apparently  vhat  happens  for  HMX  composite 
propellants . 

Figure  23  is  a  plot  of  burning  rate  versus  surface  temperature  for  A?, 

HMX,  HTF3  (a  typical  propellant  binder)  and  double  base  propellant  using 
typical  activation  energies  and  pyrolysis  parameters  for  each  ingredient. [ 140] 
For  the  range  of  burning  rates  of  interest,  the  differences  between  HMX  and 
HTPB  are  extremely  large.  Their  rates  are  quite  different,  and  one  can  expect 
an  oxidizer  to  binder  ratio  of  5-10  being  the  norm. 
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The  net  result  is  that,  statistically,  few  oxidizer  crystals  are  burning 
on  the  surface  at  any  one  time.  This  can  lead  to  a  very  unstable  surface 
configuration,  which  leads  to  totally  unrealistic  results.  It  can  also  be 
shown  that  the  original  BDP  assumption  of  a  spherical  sector  will  lead  to  very 
similar  unacceptable  results.  Use  of  the  Petite  Ensemble  model  (PEM) 
statistics  will  tend  to  smear  out  the  discontinuities  since  PEM  the  model 
examines  distributions  of  oxidizer  particles  not  just  discrete  particle  sizes. 
For  large  ignition  delay  times  and/or  large  discrepancies  in  the  oxidizer- 
binder  burning  rates,  even  the  PEM  statistics  will  not  fully  resolve  the 
numerical  problems  caused  by  having  empty  and/or  unignited  crystals  on  the 
surface.  Therefore,  the  conclusion  ultimately  was  reached  that  the  space¬ 
averaging  BDP  model  if  used  consistently  will  result  in  numerical 
instabilities  and  unrealistic  physical  limiting  conditions. 

The  approach  of  assuming  a  time-averaged  propellant  surface  is  realistic 
on  a  physical  basis  because  the  binder  and  oxidizer  appear  to  burn  in  series, 
rather  than  in  parallel  as  in  the  original  BDP  model.  The  equation  for  the 
average  propellant  burning  rate  is  of  the  form: 

Y  distances  Y  distances 

r - (75) 

Y  times  Y.  (distances/rates  +  delay  times) 

The  distances  involved  are  an  oxidizer  dimension  and  an  associated  binder 
dimension.  Figure  24  shows  what  occurs  on  a  statistical  basis. [140]  Burning 
will  be  somewhat  random.  However,  for  a  sufficiently  large  sample  it  can  be 
assumed  that  burning  occurs  approximately  in  a  straight  line  through  a  packed 
bed  of  oxidizer/binder .  The  statistical  intersection  diameter  for  straight 
line  combustion  is  the  same  as  used  in  the  previously  described  BDP  and 
Hermance  models : [ 5 , 13 ] 

D'  -  72/5  D0  (19) 

The  intercrystal  thickness  of  binder  associated  with  a  crystal  is  calculated 
assuming  that  the  binder  is  distributed  according  to  the  specific  surface  area 
of  the  individual  particle  sizes  and  that  all  of  the  binder  forms  an  annulus 
around  each  particle .[ 23 ] 

In  other  respects  the  time  averaged  approach  is  similar  to  the  Beckstead, 
Derr,  Price  multiple  flame  model  discussed  previously.  The  two  flames  assumed 
in  the  analysis  are  the  kinetic  HMX  mono -propellant  flame  and  the  primary 
diffusion  flame  between  the  decomposed  HMX  crystal  and  the  binder  (see  Figure 
19).  Unlike  the  Beckstead,  Derr,  Price  model,  separate  surface  temperatures 
were  assumed  for  both  the  HMX  and  the  binder  and  two  energy  equations  were 
written,  one  over  the  HMX  crystal  and  one  over  the  binder.  Each  energy 
equation  balances  the  heat  feedback  from  the  flames  above  the  propellant 
against  the  energy  required  to  raise  the  HMX  and  the  binder  to  their 
respective  surface  temperatures.  The  full  details  of  this  model  will  be 
discussed  later  in  this  report. 

The  model  was  very  successful  in  determining  the  burning  rates  for  a 
series  of  HMX/HTPB  composite  propellants  which  contained  particle  sizes  from 
4/im  to  400/im.[i40]  Parametric  calculations  were  performed  varying  particle 
size,  solids  loading,  and  initial  propellant  temperature.  These  calculations 
agreed  well  with  expected  experimental  and  observed  results. 
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Figure  24 

Time  Averaged  Propellant  Burning 
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CHAPTER  III 

HIGH  ENERGY  PETITE  ENSEMBLE  MODEL  (HYPEM) 


3 . 1  Theoretical  Procedure 


3.1.1  Structure  of  the  Model 


3. 1.1.1  Statistical  Analysis 


The  structure  of  the  High  Energy  Petite  Ensemble  Model,  HYPEM,  as  the  name 
implies  is  based  on  the  Petite  Ensemble  Model  (PEM)  described  previously. [29] 
The  PEM  framework  recognizes  the  fact  that  a  propellant  is  made  up  of  a 
distribution  of  oxidizer  particles.  The  burning  surface  of  a  polydisperse 
propellant  with  many  oxidizer  diameters  is  statistically  rearranged  into 
imaginary  families  of  monodisperse  propellants  which  consist  of  only  one 
particle  size.  These  monodisperse  propellants  burn  as  Individual  "pseudo" 
propellants  and  do  not  interact  physio -chemically  with  the  flames  of 
surrounding  oxidizer  particles .  The  total  burning  rate  is  computed  by 
statistically  adding  up  all  the  pseudo  monodisperse  propellant  burning  rates. 
Real  propellants  are  made  up  of  several  distributions  of  oxidizer  crystals 
each  with  their  own  oxidizer  mass  fraction.  Equations  (46)  and  (47)  presented 
above  describing  the  PEM  model  are  for  propellants  with  only  one  distribution 
of  particles.  The  equation  for  the  distribution  function,  now  called 
Fd  ox,  for  a  particular  particle  distribution  is: 


Fd,ox 


ft  „ 

- 

InD„v - InD_„ 

2- 

OX 

1 

ox  ox 

GX  p 
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Inffox 

(76) 


In  this  equation  aQX  is  the  oxidizer  mass  fraction  of  a  particular  particle 
distribution.  The  constants  DQX  which  represents  the  fifty  percent  weight 
diameter  of  a  unique  oxidizer  mode  aud  <70X  which  is  the  mode  width  parameter 
are  inputs  to  the  model  which  are  determined  by  the  properties  of  the 
particular  oxidizer  grind.  The  total  propellant  burning  rate  is  computed  by 
the  following  expression: 


r 


f(rdFd>dinDox 


EM 


ox 


(77) 
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Where  is  the  total  propellant  oxidizer  mass  fraction,  rd  is  the  pseudo 
propellant  burning  rate  and  Fd  is  the  distribution  function  given  by  equation 
(76).  This  expression  for  the  total  propellant  burning  rate  is  an  integral 
over  all  particle  sizes  for  both  AP  and  HMX  particle  size  distributions. 


3. 1.1. 2  Oxidizer  Mass  Fraction 

A  very  important  parameter  in  characterizing  the  pseudo  propellant  is  the 
amount  of  binder  that  forms  an  annulus  around  individual  oxidizer  particles. 

It  is  assumed  that  the  volume  of  the  fuel  binder  associated  with  each  oxidizer 
crystal  can  be  related  to  the  effective  diameter  of  the  assumed  spherical 
oxidizer  particle  raised  to  a  power.  This  expression  is  as  follows: 

Vb  -  CD5™  (78) 

Where  Vb  is  the  volume  of  binder,  D  is  the  oxidizer  particle  diameter,  XN  is 
assumed  to  be  between  2  and  3  and  C  is  a  constant  that  must  be  solved  for  by 
considering  the  propellants  oxidizer  distribution  along  with  other  overall 
propellant  properties.  In  order  to  develop  an  meaningful  expression  for  the 
constant  C  several  relationships  must  be  developed.  Integrating  Fd  over  InD 
yields : 


DJ 


<|>Fdd2nD  -  aj 


(79) 


Where  aT  is  the  overall  propellant  oxidizer  mass  fraction.  With  this 
knowledge  the  differential  amount  of  oxidizer  mass  about  some  diameter,  D,  per 
unit  mass  of  total  propellant,  drad,  can  be  written  as: 


d“d  “  D 


(80) 


The  differential  number  of  oxidizer  particles  of  some  diameter,  D,  per  unit 
mass  of  total  propellant,  dNd,  can  be  written  as: 


dNd 


dVox,d 

Vox,d 


(dmd//>ox> 


(81) 


Where  pox  is  the  density  of  the  oxidizer  and  VQX  d  is  the  volume  of  one 
oxidizer  particle  of  diameter,  D.  dVQX  d  is  the ’ different ial  volume  of 
oxidizer  of  diameter,  D,  per  unit  mass  of  total  propellant.  From  equation 
(78)  and  (31)  the  differential  amount  of  binder  associated  with  oxidizer 
particles  of  diameter,  D,  per  unit  mass  of  total  propellant,  dVb ,  can  be 
written  as: 


dVb  -  VbdNd  -  CD^dNd  (82) 

Integrating  this  expression  over  all  diameters  yields  the  total  binder  volume 
per  unit  mass  of  total  propellant.  By  introducing  the  overall  propellant 
density,  p p,  and  the  overall  propellant  oxidizer  volume  fraction,  fj,  the 
following  expression  can  be  written: 
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Ppjdvb  -  (i-rT) 


Utilizing  expressions  (80),  (81),  and  (82);  equation  (83)  becomes: 


(83) 


CD^FjjdlnD 

0ox<iD3> 


(i*fx) 


(84) 


Since  C  is  a  constant  in  the  above  expression,  it  can  be  removed  from  the 
integral  and  solved  for: 
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(85) 


Now  that  the  constant  C  is  known  it  will  be  possible  to  write  expressions  for 
the  various  pseudo  propellant  properties  in  terms  of  C  and  known  propellant 
properties.  The  most  important  of  these  is  the  oxidizer  mass  fraction  for 
each  pseudo  propellant. 


Qd  “  - 

%  +  mox 

Where  mox  and  mb  are  the  mass  of  oxidizer  and  fuel  in  the  pseudo  propellant. 
Defining  these  variables  further  yields: 

mb  -  phMb  -  (87) 


1  + 


% 


m. 


(86) 


m 


ox 


<’ox<!'=3> 


(88) 


Finally,  substituting  equation  (85),  (87),  and  (88)  into  equation  (86)  the 
pseudo  propellant  oxidizer  mass  fraction  is: 
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(89) 


The  only  value  that  remains  to  be  defined  to  characterize  the  pseudo 
propellants  is  the  value  for  XN.  If  XN  is  assumed  to  be  equal  to  two,  it  is 
the  same  as  stating  that  the  volume,  or  mass,  of  fuel  binder  associated  with  a 
given  oxidizer  particle  is  related  to  the  surface  area  of  the  oxidizer 
particle.  In  a  like  fashion,  if  XN  is  assumed  to  be  three,  then  the  volume  of 
fuel  surrounding  a  particle  of  oxidizer  is  proportional  to  the  volume  of  the 
oxidizer.  When  the  value  of  XN  is  equal  to  three,  the  analysis  indicates  that 
each  of  the  oxidizer  particle/fuel  binder  pairs,  regardless  of  the  size  of  the 
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oxidizer  particle,  burn  at  the  same  value  of  oxidizer-to-fuel  ratio;  this 
ratio  being  the  overall  oxidizer-to-fuel  ratio  of  the  propellant  formulation. 
Any  value  of  XN  less  than  three  dictates  that  the  small  oxidizer  particles 
within  the  propellant  distribution  burn  rich  while  the  larger  oxidizer 
particles  burn  fuel  lean,  all  with  respect  to  the  overall  oxidizer-to-fuel 
ratio  of  the  propellant.  A  detailed  analysis  of  this  variable  oxidizer-to- 
fuel  ratio  burning  of  the  individual  oxidizer  particles  showed  that  the 
predictive  capability  of  bimodal  propellants  was  enhanced. [145]  Since  there 
is  no  ideal  value  for  this  parameter  it  is  an  HYPEM  model  input  which  can  be 
set  to  any  value  between  two  and  three. 


3. 1.1. 3  HYPEM  Model  Procedure 

From  the  above  analysis  each  pseudo  propellant  is  now  completely 
characterized.  The  particle  size,  volume  of  fuel  allocated  to  each  particle, 
fuel  density,  particle  density,  particle  type  and  input  parameters  such  as 
pressure  and  initial  propellant  temperature  are  the  only  parameters  needed  to 
compute  the  burning  rate  of  the  individual  pseudo  propellants.  These  burning 
rates  of  the  individual  pseudo  propellants  along  with  their  respective 
distribution  functions  will  be  used  in  equation  (77)  to  compute  the  total 
propellant  burning  ’-ate. 

Because  of  the  many  differences  in  the  manner  in  which  AP  and  HMX  crystals 
burn  and  especially  differences  in  the  way  the  crystals  burn  in  a  binder 
medium,  it  was  decided  to  determine  the  burning  rates  of  the  two  types  of 
pseudo  propellants  by  separate  mathematical  approaches.  The  AP  propellants 
are  computed  by  a  modified  Beckstead,  Derr,  and  Price  (BDP)  multiple  flame 
surface  averaging  technique  while  the  HMX  propellants  are  computed  using  a 
time  averaged  approach  developed  by  Beckstead  and  McCarty .[ 13 , 140]  Each 
method,  although  different  in  many  of  the  details,  is  similar  in  structure. 

For  both  methods,  separate  surface  temperatures  must  be  initially  chosen 
for  the  binder  and  respective  oxidizer  crystal.  Then,  the  mono-propellant 
mass  burning  rates  are  computed  by  the  Arrhenius  expression: 

.  ^’^ox^s.ox^ 

mox  “  rox^ox  ”  Aoxe  (90) 


(-VRTs.b) 

®b  “  rb'b  -  Abe  (91) 

where  A  is  the  Arrhenius  frequency  factor,  E  is  the  activation  energy,  R  is 
the  gas  constant  and  Tg  is  the  surface  temperature.  The  subscripts  ox  and  b 
represent  the  oxidizer  and  binder,  respectively. 

Each  pseudo  propellant  is  postulated  to  have  certain  flames  over  the 
oxidizer  crystal  surrounded  by  the  fuel  or  binder  annulus.  For  the  AP 
propellants  their  are  three  flames;  the  AP  mono-propellant  flame,  the  primary 
flame  where  decomposed  products  of  AP  and  binder  react,  and  the  final 
diffusion  flame  in  which  binder  pyrolysis  products  react  with  AP  mono¬ 
propellant  flame  products  (see  Figure  4).  For  the  HMX  propellants  their  are 
only  two  flames;  the  HMX  mono-propellant  flame  and  the  primary  flame  between 
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the  HMX  and  binder  (see  Figure  19).  There  is  no  final  flame  since  the 
products  of  the  HMX  mono-propellant  flame  are  assumed  inert.  For  both  types 
of  pseudo  propellants,  an  energy  equation  is  written  over  the  propellant.  One 
side  of  the  equation  is  the  heat  transferred  from  the  various  flames  to  the 
surface.  This  is  compared  to  the  other  side  of  the  energy  equation  which  is 
the  amount  of  heat  required  to  raise  the  propellant  from  its  initial 
temperature  to  the  surface  temperature.  The  solution  to  the  energy  balance  is 
obtained  by  iterating  on  the  various  surface  temperatures  until  the  two  heat 
transfers  are  equivalent.  As  mentioned  above  the  details  in  computing  the 
pseudo  propellant  burning  rates  different  in  many  ways.  The  two  methods  will 
be  described  in  detail  and  their  differences  and  similarities  examined. 


3.1.2  HMX  Portion 


3. 1.2.1  Time  Averaged  Burning  Rate 


As  mentioned  previously,  very  few  HMX  crystals  burn  at  the  same  time  on 
the  propellant  surface  and  the  use  of  surface  averaging  over  the  total 
propellant  surface  leads  in  unrealistic  results.  A  more  accurate  approach  is 
to  average  the  burning  surface  in  time.  The  HMX/binder  burn  rate  will  be 
computed  using  a  time  averaged  burning  rate  approach  developed  by  Beckstead 
and  McCarty. [31]  The  burning  rate  can  be  defined  as: 

distances  X  distances 

- - (75) 

X  times  X  (distances/rates  +  delay  times) 


rewritten: 
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(92) 


whe  re : 


D'  -  statistical  intersection  diameter  of  HMX  crystal 
A  -  distance  from  HMX  crystal  edge  to  edge  of  assumed  binder 
annulus 

tj^n  -  ignition  delay  time  of  HMX  crystal 

-  binder  burn  through  time  for  a  given  HMX  crystal 

-  burning  rate  of  HMX  crystal 
ijj  -  average  binder  burning  rate 


Dividing  through  by  D'  yields  the  following: 


(1  +  S ) 


(93) 


where:  S  -  2A/D' 
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Incomplete  reaction  of  the  binder  affects  the  burn  rate,  especially  in  HMX 
composite  propellants,  Windowbomb  movies  show  shedding  of  layers  of  binder 
from  the  burning  surface  of  HMX  composite  propellants .( 143 , 144]  Stoichiometry 
indicates  that  HMX  composite  propellants  burn  very  fuel  rich.  It  is  probable 
that  all  of  the  fuel  is  only  partially  oxidized  or,  more  likely,  that  part  of 
the  fuel  does  not  even  react.  The  windowbomb  movies  seem  to  indicate  the 
latter  is  more  accurate;  that  sheets  of  binder  lift  from  the  surface 
unreacted. 

To  simulate  the  effect  of  unreacted  binder  within  the  time  averaged 
burning  rate  equation,  two  terms  will  be  introduced.  The  first  is  a  fraction 
binder  reacted  (fr^)  term  which  will  be  incorporated  into  the  surface  energy 
balances  for  binder  and  oxidizer  and  must  be  incorporated  into  the  time- 
averaged  equation.  The  second,  the  binder  burnthrough  time  (t^) ,  is 
incorporated  into  the  time -averaged  equation  above.  The  final  equation  for 
the  time-averaged  rate  is: 
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3. 1.2. 2  Fraction  of  Binder  Reacting 

Quantitative  data  are  not  available  to  allow  estimation  of  the  fraction  of 
binder  that  does  react  as  a  function  of  pertinent  variables.  Based  on  the 
windowbomb  movies,  binder  shedding  appears  to  be  much  worse  at  low  pressures 
than  at  high  but  did  not  appear  too  dependent  on  particle  size  or  solids 
loading. [ 143 , 144]  Logically,  it  would  seem  that  the  amount  of  binder  reacting 
should  depend  principally  on  the  energy  feedback  from  the  primary  flame  over 
the  binder.  It  would  appear  to  be  related  to  the  initial  temperature  and 
possibly  the  stoichiometry.  Therefore,  the  following  function  is 
assumed: [140] 


f*b  “ 


(95) 


Where  fra  is  an  arbitrary  constant,  0/F  is  the  oxidizer  to  fuel  ratio,  <tst  is 
the  stoichiometric  oxidizer  to  fuel  ratio,  TQ  is  the  initial  propellant 
temperature  and  £pp  is  the  non-dimensional  height  of  the  primary  flame 
associated  with  HMX.  By  matching  the  data  with  frfl  as  an  arbitrary  constant 
it  was  found  to  require  a  value  of  1.2.(140]  The  fraction  of  binder  that 
reacts  increases  with  higher  initial  propellant  temperatures  and  oxygen  rich 
mixtures  and  decreases  with  higher  primary  flame  heights.  The  assumed  form 
seems  reasonable  and  appears  to  give  the  appropriate  dependencies. 


3. 1.2. 3  Binder  Burnthrough  Time 

The  binder  burnthrough  time  is  meant  to  represent  the  unsteady  nature  of 
binder  pyrolysis  in  a  propellant  environment,  particularly  in  fuel-rich 
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situations.  As  an  oxidizer  burns,  the  binder  adjacent  to  the  oxidizer  is 
heated  by  the  primary  flame,  with  the  resultant  pyrolysis  products 
participating  in  the  flame.  However,  as  the  oxidizer  burns  out,  particularly 
in  cases  where  the  oxidizer  rate  is  much  faster  than  the  binder  rate,  there  is 
no  direct  heat  source  to  continue  heating  the  binder,  except  the  energy  stored 
in  the  thermal  wave.  At  this  point,  the  pyrolysis  will  slow  down  and  possibly 
stop.  For  the  remaining  binder  to  burn  or  slough  off,  energy  must  be  robbed 
from  an  adjacent  particle.  Alternatively,  if  underlying  particle  gets  ignited 
from  another  direction,  it  could  burn  under  the  binder,  resulting  in  the 
binder  sloughing  off.  This  effect  has  been  seen  in  the  windowbomb  movies. 
Strahle  proposed  a  format  for  the  binder  burnthrough  concept. [146]  The 
equation  for  the  burnthrough  time  is: 


cb 


(6 
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where  r  is  the  overall  burning  rate  of  the  propellant,  r^  is  the  binder 
burning  rate,  S  is  the  non-dimensional  ratio  of  twice  the  binder  annulus 
thickness  around  the  oxidizer  crystal  to  the  statistical  intersection  diameter 
of  the  oxidizer  crystal,  a  is  the  thermal  diffusivity  and  t^a  is  a 
proportionality  constant.  Thus,  the  burnthrough  time  increases  with 
increasing  interparticle  distance  and  decreases  with  increasing  average  rate. 
The  exponential  terra  approaches  unity  for  most  cases  that  have  been 
considered. 


3. 1.2. 4  Ignition  Delay  Time 

The  postulated  ignition  delay  time  results  from  the  discontinuous  nature 
of  composite  propellants  (i.e.,  the  crystals  are  not  necessarily  in  good 
thermal  contact  with  the  binder) .  The  ignition  delay  time  is  the  time  it 
takes  an  oxidizer  crystal  to  ignite,  assuming  that  the  ignition  process  starts 
when  a  crystal  first  becomes  exposed  to  the  burning  surface.  As  the  burning 
surface  sweeps  over  the  crystal,  heat  is  fed  to  the  crystal,  which  heats  up 
and  starts  to  decompose.  The  time  required  for  this  process  to  occur  is  the 
ignition  delay  time. 

The  first  attempt  to  incorporate  this  concept  into  a  combustion  model  was 
made  by  Hermance  using  data  generated  by  Shannon  for  weakly  pressed  AP 
strands. [5,147]  The  Beckstead,  Derr,  Price  multiple  flame  model  and  other 
derivative  models  have  all  incorporated  the  same  idea  and  form  of  equation, 
which  is: 


cign  "  Cign  p0.721  (97) 

In  most  of  these  models,  varying  the  ignition  delay  time  parameters  has  little 
effect  on  the  net  burning  rate.  However,  the  ignition  delay  of  the  HMX 
crystals  can  be  a  dominant  mechanism  in  determining  the  burning  rate  of  the 
propellant. [140]  Furthermore,  the  previously  used  equation  for  ignition  delay 
was  developed  for  AP,  and  there  are  no  corresponding  data  for  HMX  (pressed 
strands  of  HMX  burn  with  a  relatively  thick  molten  layer  on  the  surface) . 
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To  incorporate  a  diameter  dependence,  Beckstead  and  McCarty  tried  to  apply 
transient  two-dimensional  heat  transfer  analysis.  A  solution  to  the  two- 
dimensional  differential  heat  conduction  equation  was  obtained  from  Cars  law 
and  Jaeger .[ 140 , 148 ]  However,  it  was  concluded  that  the  Carslaw  and  Jaeger 
solution  could  not  be  applied  directly  because  of  nonuniform  boundary 
conditions,  (i.e.,  heating  occurs  on  only  one  side  of  the  crystal).  Beckstead 
and  McCarty  decided  to  use  an  empirical  data  correlation  for  the  ignition 
delay  time.  From  the  derivation  and  from  the  data,  it  was  apparent  that  an 
ignition  delay  time  should  be  proportional  to  the  diameter  raised  to  a  power 
and  inversely  proportional  to  the  rate. [140] 


tign 


(98) 


A  linear  regression  analysis  of  Shannon's  data  for  AP  gives: 


a. 7 


t^gn  -  4. 321sec/cm® ■ ^ 


(99) 


The  similarity  between  equations  (97)  and  (99)  is  readily  apparent.  The 
difference  is  that  use  of  burning  rate  in  place  of  pressure  as  a  variable  is 
more  accurate  considering  the  transient  heat- transfer  equations.  Use  of 
burning  rate  also  permits  application  of  the  equation  to  HMX  crystals.  For 
HMX  it  was  determined  that  a  diameter  exponent  of  around  1.0  gave  the  best 
results . [140]  The  ignition  delay  time  should  also  be  proportional  to  the 
chemical  activity  of  the  oxidizer.  In  other  words,  the  higher  the  activation 
energy  of  the  initial  decomposition  reaction,  the  longer  the  ignition  delay 
time.  After  the  crystal  has  been  heated,  there  is  probable  a  chemical 
induction  time  before  ignition  is  established  completely.  Therefore,  the 
ignition  delay  time  should  be  proportional  to  the  chemical  activity  in  some 
form  or  another  and  to  the  ignition  temperature.  To  simulate  these  effects, 
C£gn  was  defined  as: [140] 


Cign 


Eox^melt‘^o) 

Eref<TrerTo> 


(100) 


where  EQX  is  the  surface  reaction  activation  energy,  Tme^t  the  melting 
temperature  of  the  oxidizer  crystal,  and  Eref  and  Tre£  arbitrary  reference 
values.  In  actuality,  a  reference  Cign  ref  is  input  instead  of  a  Ere£  value. 
Eref  is  then  computed  by: 


Eox^raelt'^o^ 

Eref  “ 

Cign,  ref  ^ref'^o^ 


(101) 


and  equation  (100) 


is 


used  in  to  compute  C^gn 


which  is  used  in  equation  (98) . 
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3 . 1 . 2 . 5  Remaining  Terms 

All  the  terms  in  the  time  averaged  burning  rate  expression  have  been 
defined  except  the  non-dimensional  primary  flame  standoff  distance,  £pp,  and 
determination  of  the  surface  temperatures  which  are  used  by  the  Arrhenius 
expressions  in  equation  (90)  and  (91)  to  determine  the  burning  rates  of  the 
mono -propellant  oxidizer  and  binder,  respectively.  Since  the  methods  to 
compute  these  values  is  very  similar  to  the  AP  portion  of  the  model  they  will 
be  described  later  together  with  the  AP  explanation. 


3.1.3  AP  Portion 


3 . 1 . 3 . 1  Surface  Averaged  Burning  Rate 


The  computation  of  the  AP  pseudo  propellant  burning  rate  is  based  on  many 
of  the  assumptions  evoked  in  the  Beckstead,  Derr,  and  Price  (BDP)  model 
discussed  earlier.  This  especially  true  in  the  description  of  the  surface 
geometry  of  the  oxidizer  crystal  embedded  in  the  fuel  binder.  Unlike  the  HMX 
pseudo  propellants  where  the  burning  rate  was  computed  by  averaging  in  time, 
the  AP  pseudo  propellant  burning  rate  is  averaged  in  space  over  the  surface  of 
the  propellant.  In  describing  the  propellant  surface,  it  is  assumed  that  the 
oxidizer  crystals  are  randomly  mixed  within  the  fuel  binder.  The  statistical 
average  intersection  diameter,  DqX,  is  the  same  as  defined  previously: 

Dqx  -  ^73  Dox  (19) 


As  in  the  BDP  model  and  the  HMX  section  described  above  the  AP  oxidizer  mass 
burning  rate  is  assumed  to  be  represented  adequately  by  an  Arrhenius 


expression  based  on  the  temperature  of  the  oxidizer  surface,  T 


s,ox" 


m 


ox 


-  A 


ox 


exp ( - E 


ox 


/RTs.ox) 


(102) 


In  this  expression,  AQX  and  EQX  are  the  pre- exponential  frequency  factor  and 
the  activation  energy,  respectively,  of  the  oxidizer  surface  decomposition 
reaction.  The  oxidizer  mass  flux,  m^,  has  a  superscript,  T,  denoting  that 
the  flux  is  based  on  total  exposed  oxidizer  surface  area,  SqX.  Since  the 
oxidizer  crystal  can  be  either  protruding  or  recessed,  the  total  oxidizer 


surface  area  can  be  greater  than  the  planar  oxidizer  surface  area,  S 
calculated  from  the  average  intersection  oxidizer  diameter,  D 


OX" 


OX’ 

The  linear 


burning  rate  of  the  oxidizer  surface,  rQX  can  be  determined  from  this  mass 
flux  and  the  density  of  the  oxidizer,  Pox,  as  such: 


m, 


ox 


"OX 


(103) 


ox 


In  a  similar  fashion,  die  fuel  surface  is  assumed  to  pyrolyze  in  a  fashion 
also  represented  by  an  Arrhenius  expression: 


®b  “  “b  “  “b  “  Ab  exp(-Eb/RTs,b> 


(104) 
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In  this  expression,  A^  and  are  the  pre- exponential  frequency  factor  and  the 
activation  energy,  respectively,  of  the  fuel  pyrolysis  reaction.  Since  the 
fuel  surface  is  assumed  to  be  planar,  then  the  mass  flux  based  on  total  fuel 
surface  area,  m£,  can  be  used  interchangeable  with  the  fuel  mass  flux  based  on 
planar  fuel  area,  m£.  In  this  fuel  pyrolysis  equation,  Ts  is  the  surface 
temperature  of  the  fuel  binder.  This  temperature  can  be  either  equal  to  or 
not  equal  to  the  oxidizer  surface  temperature,  Tg  ,  depending  upon  whether  a 
single  or  dual  surface  temperature  analysis  is  employed.  As  was  the  case  with 
the  oxidizer,  the  linear  burning  rate  of  the  fuel,  r^,  can  be  determined  from 
the  fuel  mass  flux  and  the  fuel  density, 


mb 

Pb 


(105) 


The  equation  of  mass  conservation  can  be  written  for  the  burning  propellant 
based  on  the  mass  fluxes  and  surface  areas  of  the  oxidizer  and  binder: 


•PoP 


•P  cP 
raoxSox 


+  "1bsb 


(106) 


Note  that  all  the  mass  fluxes  are  based  on  planar  areas.  The  total  planar 
mass  flux,  m£,  is  desired  since  it  can  be  directly  related  to  the  propellant 
burning  rate,  r,  and  the  total  propellant  density,  p  : 


r  - 


•P 

nix 


(107) 


where : 


PP 


a  (1-a) 

-  +  - 


ox 


P  b 


-1 


(108) 


•  P 

The  mass  conservation  equation  can  be  written  for  m^: 


sL 

.p  .p  °x  . 

mT  -  mox~ +  mb' 

sS 


(109) 


The  total  planar  mass  flux  is  related  to  the  oxidizer  and  binder  planar  mass 
fluxes  by: 


•P  *P  PP  • 

mT  "  moxp  “b" 

ox  b 


(110) 


and  the  density  ratios  are  related  to  the  oxidizer  mass  fraction,  a,  and 
volume  fraction,  f ,  by: 

Pn 


ox  _  °  P ox 


(111) 


Combining  equation  (110)  and  (111)  the  following  can  be  derived: 
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,p  .p  S’ox 


nix  ~  m 


ox 
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(1-Q) 


(112) 


•  P 

The  planar  mass  flux  of  the  oxidizer,  mox,  can  be  related  to  the  total  surface 
mass  flux  of  oxidizer,  m^,  represented  by  equation  (102)  by  the  following: 


r  cT  n 
sox 
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ox 
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ox 


Sp 

•-  sox  -I 


(113) 


The  goal  of  this  section  is  to  find  an  expression  for  the  area  ratio: 


Referring  to  Figure  25,  the  surface  area  of  the  spherical  segment 
representing  the  burning  AP  surface  is : 

Sjx  -  2*R  |h|  (114) 

From  triangle  I: 

R2  -  a2  +  (R  -  h)2 

2R  |h|  -  a2  +  h2  (115) 

From  triangle  II: 

(D0/2)2  -  a2  +  <D0/2  -  L)2 

a2  -  DqL  -  L2  -  L(Dq  -  L)  (116) 

Substituting  equations  (114)  and  (115)  into  equation  (116): 

Sjx  -  w[h2  +  L<Do  -  L>]  (U7> 

where  h  is  a  measure  of  the  non-planarity  of  an  oxidizer  crystal  and  L  is  the 
distance  from  the  top  of  the  original  crystal  surface  to  the  fuel  surface. 

From  equation  (19)  describing  the  statistical  average  intersection  diameter, 

D'  or  2a,  equation  (116)  can  be  rewritten: 


and: 


■rui  D0 
2 


L 


— [  i  ±  jm  ] 
2  L  J 


(118) 


(119) 


Combining  equation  (117)  with  equation  (118)  gives: 
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(120) 


The  terra  involving  h/DQ  can  be  evaluated  using  the  same  ignition  delay  concept 
employed  by  the  HMX  portion  described  previously  and  developed  by  Beckstead 
and  McCarty .[ 140]  Thus,  referring  to  Figure  25: 

h  -  L  -  yQX  (121) 

where  yQX  is  the  linear  distance  that  the  AP  crystal  regresses  and  is  equal  to 
the  product  of  the  linear  burning  rate  of  AP  under  the  given  conditions  and 
the  time  that  the  crystal  burns.  The  burning  time  of  the  crystal  is  equal  to 
the  total  time  necessary  for  the  surface  to  regress  the  distance  L  minus  the 
ignition  delay  time.  Thus: 


^  ^  rox^tb  ‘  Cign^ 

^  ”  L  rox(  r£  *  tign) 

(122) 

where  the  average  regression  rate  of  the  fuel  binder  has  been  used  to  evaluate 
the  total  burning  time.  Equation  (122)  can  be  divided  by  DQ  and  combined  with 
equation  (119)  to  give: 


L 

D„ 


1  - 
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-  J(i  ±  7173) 
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ox 


(123) 


where  the  ignition  delay  time  is  assumed  to  have  the  form  derived  by  Beckstead 
and  McCarty  and  discussed  previously . (140] 


cign 


(98) 


The  +  or  -  refers  to  the  positive  of  negative  particles,  that  is,  whether  the 
lower  hemisphere  of  the  upper  hemisphere  is  cut  by  the  planar  surface.  It  is 
assumed  that  there  are  an  equal  number  of  positive  and  negative  crystals. 

The  planar  surface  area  of  the  oxidizer  crystal  can  be  written  as: 

sox  “  ?D'2  “  Do>2  -  £Do2  <124> 

Dividing  equation  (120)  by  (124)  yields  the  final  relationship  for  the  total 
to  planar  area  ratio  of  the  exposed  oxidizer  crystal: 
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(125) 


This  can  now  be  used  to  determine  the  planar  mass  flux  of  the  oxidizer 
given  by  equation  (113),  and  then  the  total  mass  flux  off  the  propellant 
surface  can  be  determined  with  equation  (112).  The  solution  to  this  set  of 
equations  is  iterative  in  nature.  Given  a  surface  temperature  of  the 
oxidizer,  a  initial  guess  of  the  binder  surface  temperature  is  used  to  compute 
the  binder  burning  rate  used  in  equation  (105).  In  order  for  mass  continuity 
to  be  maintained  equation  (112)  must  be  satisfied.  This  equation  relates  the 
three  planar  mass  fluxes  off  the  surface;  m|,  ix1qX,  and  m^.  The  binder  mass 
flux  is  recomputed  and  the  binder  surface  temperature  computed  that  would  be 
required  to  yield  this  mass  flux.  The  procedure  is  to  iterate  on  binder 
surface  temperature  until  all  equations  are  satisfied.  If  the  oxidizer  and 
binder  are  assumed  to  have  like  surface  temperatures,  mass  continuity  cannot 
be  maintained. 


3. 1.3. 2  Surface  Geometry 


In  the  following  analysis,  it  is  assumed  that  the  diffusion  flames  above 
the  propellant  surface  are  calculated  in  a  manner  based  upon  the  original 
Burke -Schumann  diffusion  flame  analysis .[ 19 ]  In  such  an  analysis,  the  two 
reactant  constituents  are  assumed  to  issue  from  the  surface  in  two  separate 
streams.  The  oxidizer  decomposition  products  flow  from  the  surface  in  a 
circular  duct  centered  above  the  oxidizer  crystal.  The  fuel  binder 
decomposition  products  flow  in  a  concentric  annulus  around  the  oxidizer 
stream.  The  physical  dimensions  of  these  tubes  are  important  in  calculating 
the  diffusion  distances  and  can  be  related  directly  to  the  propellant  surface 
geometry.  The  inner  tube  radius,  a,  is  equal  to  half  the  intersection 
diameter,  DgX,  whereas  the  outer  fuel  radius,  b,  can  be  obtained  from  D^x  and 
the  oxidizer  volume  fraction,  fox.  More  precisely; 


a 


(126) 


also ; 
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'ox 


ox 
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(128) 


•  P 

There  are  now  three  planar  mass  fluxes  which  have  been  calculated;  m£x, 
m£,  and  m£.  Since  the  Burke -Schumann  analysis  requires  that  the  mass  flux 
from  each  tube  be  equal,  the  analysis  holds  the  outer  tube  radius,  b,  constant 
and  permits  the  inner  tube  radius,  a,  to  vary,  via  a  flow  slipline,  so  that 
the  following  conditions  exists: 


•  P  •  P  •  P 

®b  “  mox  “  mT 


(129) 
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The  new  inner  tube  radius,  a',  can  be  shown  to  be  related  to  a  and 
propellant  formulation  variables  as  such: [140] 


a'  -  a 


ox 


1/2 


(130) 


3.1.4  Energy  Balance 


3. 1.4.1  HMX  Energy  Balance 

The  fundamental  equation  which  both  the  HMX  and  AP  portions  of  the  HYPEM 
model  revolve  around  is  the  surface  energy  balance.  The  energy  transmitted  to 
the  surface  from  the  various  flames  situated  above  the  propellant  constituents 
must  equal  the  amount  of  energy  needed  to  bring  the  ingredients  from  their 
initial  temperature,  T  ,  to  their  surface  temperature.  For  the  HMX  pseudo 
propellants,  two  separate  energy  balances  are  used,  one  over  the  HMX  crystal 
and  one  over  the  binder.  The  energy  equation  over  the  HMX  crystal  is: 

raXCs,X(Ts,X'To)  ”  ‘  mXQDX  (131) 

*^XPF 

+  £p<mX,PF  +  “b ,  PF^XPFe 

.  . 

+  (mx  -  mx,PF)QXe 

-  non-dimensional  flame  heights  of  primary  and  HMX  flames 

-  mass  flux  rate  of  HMX  leaving  the  HMX  crystal 

-  mass  flux  rate  of  HMX  reacting  in  primary  flame 

-  mass  flux  rate  of  binder  reacting  in  primary  flame 

-  decomposition  energy  of  HMX 

-  energy  release  in  HMX  primary  flame 

-  energy  release  in  HMX  flame 

-  solid  phase  specific  of  HMX 

-  surface  temperature  of  HMX  crystal 

-  initial  propellant  temperature 

-  fraction  of  heat  from  primary  flame  to  oxidizer 

side  of  this  equation  represents  the  heat  required  to  raise  the 
HMX  crystal  from  its  initial  temperature,  TQ,  to  its  surface  temperature.  The 
first  term  cn  the  right  hand  side  of  equation  (131)  is  the  energy  either 
absorbed  or  released  at  the  surface.  The  final  two  terms  represent  the  heat 
feedback  to  the  propellant  from  the  primary  flame  and  HMX  mono-propellant 
flame.  The  actual  quantity  of  energy  transferred  back  to  the  crystal's 
surface  is  equal  to  the  energy  released  in  the  flame  times  the  exp(-£  ),  where 
£  is  the  non-dimersional  flame  standoff  distance. 

An  assumption  is  made  that  the  binder  and  oxidizer  decomposition  products 
react  in  the  primary  flame  in  stoichiometric  proportions: 

“X,PF  “  ^ST^rb"lb  (132) 


where: 

Tx 

™X,  PF 
“b.PF 
QDX 
%PF 

Qx 

Cs,X 

ls,X 

To 

The  left  hand 


83 


NWC  TP  6992 


Where  is  the  stoichiometric  oxidizer  to  fuel  ratio  and  frb  is  the  fraction 

of  binder  that  reacts  in  the  primary  flame.  In  addition,  the  following 
relationships  are  known: 


mbfrb  “  ^b.PF  and  “x^f  “  ™X,PF  (133) 

is  the  fraction  of  HMX  that  reacts  in  the  primary  flame.  Solving  equation 
(131)  for  Ts  x  and  utilizing  relationships  (132)  and  (133)  the  follow  equation 
results: 


Ls  ,X 


%X 

:s,X 


+  Vft1  +  4?) 


Qxpf 
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Qx  < 


STJ^s.X 


+(1^f)c73(e 


(134) 


In  the  above  equations,  and  in  the  ones  that  follow,  the  solid  phase 
specific  heat  of  HMX,  AP  and  binder  are  computed  from  measured  solid  phase 
specific  heat  versus  temperature  data. [77]  The  temperature  used  to  determine 
the  specific  heat  from  the  data  is  the  average  of  the  surface  temperature  and 
the  initial  propellant  temperature,  TQ .  The  use  of  variable  solid  phase 
thermal  properties  have  improved  the  results  in  similar  multiple  flame 
combustion  models. [149] 

Next,  an  equation  for  the  binder  is  developed.  The  surface  temperature 
equation  for  the  binder  fuel  is  analogous  to  equation  (131)  for  the  HMX 
crystal  and  is: 


”bCSlb<Ts,b-To>  -  -fb%QDb  „  (135) 

'^XPF 

+  (l-/9p)  (mX,  PF  +  ®b ,  PF^XPFe 

Where  m^  and  Cg  b  are  the  mass  rate  of  binder  leaving  the  surface  and  the 
solid  phase  specific  heat  of  the  binder.  The  left  hand  side  of  this  equation 
represents  the  heat  required  to  raise  the  binder  from  its  initial  temperature , 
Tq,  to  its  surface  temperature.  The  first  term  on  the  right  hand  side  of 
equation  (135)  is  the  energy  either  absorbed  or  released  at  the  surface.  The 
final  term  represent  the  heat  feedback  to  the  propellant  from  the  primary 
flame.  As  before,  the  use  of  equations  (132)  and  (133)  into  (135)  will  yield 
a  expression  for  the  binder  surface  temperature,  Ts  b: 

frbQDb  QpF  ‘^XPF 

Ts,b  “  *To  '  ~cs  b  +  <1-0p)frb<*sr*'1>c^e  (136) 

The  energy  release  terms  are  defined  as  follows: 

Qx?f  “  QtCs,X(Ts,X'To)  +  Cp(TXPF'Ts,X>  +  QDxl  (137) 

+  (1-q) (Csb(Tsb-T0)  +  Cp(TxpF-Tsb)  +  QDb] 

%  “  a^Cs  ,X^Ts  ,X"To^  +  Cp(TX'Ts,X)  +  °*DX  (138) 
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In  these  expressions  a  is  the  mass  fraction  of  HMX  in  the  pseudo  propellant, 

Cp  is  the  gas  phase  specific  heat,  T^pp  is  the  HMX  adiabatic  primary  flame 
temperature,  is  the  HMX  mono -propellant  adiabatic  flame  temperature. 

The  decomposition  energy,  QDX ,  is  computed  by  assuming  that  the 
decomposition  involves  two  steps.  First  a  solid  to  liquid  phase  transition 
and  then  a  liquid  to  gas  phase  transition.  The  solid  to  liquid  phase 
transition  energy  is  assumed  constant.  The  liquid  to  gas  phase  transition 
energy  is  a  function  of  pressure  which  has  been  empirically  computed. [ 132 ] 

The  decomposition  energy  for  AP  is  computed  in  a  similar  manner. [18] 

Equations  (134)  and  (136)  for  Tg  ^  and  Ts  ^  are  non-linear  equations  and 
are  solved  for  simultaneously.  This ’adds  considerable  numerical  complexity  to 
the  model.  The  constants  and  p£  which  represent  the  fraction  of  heat  from 
the  HMX  primary  flame  to  the  oxidizer  and  the  fraction  of  HMX  entering  the 
primary  flame  along  with  the  non-dimensional  flame  heights,  fypp  and  will 
be  discussed  in  a  later  section. 


3. 1.4. 2  AP  Energy  Balance 

The  energy  equation  for  the  AP  pseudo  propellants  consists  of  only  a 
single  expression  since  the  AP  portion  of  the  model  averages  the  propellant 
properties  over  the  entire  surface  of  the  pseudo  propellant.  This  equation 
performs  an  energy  balance  over  both  the  AP  oxidizer  crystal  and  the  binder 
simultaneously: 


°Cs,AP<Ts.AP-To>  +  <1-“)Cs,b(Ts,b'To> 


(139) 


"  '“^DAP  *  (1*a)QDb 


-£*  -£*  -{* 

+  (l-/Jf)QApe  AP  +  /JfQppe  PF  +  (l-£f)QFpe  FF 


where  Cs  AP  ia  the  solid  phase  specific  heat  of  AP;  Tg  Ap  is  the  AP  crystal's 
surface  temperature;  Qj)Ap.  is  the  surface  decomposition  energy  of  AP;  QAP. 

Qpp,  and  QPp  are  the  heat  release  terms  associated  with  the  AP  mono-propellant 
flame,  AP  primary  flame  between  decomposed  AP  and  binder  pyrolysis  products 
and  the  final  flame  between  AP  mono-propellant  flame  products  and  unreacted 
binder;  and,  finally,  £AP,  £pp.  and  £pp  are  the  non-dimensional  flame  heights 
associated  with  the  three  flames.  Like  in  the  previous  energy  equations,  the 
left  hand  side  represents  the  sensible  energy  required  to  raise  the  propellant 
ingredients  to  the  surface  temperature.  The  first  term  on  the  right  side 
represents  surface  decomposition  or  pyrolysis  energy  released  of  absorbed  at 
the  surface.  The  remaining  three  terms  include  the  energy  transferred  back  to 
the  propellant  surface  by  the  AP  mono-nropellant  flame,  primary  flame  and  the 
final  flame.  Solving  for  the  AP  surface  temperature  yields. 
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(l-a)cs ,b  “QDAP  (l'a)QDb 

, AP  “  aCs  Ap  ^s.b'^o^  * 


's ,  AP  *  aCs,AP 


a^f)QAp  -e*  -^pF  ci-/*f)QFF  -C 

+  +  «cs;ap  e 


a^s , AP 


The  energy  release  terms  are  as  follows: 

QpF  "  atCs , AP^Ts , AP‘To)  +  Cp(TPF-Ts,AP>] 

+  (l-a)[CSib(Ts,b-T0)  +  Cp(TpF-TSjb  +  QDb) 

QFF  “  aCp , FF<TFF"TAP> 

-f  (l-a)[CSfb(TSfb-T0)  +  Cp(TFF-Tsb)  +  QDb] 
Qap  -  af^s  ,AP^Ts  , AP*To^  +  CP(TAP-Ts,AP>  +  %Apl 


(140) 


(141) 


(142) 


(143) 


In  these  equations  a  is  the  mass  fraction  of  AP  in  the  pseudo  propellant, 
TpF  is  the  adiabatic  primary  flame  temperature,  TFF  is  the  adiabatic  final 
flame  temperature ,  band  TAP  is  the  AP  mono -propellant  adiabatic  flame 
temperature .  These  flame  temperatures  and  those  used  in  the  HMX  energy 
equations  are  calculated  by  using  a  aerothermo-chemistry  program  with  the 
various  propellant  ingredients.  In  the  model  a  table  look-up  routine  is  used 
to  compute  the  flame  temperature  as  a  function  of  pressure  and  oxidizer  mass 
fraction. 

Unlike  the  HMX  section  where  the  two  temperatures  must  be  solved  for  by 
simultaneously  solving  two  energy  equations  for  Tg  ^  and  Ts  b,  The  two  surface 
temperatures  are  coupled  through  the  continuity  of’ mass  relationships 
discussed  previously. 


3.1.5  Flame  Standoff  Distances 

Both  the  HMX  and  AP  mono-propellant  flames  are  kinetic  premixed  flames 
whose  flame  standoff  distance  is  purely  kinetically  controlled.  Both  the  AP 
and  HMX  primary  flames  are  assumed  to  be  partially  kinetic  controlled  and 
partially  diffusion  controlled.  The  total  flame  height  is  assumed  to  the  sum 
of  the  kinetic  and  diffusion  heights.  The  final  flame  over  the  AP  crystal  is 
assumed  to  be  purely  diffusion  controlled  since  the  reacting  species  are 
preheated  by  both  the  AP  mono -propellant  flame  and  the  primary  flame  beneath, 
thus,  the  reaction  kinetics  are  assumed  to  be  very  fast. 
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3. 1.5.1  Kinetic  Flames 

The  general  form  of  the  equations  describing  the  kinetics  and  combustion 
of  a  premixed  flame  can  be  derived  in  a  simple  manner  considering  global 
kinetics.  For  a  uniform  flow  rate,  the  distance  from  the  propellant  surface 
to  the  flame  can  be  approximated  as  the  linear  gas  velocity  multiplied  by  the 
time  it  takes  the  reactants  to  react: 

x*  -  vr  (144) 

where  both  the  velocity  and  the  time  are  averaged  quantities ,  and  it  is 
assumed  that  the  reaction  occurs  in  one  simple  step.  From  continuity 
the  gas  velocity  can  be  evaluated  as: 


“solid 


(145) 


where  p  is  the  gas  density.  For  an  arbitrary  reaction  order  the  reaction 
rate  is : 


w  -  kPS 


(146) 


where  k  is  the  rate  constant,  and  S  the  order.  Integration  of  the  reaction 
rate  to  give  the  reaction  time  results  in  the  proportionality: 


1 

kP7* 


(147) 


Combining  equations  (145)  and  (147)  with  (144)  and  applying  the  ideal  gas  law 
yields  the  standoff  distance  as: 


x* - —  (148) 

kP6 

where  k  now  represents  a  pseudo  rate  constant,  i.e.,  the  rate  constant 
combined  with  the  proportionality  constant  of  equation  (147). 

The  flame  standoff  distance  is  normally  applied  in  a  non-dimensional  form 
as  in  equation  (134),  (136)  and  (140)  which  define  the  various  surface 
temperatures.  From  the  heat  conduction  equation  (i.e.,  the  energy  equation 
written  as  a  second  order  differential  equation)  the  non-dimensional  standoff 
distance  is: 


*  Cpm  * 

i - : — x 


(43) 


where  cp  is  the  specific  heat  and  A  is  the  thermal  conductivity  of  the  gas.  A 
is  defined  as: 
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A 


(149) 


where  Are£  is  the  gas  phase  thermal  conductivity  at  some  reference 
temperature,  Tref. 


3. 1.5. 2  Diffusion  Flames 

Diffusional  processes  are  important  in  the  model  because  of  the  physical 
separation  of  the  binder  and  oxidizer.  As  reaction  times  become  short  with 
increasing  pressure,  the  diffusional  processes  become  significant  and  even 
dominant.  The  diffusional  mixing  distance  is  the  quantity  that  will  be 
applied  in  the  model,  and  will  be  discussed  here. 

A  successful  detailed  analysis  of  a  diffusion  flame  was  reported  by  Burke 
and  Schumann  in  1928. [19]  They  solved  the  mass  conservation  equation  for  a 
Bunsen  burner  type  of  flame.  Because  of  the  circular  symmetry  involved  in 
their  analysis,  it  appears  to  be  very  appropriate  for  the  present  problem  of 
an  assumed  spherical  oxidizer  crystal  surrounded  by  a  fuel.  A  complete  review 
of  the  analysis  and  assumptions  was  made  by  Williams. [25]  One  of  the 
assumptions  involved  neglecting  one  of  the  derivatives  along  the  axis  of  the 
flame,  which  in  effect,  limits  the  assumption  to  tall  flames.  Initial 
calculations  indicated  that  the  flame  height  for  the  propellant  configuration 
and  typical  mass  burning  rates  was  on  the  order  of  the  oxidizer  dimensions. 
Therefore,  the  Burke -Schumann  analysis  was  modified  to  include  short  flames. 


Utilizing  the  nomenclature  of  Williams,  the  diffusion  equation  can  be 
written  as: [25] 

2 


U  +  I  a 

o  r  Jr 
3xz 


v 

D 


(150) 


where  0  is  a  concentration  term  and  the  second  derivative  with  respect  to  x  is 
the  term  that  is  usually  neglected.  D  is  the  gaseous  diffusivity  and  is 
defined  as : 

(T/Tref )X ' 75 

D  ~  Dq -  (151) 

(P/Pref) 

where,  T  is  the  gas  temperature,  F  is  the  pressure,  and  DQ  is  a  reference 
value  of  diffusivity  at  Tref  and  Pref- 

The  Burke -Schumann  solution  of  the  equation  for  the  flame  surface  location 
neglecting  the  second  derivative  in  equation  (150)  is: 


v-  (,l+i/)c^ 
2(l+i/)c 


00 


n-1 


- V*n*>  exP[^h] 

(Wl2  * 


(152) 
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where  v  is  a  stoichiometric  ratio,  J  represents  a  Bessel  function,  <£ 
represents  roots  of  the  first  order  Bessel  function,  £  is  the  non-dimensional 
radial  distance,  and  rj  is  the  non-dimensional  axial  distance.  The  constant  c 
is  related  to  the  burner  geometry  and  is  equal  to  a/b  where  a  is  the  radial 
distance  from  the  center  of  the  flame  to  the  edge  of  the  oxidizer  crystal  and 
b  is  the  radial  distance  from  the  center  of  the  flame  to  the  outer  edge  of  the 
binder  annulus  around  the  crystal  (See  Figure  25).  The  constants  a  and  b  were 
defined  previously  in  equations  (126)  and  (127).  The  non-dimensional 
distances  £  and  ij  are  defined  as: 

€  -  §  (153) 

V  -  (154) 

where  \j>  is  defined  as  D/(vb) .  The  stoichiometric  ratio,  u,  is  the 
stoichiometric  oxygen  to  fuel  ratio  divided  by  the  actual  oxygen  to  fuel 
ratio. 


Solving  this  equation  at  the  center  of  the  flame,  f  equal  to  zero,  and  only 
using  the  first  term  in  the  summation,  the  flame  height  is: 


* 


x 


In 


2(l+i/)cJ1(c*1) 

.  ^-(l+i/)c2^1[J0(^1)]2  . 


(155) 


This  is  the  form  of  the  equation  used  in  the  original  Beckstead,  Derr,  Price 
multiple  flame  model  discussed  previously. 


In  the  modified  solution  to  equation  (150)  including  the  second  derivative 
with  respect  to  the  axial  direction,  the  equation  for  the  location  of  the 
flame  surface  becomes : 


v-  (l+«/)c2 
2(l+v)c 


Ji  <<=<>„> 


Jo(0nOexp 


'-[[l+(2^n)2]1/2-l]r,' 

2  V"2 


(54) 


which  is  identical  to  equation  (152)  with  the  exception  of  the  term  in  the 
exponent  where  $  -  D/(vb) .  This  is  the  equation  which  is  solved  in  the  HYPEM 
combustion  model  and  is  also  used  in  the  Petite  Ensemble  Model. 


There  are  two  possible  solutions  to  this  equation,  one  when  the  flame 
closes  over  the  oxidizer  crystal  and  one  when  it  closes  over  the  binder. 

These  conditions  correspond  to  the  fuel  rich  and  oxygen  rich  conditions, 
respectively.  The  first  solution  is  obtained  by  setting  £  equal  to  zero,  the 
center  of  the  flame.  The  second  solution  is  obtained  by  setting  (  to  1  which 
represents  the  outer  edge  of  the  binder  annulus.  To  determine  which  condition 
to  solve  for,  the  left  hand  side  of  the  equation  (54)  is  examined: 


./-(l-H/)c2 

2(l+«/)c 


(156) 
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If  this  expression  is  positive  equation  (54)  can  only  be  solved  when  |  is  zero 
(flame  closes  over  oxidizer)  and  if  it  is  negative  £  must  be  set  to  1  (flame 
closes  over  binder)  for  a  solution  to  exist.  Either  way,  the  solution  is 
iterative  in  nature  and  terms  in  the  series  are  used  until  they  become  very 
small . 

If  the  assumption  is  made  that  the  flame  has  a  parabolic  shape,  it  can  be 
shown  that  a  distributed  heat  release  all  along  the  parabolic  flame  front 
surface  can  be  represented  by  a  planar  flame  sheet  at  some  fraction  of  the 
diffusion  flame  height,  t.J^29]  This  parameter  is  calculated  from  the  non- 
dimensional  flame  height  £  ,  by  the  relation: 

e  - l*[ln(£*exp(£*))  -  In(exp(£*)  -  1)]  (57) 

The  diffusion  height,  xD,  calculated  from  equation  (54)  is  multiplied  by  e  to 
determine  the  effective  planar  flame  height,  x^: 

xj  -  (157) 


3. 1.5. 3  Flame  Standoff  Distances 


The  various  flame  heights  will  now  be  presented  in  their  final  form  which 
can  be  used  in  the  surface  temperature  equations  (134,  136,  and  140).  Since 
the  final  flame  sits  over  the  AP  crystal  its  flame  height  is  the  sum  of  two 
terms,  the  AP  kinetic  flame  height,  xAp,  given  by  equation  (148)  and  the  final 
flame  effective  diffusion  height,  xFD,  given  by  equations  (54)  and  (157)  and 
normalized  by  equation  (43): 


cpmox 
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cpmox 

AkApP 


c„m 

P  ox  _* 
+  - x 
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(42) 


Similarly,  the  primary  flame  standoff  distance  is  the  sum  of  two  terms.  In 
the  primary  flames,  the  components  must  mix  together  and  then  react  so  that 
there  is  a  diffusion  distance,  Xpp  for  the  AP  primary  flame  and  x^pp  for  the 
HMX  primary  flame,  followed  by  the  kinetic  distance  required  for  the 
ingredients  to  react,  xpp  for  AP  and  x^pp  for  HMX.  Again  in  the  non- 
dimensional  form: 
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(158) 


where  £pp  and  £jrpp  are  the  non-dimensional  distances  of  the  AP  and  HMX  primary 
flames.  The  final  non-dimensional  distances  needed  to  evaluate  the  surface 
temperature  is  that  of  the  kinetically  controlled  oxidizer  flames.  These 
kinetic  flame  heights  are  from  equations  (148)  and  (43)  : 
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3. 1.5.4  The  Competing  Flames 

The  fraction  of  the  total  oxidizer  that  reacts  with  the  primary  flame,  0^, 
can  be  determined  from  the  foregoing  flame  heights.  This  term  is  required 
since  there  is  competition  for  the  oxidizer  decomposition  products  between  the 
primary  flame  and  the  oxidizer  monopropellant  flame.  Based  solely  on 
geometry,  again  assuming  that  the  primary  diffusion  flame  is  parabolic  in 
shape,  this  term  can  be  calculated  with  the  following  expression: [ 29 ] 

<XAP  '  XPF> 

H  -  - * -  (*5) 

XPD 

This  is  the  equation  for  the  AP  gseudo  propellants ,  by  substituting  x^  for 
xAp .  Xxpp  f°r  XPF  anci  xXPD  f°r  XPD  equation  will  determine  0£  for  the  HMX 

pseudo  propellant. 

The  final  parameter  that  must  be  defined  is  0^,  the  fraction  of  heat  from 
the  HMX  primary  flame  to  the  HMX  crystal.  This  parameter  is  needed  in  the 
energy  equations  to  determine  how  the  heat  transfer  from  the  primary  flame  is 
distributed  between  the  HMX  crystal  and  the  binder.  It  is  to  assumed  that  0 
is  proportional  to  0£,  the  fraction  of  HMX  entering  the  primary  flame,  the 
mass  fluxes  at  the  surface,  the  properties  of  the  gas,  and  inversely 
proportional  to  the  primary  flame  height.  Therefore,  the  following  function 
is  assumed. [ 140] 
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xXPF 


(163) 


In  this  equation  a  high  binder  to  oxidizer  flow  rate,  a  high  thermal 
conductivity,  and  a  large  value  of  increases  the  heat  transfer  to  the 
oxidizer  from  the  primary  flame  while  a  increase  in  the  primary  flame  height, 
xXPF>  and  average  specific  heat,  Cpav,  decreases  £p.  £pa  is  a  constant  that 
has  a  value  of  0.2.  In  addition,  equation  (163)  is  constrained  to  never  be 
greater  than  0.8  and  never  less  than  0.05.  Despite  the  complexity  of  the 
expression,  /0p  is  usually  about  1/3  the  value  of 


3.1.6  Pseudo  Propellant  Interactions 

As  will  be  described  later  in  the  model  validation  section,  the  HYPEM 
model  successfully  predicts  the  burn  rate  of  a  wide  variety  of  HMX/HTPB 
propellants  and  AP/HTPB  propellants.  However,  when  the  model  was  used  with 
propellants  containing  both  HMX  and  AP  crystals  the  agreement  with 
experimental  data  was  poor  when  compared  with  results  of  propellants 
containing  only  one  type  of  oxidizer .[ 150]  Two  modif ications  were  made  to  the 
model  to  help  correct  the  problem. 


3. 1.6.1  Mass  Interaction  Parameter 

First,  a  change  was  made  in  the  way  the  pseudo  propellant  mass  fraction 
was  determined.  AP  is  a  true  oxidizer  and  requires  binder  to  burn 
stoichiometric  so  that  all  oxygen  is  reacted  and  a  high  flame  temperature  can 
be  achieved.  HMX,  on  the  other  hand,  burns  stoichiometric  as  a  mono¬ 
propellant  and  the  introduction  of  binder  reduces  the  flame  temperature.  It 
is  reasonable  to  assume  that  AP  crystals  of  comparable  size  will  rob  binder 
from  corresponding  equal  sized  HMX  crystals.  It  was  decided  to  lower  the  AP 
pseudo  propellant  mass  fraction,  thereby,  enriching  the  AP  pseudo  propellants, 
by  introducing  the  mass  interaction  parameter.  This  parameter  multiplies  the 
AP  mass  fraction,  a^,  by  a  constant  less  than  or  equal  to  one.  In  order  for 
continuity  to  be  maintained  the  HMX  mass  fractions  would  have  to  be  multiplied 
by  a  factor  greater  than  one.  These  two  constants  are  obviously  related  and 
depend  on  the  overall  mass  fractions  of  both  oxidizer  types.  Given  the  AP 
mass  fraction  multiplier,  it  can  be  shown  that  the  factor  for  the  HMX  crystals 
is : 


nAPaX 

nx  -  -  (164) 

nApOT  -  aAP 

where  nAp  is  the  AP  mass  interaction  parameter,  n^  is  the  HMX  mass  interaction 
parameter,  q'^  is  the  total  propellant  oxidizer  mass  fraction,  aAp  is  the  total 
AP  mass  fraction  and  a %  is  the  total  HMX  mass  fraction.  In  order  for  the 
modified  pseudo  propellant  mass  fractions  to  remain  less  than  one,  limits  must 
be  placed  on  the  mass  interaction  parameter  such  that  no  pseudo  propellant 
mass  oxidizer  mass  fraction  exceeds  one.  A  value  of  approximately  0.98 
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significantly  improved  the  burn  rate  predictions  for  propellants  containing 
both  AP  and  HMX  crystals. 


3. 1.6. 2  Flame  Temperature  Modifications 

A  second  source  of  error  in  computing  the  burn  rate  for  propellants 
containing  both  AP  and  HMX  crystals  is  the  use  of  discrete  values  for  the 
various  flame  temperatures  used  in  the  energy  equations.  As  discussed  in 
section  2. 3. 2. 2  describing  the  differences  between  AP  and  HMX  combustion,  the 
flame  temperature  of  the  AP/HTPB  mixture  can  be  over  1000°K  higher  that  that 
for  the  HMX/HTPB  mixture.  In  reality  there  is  a  smearing  the  propellant  flame 
temperatures  since  the  flames  are  not  actually  discrete  flames  but  are 
combined  over  the  propellant  surface.  To  account  for  this  smearing,  the 
actual  flame  temperatures  used  in  the  model  are  mass  weighted  by  the  amount  of 
AP  and  HMX  in  the  propellant.  The  mono -propellant  flame  temperatures  for  AP 
and  HMX  were  not  modified  since  these  flames  sit  directly  over  the  crystal  and 
should  be  less  affected  by  other  flames.  The  primary  flames,  on  the  other 
hand,  are  assumed  to  be  at  the  outer  edges  of  pseudo  propellant  and  are  more 
easily  influenced  by  surrounding  flames.  As  with  the  mass  interaction 
parameter,  a  noticeable  improvement  resulted  in  predicting  the  burning  rate  of 
propellants  containing  both  types  of  oxidizer  crystals. 
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3.2  Summary  of  Equations 


The  basic  equations  of  the  HYPEM  composite  combustion  model  have  all  been 
derived.  A  summary  of  the  primary  equations  in  the  model  follows.  First  is 
the  overall  controlling  equations  that  determine  the  overall  burning  rate  of 
the  propellant  and  the  pseudo  propellant  properties: 
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Next  are  the  Arrhenius  equations  which  determine  the  burn  rate  of  the 
individual  propellant  ingredients. 
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The  time  averaged  burning  rate  expression  for  the  HMX  pseudo  propellant: 
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AP  pseudo  propellant  burning  rate  and  surface  geome  ry  equations: 
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The  three  surface  temperature  relationships: 
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Ts,b  “  -To  -  C_  b  +  <1-^p>frb($ST+1)C7Te 


Qpp  -^XPF 
s,b 


(136) 


(l-o)Cs  ,b  “^DAP  (1*Q;)QDb 

Ts,AP  “  ^o  ctCs  Ap  ^s,b'To^  ‘  qCs 


's,AP  Q<^s , AP 


(140) 


(l-^f)QAp  -€*p  ^PF  -C  d-0f)QpF  ■€* 


”'«^P  6 


PF 


+  +  ~^Tt 


FF 


AP 


Equations  describing  the  kinetic  flame  height,  the  non-dimensionalization  of 
the  flame  heights  and  the  diffusional  flame  height  equation: 


* 

x  - 


m 


kP° 


(1*8) 


*  V  * 
(  — — x 


(43) 
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v-(l+v)c2  Ji(c*n)  r-f[l+(2^n)2]1/2.1'|^ 

-  -  )  - r  Jo(^n ^exP  - o -  (5A) 

2(l+u)c  /__  4>n[ J0Wn)]Z  202  J 

n-1 


The  five  flame  standoff  distances  are: 


cpmox  *  _* 

X  (XAP  +  XFD> 


cpmox  cpmox  _* 

+  Xrr 

Ak.nP5  A 


:pmT  *  *  cpmT  _*  cpmT 

-^(x*D  +  x*F)  -  -V“X?D  + 


AkppP0 


__  in»j'  c  _  iii'T* 

P  T  *  .  P  T  _* 


~  <  xXPD  +  xXPF^  " 


~xXPD  + 


A1cydi?P 


(158) 


«Jp  ~ 


cnmnv 
p  OX  -fc 

X  XAP 


cpmox 

AkApP6 


cpmox  * 
— jc — xX 


cpmox 


(159) 


Finally,  the  fraction  of  oxidizer  that  reacts  in  the  primary  flame  and  the 
fraction  of  heat  from  the  HMX  primary  flame  that  is  feed  back  to  the  HMX 
crystal : 

/  *  *  \ 

(xAP  -  xpF) 


av  _* 
Xvt 


®b  l  XPf 
PpaSPl  —  6 


(16?) 
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3.3  HYP EM  Model  Options 

3.3.1  Model  Inputs 

The  HYPEM  composite  combustion  model  is  written  in  standard  FORTRAN  and 
will  run  on  most  computers  supporting  FORTRAN.  Inputs  to  the  model  can  be 
grouped  in  three  categories.  The  first  type  of  inputs  are  the  many  propellant 
constants  such  as  densities,  reaction  orders,  and  specific  heats.  These 
constants  default  to  known  values  and  seldom  need  to  be  modified.  The  second 
type  are  the  propellant  characterization  variables  such  as  oxidizer  particle 
size,  type  (either  AP  and/or  HMX) ,  and  the  overall  mass  fraction  of  the 
various  particles.  Up  to  four  particle  distributions  may  be  specified.  The 
binder  is  assumed  to  be  HTPB.  The  final  type  of  inputs  are  the  control 
variables  which  specify  the  propellant  conditions  to  evaluated  the  burning 
rate.  These  include  pressure  and  initial  propellant  temperature.  Up  to  ten 
values  of  both  pressure  and  initial  propellant  temperature  may  be  input.  In 
addition,  the  model  can  determine  the  burning  characteristics  of  more  that  one 
propellant  per  execution  of  the  computer  code.  A  sample  computer  output  from 
the  HYPEM  model  is  located  in  the  appendix. 


3.3.2  Pressure  Exponent 

If  more  than  one  pressure  is  input  the  burning  rate  pressure  exponent  is 
calculated  by  using  a  least  squares  fit  of  the  equation: 

r  -  cPn  (18) 

where  r  is  the  burning  rate  of  the  entire  propellant  computed  by  equation 
(77),  c  is  a  calculated  constant,  P  is  the  pressure  and  n  is  the  pressure 
exponent.  An  important  factor  in  determining  if  a  combustion  model  is 
correctly  understanding  propellant  behavior  is  if  the  model  not  only  can 
predict  the  correct  burning  rate ,  but  also  the  burning  rate  dependence  on 
pressure,  hence,  predict  the  pressure  exponent. 


3.3.3  Temperature  Sensitivity 


If  more  than  one  initial  propellant  temperature  is  input  the  temperature 
sensitivity  of  the  propellant  is  calculated.  The  temperature  sensitivity  is 
the  sensitivity  of  the  burning  rate  to  changes  in  initial  propellant 
temperature.  This  is  very  important  in  determining  how  a  solid  rocket  motor 
will  behave  if  it  is  conditioned  to  a  initial  temperature.  This  temperature 
can  range  from  -40*F  in  the  Arctic  to  140*F  in  desert  areas.  The  temperature 
sensitivity  is  defined  as: 


a 


P 


(165) 


where  a  is  the  temperature  sensitivity,  r  is  the  burn  rate,  TQ  is  the  initial 
propellant  temperature  and  the  subscript  P  indicates  that  a p  is  computed  at  a 
constant  pressure.  This  parameter  is  automatically  computed  and  averaged  over 
the  initial  temperatures  that  are  input  to  the  model. 
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CHAPTER  IV 

COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


4 . I  Determination  of  Propellant  Constants 

Many  constants  are  required  to  successfully  use  the  HYPEM  combustion 
model.  These  constants  include  physical  properties  of  HMX,  AP  and  HTPB  binder 
such  as  density,  solid  phase  specific  heat  and,  phase  transition  temperatures. 
Constants  such  as  these  are  easily  found  in  chemistry  and  physics  handbooks. 
The  energy  equation  for  both  the  HMX  and  AP  sections  of  the  model  require 
monopropellant  flame  temperatures,  primary  flame  temperatures  and  final  flame 
temperatures  and  the  respective  molecular  weights  of  the  products.  Tables 
were  created  for  all  flames  and  molecular  weights  as  functions  of  oxidizer  to 
fuel  ratio  and  pressure  from  calculations  performed  with  an 
aero thermochemistry  program. [ 151]  These  tables  are  incorporated  into  the 
HYPEM  combustion  model  and  appropriate  values  are  automatically  determined  as 
functions  of  pressure  and  air  to  fuel  ratio  by  the  computer  program. 

Some  values  of  required  constants  are  not  as  readily  known  but  upper  and 
lower  limits  can  be  placed  on  their  values  from  known  values  of  similar 
compounds.  For  example,  the  gas  phase  thermal  conductivity  of  the  product 
gases  has  not  been  measured  but  handbooks  of  physical  constants  do  have 
conductivity  values  for  the  primary  species  found  in  the  products.  For  other 
constants,  the  range  of  measured  values  is  considerably  large,  and  specific 
values  cannot  be  accurately  determined.  For  -Xc-mple,  the  measured  range  of 
the  HMX  decomposition  frequency  factor  ranges  ..  l.OxlO^  to  2.7x10^ 

1/sec. [67-73]  In  this  case  the  many  different  s.iues  were  averaged  and  a 
value  of  near  l.OxlO1  1/sec  was  used. 

As  in  the  case  with  other  steady  state  burning  rate  models,  the  HYPEM 
combustion  model  requires  the  use  of  several  input  constants  such  as 
activation  energies  and  pre-exponential  frequency  factors  for  the  gas  phase 
reactions  considered,  the  heats  of  fuel  pyrolysis  and  oxidizer  decomposition, 
and  the  gas  phase  specific  heats.  Some  of  these  constants  are  known  only  to  a 
small  degree  of  precision.  Many  of  these  parameters  can  only  be  estimated  and 
experimental  burning  rate  and  pressure  exponent  data  for  known  propellant 
formulations  must  be  used. 

An  optimization  program  was  utilized  to  evaluate  some  of  these  lessor 
known  constants.  In  this  manner,  some  of  the  lesser  known  numerical  input 
constants  can  be  varied,  each  within  prescribed  physical  limits,  until  the 
best  fit  to  the  experimental  burning  rate  and  exponent  data  is  obtained.  By 
computing  the  burning  rate  and  pressure  exponent  for  61  propellants,  the 
optimization  program  was  able  to  minimize  the  errors  between  experimental  and 
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theoretically  determined  values.  (These  propellants  will  be  explained  in  the 
next  section.)  This  optimization  process  required  hours  and  sometimes  days  of 
computer  time  to  satisfactorily  determine  a  value  for  a  constant.  Once  a 
satisfactory  value  was  determined,  however,  it  was  fixed  and  was  not  allowed 
to  vary  between  propellants  or  propellant  types  to  improve  the  performance  of 
the  model.  In  the  following  data  comparison  plots  the  model's  input  constants 
are  fixed  for  all  different  propellant  types  and  are  not  allowed  to  vary. 

Table  3  presents  all  the  physical  constants  used  in  the  HYPEM  model. 

Constants  with  a  "yes"  in  the  OPT  column  are  ones  that  were  determined  by  the 
optimization  process. 
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Table  3 

HYPEM  Combustion  Model  Physical  Constants 


Constant  Name 

Value 

Opt 

HMX  crystal  density 

1.900  g/cmA3 

AP  crystal  density 

1.950  g/cmA3 

HTPB  binder  density 

0.930  g/cmA3 

Mass  proportionality  exponent 

3.0 

yes 

HMX  decomposition  frequency  factor 

1.17E+19  g/cmA2-sec 

yes 

HMX  decomposition  activation  energy 

50000.0  cal/mole 

AP  decomposition  frequency  factor 

173359.0  g/cmA2-sec 

yes 

AP  decomposition  activation  energy 

22000.0  cal/mole 

HTPB  pyrolysis  frequency  factor 

300.0  g/cmA2-sec 

HTPB  pyrolysis  activation  energy 

16900.0  cal/mole 

HMX  solid  phase  specific  heat  (298°K) 

0.246  cal/g-K 

change  in  HMX  specific  heat  with  temp. 

.000652  cal/g-KA2 

AP  solid  phase  specific  heat-298°K 

0.2615  cal/g-K 

change  in  AP  specific  heat  with  temp. 

.0005194  cal/g-KA2 

HTPB  solid  phase  specific  heat-298°K 

0.4549  cal/g-K 

change  in  HTPB  specific  heat  with  temp. 

.000593  cal/g-KA2 

HMX  gas  phase  specific  heat 

0.25333  cal/g-K 

yes 

AP  gas  phase  specific  heat 

0.23218  cal/g-K 

yes 

HMX  gas  phase  thermal  conductivity 

.0036112  cal/g-sec-K 

yes 

AP  gas  phase  thermal  conductivity 

.0003990  cal/g-sec-K 

yes 

HMX  or  AP  gas  phase  thermal  diffusivity 

.16244  cmA2-atm/s-K 

HMX  flame  frequency  factor 

13.8  g/cmA3-sec-K 

yes 

HMX  flame  activation  energy 

22021.0  cal/mole 

yes 

HMX  flame  reaction  order 

1.39 

yes 

AP  flame  frequency  factor 

278.5  g/cmA3-sec-K 

yes 

AP  flame  activation  energy 

7929.2  cal/mole 

yes 

AP  flame  reaction  order 

1.64 

yes 

HMX  primary  flame  frequency  factor 

77.3  g/cmA3-sec-K 

yes 

HMX  primary  flame  activation  energy 

16746.0  cal/mole 

yes 

HMX  primary  flame  reaction  order 

1.15 

yes 

AP  primary  flame  frequency  factor 

179.0  g/cmA3-sec-K 

yes 

AP  primary  flame  activation  energy 

20594.0  cal/mole 

yes 

AP  primary  flame  reaction  order 

2.02 

yes 

HMX  total  heat  of  decomposition  @1000psi 

55.0  cal/g 

AP  total  heat  of  decomposition  @1000psi 

260.0  cal/g 

HTPB  total  heat  of  decomposition 

500.0  cal/g 

HMX  primary  flame  stoich.  oxidizer/fuel 

10.143 

yes 

AP  primary  flame  stoich.  oxidizer/fuel 

5.211 

yes 

AF  final  flame  stoich.  oxidizer/fuel 

12.129 

yes 

HMX  ref.  ignition  delay  constant 

0.78 

HMX  ref.  temperature  for  ignition  delay 

675.0  K 

HMX  ignition  delay  diameter  exponent 

1.10 

HMX  melting  point 

553.0  K 

AP  ref.  ignition  delay  constant 

4.321 

AP  ref.  temperature  for  ignition  delay 

480.0  K 

AP  ignition  delay  diameter  exponent 

1.70 

AP  melting  point 

865.0  K 
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4.2  Burning  Rate .  Pressure  Exponent  and  Temperature  Sensitivity 
Prediction 


The  HYPEM  model  has  been  used  to  predict  the  burning  rate,  pressure 
exponent  and  temperature  sensitivity  of  four  types  of  HMX/AP/HTPB  propellants, 
a  total  of  61  different  propellants.  While  all  of  the  following  propellant 
types  contain  HTPB  as  a  binder,  the  type  and  relative  amount  of.  oxidizer 
changes.  Type  I  propellants  contain  only  HMX,  Type  II  propellants  contain 
only  AP,  Type  III  propellants  contain  mostly  HMX  with  5  to  10  percent  AP,  and 
Type  IV  propellants  contain  mostly  AP  with  10  to  15  percent  HMX.  These 
propellants  contain  HMX  particles  from  4pm  to  400pm  and  AP  particles  from 
0.7pm  to  450pm.  The  solids  loading  (amount  of  oxidizer)  ranges  from  70  to  88 
percent.  For  all  of  these  propellants,  the  oxidizer  particle  distribution 
data  are  known  and  will  be  presented  along  with  the  propellant  formulation 
data. 

The  burning  rate  was  computed  for  each  propellant  at  three  different 
pressures  and  at  an  initial  temperature  of  298 °K.  The  pressures  used  were 
dependent  on  where  the  available  experimental  data  was  taken  for  a  given  type 
of  propellant.  The  pressure  exponent,  n,  was  then  determined  by  performing  a 
least  squares  fit  of  equation  (18): 

r  -  cPn  (18) 

The  pressure  exponent  was  then  compared  with  the  experimentally  determined 
value.  The  burning  rate  was  also  computed  at  initial  propellant  temperatures 
of  248,  298  and  348°K  and  at  a  chamber  pressure  of  1000  psi.  The  temperature 
sensitivity,  a  ,  was  then  determined  by  computed  a  between  248  and  298°K  and 
between  298  and  348°K  with  equation  (165): 


and  averaging  the  two  results.  The  temperature  sensitivity  was  then  compared 
with  experimental  values. 

The  HYPEM  model  results  of  each  class  of  propellant  will  now  be  presented 
along  with  comparisons  to  the  actual  measured  burning  rate,  pressure  exponent 
and  temperature  sensitivity. 


4.2.1  HMX  Composite  Propellants:  Type  I 

The  first  class  of  propellants  to  be  examined,  Type  I,  contain  only  HMX 
crystals  dispersed  in  a  HTPB  binder  medium.  These  propellants  were  formulated 
by  Hercules,  Inc. [118]  For  the  series  of  17  propellants,  the  solids  loading 
(amount  of  HMX)  ranges  from  70  percent  to  85.5  percent  and  the  particle  size 
ranges  from  4  to  400  microns. 

Table  4  presents  the  formulation  data  for  these  HMX/HTPB  propellants. 

Table  5  presents  the  experimental  data  and  the  theoretical  prediction  for  the 
burning  rate  at  600,  1000  and  1600  psi,  the  pressure  exponent  at  an  initial 
temperature  of  298°K  and  the  temperature  sensitivity  at  1000  psia.  Figures 
26,  27  and  28  plot  the  theoretical  HYPEM  values  on  the  vertical  axis  versus 
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the  experimental  values  on  the  horizontal  axis  for  the  burning  rates  at  1000 
psi,  the  pressure  exponents,  and  the  temperature  sensitivities,  respectively, 
shown  in  Table  5.  On  each  plot  there  is  a  number  by  each  data  point  which 
corresponds  to  the  propellant  identification  number  in  Table  4  and  5.  Error 
bands  shown  by  the  dashed  lines  are  included  in  the  Figures  26,  27  and  28. 

The  burning  rate  and  pressure  exponent  plots  have  10  percent  error  bands  while 
the  temperature  sensitivity  plots  have  20  percent  error  bands.  The  solid  line 
represents  perfect  agreement  between  experiment  and  model  prediction. 

As  can  be  seen  in  Figures  26  and  27,  the  HYPEM  model  predicts  the  burning 
rates  and  exponents  for  the  majority  of  these  propellants  within  ±10  percent. 
Figure  28  indicates  that  the  HYPEM  model's  burning  rate  temperature 
sensitivity  prediction  consistently  over  estimates  the  experimental  prediction 
by  more  than  20  percent. 


Table  4 

HMX/HTPB  Propellant  Formulations:  Type  I 


No 

%HTPB 

D 

a 

%HMX 

■ 

200/i 

147  /i 
1.86 

58n 

42.4/x 

2.68 

4/i 

4 . 46fi 
1.99 

1 

30.00 

- 

_ 

70.00 

_ 

2 

30.00 

70.00 

- 

- 

- 

3 

30.00 

- 

- 

35.00 

35.00 

4 

30.00 

35.00 

- 

35.00 

- 

5 

30.00 

- 

23.33 

23.33 

23.33 

6 

22.00 

- 

- 

- 

78.00 

7 

22.00 

- 

- 

78.00 

- 

8 

22.00 

- 

78.00 

- 

- 

9 

22.00 

78.00 

- 

- 

- 

10 

22.00 

- 

39.00 

39.00 

- 

11 

22.00 

- 

26.00 

26.00 

26.00 

12 

14.50 

- 

28.50 

28.50 

28.50 

13 

14.50 

- 

38.50 

29.90 

17.10 

14 

14.50 

- 

42.75 

- 

42.75 

15 

14.50 

- 

64.10 

- 

21.40 

16 

14.50 

- 

42.75 

42.75 

- 

17 

30.00 

- 

- 

7C  .00 

Distribution  data  unavailable,  the  values  shown  are  assumed. 
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Table  5 

Predicted  and  Experimental  Burning  Rates ,  Exponents  and 
Temperature  Sensitivities  for  HMX/HTPB  Propellants:  Type  I 


No 

l 

Burning  Rate  -  in/sec 

600psia  lOOOpsia  1600psia 

ex  rth  rex  rth  rex  rth 

Exponent 

nex  nth 

Temperature 

Sensitivity 

%/°K** 

CTP . ex  CTp , th 

1 

0. 

.048 

0. 

,050 

0. 

.066 

0, 

,064 

0. 

,088 

0. 

.086 

0. 

.618 

0. 

545 

_ 

2 

0. 

.030 

0, 

.040 

0. 

.041 

0. 

.056 

0. 

,067 

0. 

.074 

0. 

.816 

0. 

,630 

0. 

26 

0. 

32 

3 

0. 

.070 

0. 

.072 

0, 

.095 

0. 

.095 

0. 

.125 

0. 

,127 

0. 

,591 

0. 

,578 

- 

- 

4 

0. 

.042 

0, 

.045 

0. 

.057 

0. 

,060 

0. 

,077 

0. 

,080 

0. 

.618 

0. 

.584 

0. 

22 

0. 

34 

5 

0. 

.060 

0, 

.060 

0. 

.083 

0. 

.081 

0. 

.113 

0. 

,109 

0. 

.645 

0. 

603 

• 

- 

6 

0. 

.093 

0. 

.120 

0. 

.128 

0. 

.160 

0. 

,173 

0. 

,214 

0. 

633 

0. 

591 

- 

- 

7 

0 

.074 

0, 

.073 

0. 

.093 

0. 

.095 

0. 

.116 

0. 

,120 

0. 

.458 

0. 

,509 

- 

- 

8 

0. 

.052 

0. 

.053 

0 

.073 

0. 

.075 

0. 

,099 

0. 

,104 

0. 

657 

0. 

690 

- 

- 

9 

0. 

.054 

0. 

.054 

0. 

.073 

0. 

.075 

0. 

.112 

0. 

,100 

0. 

.742 

0. 

.619 

- 

- 

10 

0. 

.070 

0. 

.063 

0 

.092 

0. 

.085 

0. 

,119 

0. 

.112 

0. 

.541 

0. 

589 

- 

- 

11 

0. 

.084 

0. 

.082 

0. 

.120 

0. 

.110 

0. 

,166 

0. 

.146 

0. 

,695 

0. 

,590 

- 

- 

12 

0. 

.103 

0. 

.122 

0 

.143 

0. 

.163 

0. 

.195 

0. 

.213 

0. 

.651 

0. 

.567 

0. 

18 

0. 

25 

13 

0 

.097 

0. 

.112 

0. 

.134 

0. 

.149 

0, 

,180 

0. 

,195 

0. 

,630 

0. 

567 

0. 

13 

0. 

26 

14 

0, 

.104 

0. 

.125 

0 

.144 

0. 

.169 

0, 

.195 

0. 

.225 

0. 

.641 

0. 

.601 

0. 

12 

0. 

24 

15 

0 

.095 

0 

.105 

0 

.129 

0. 

.143 

0. 

.170 

0. 

.191 

0. 

.593 

0. 

.612 

- 

- 

16 

0, 

.092 

0 

.100 

0 

.127 

0. 

.133 

0. 

.177 

0. 

.172 

0. 

,667 

0. 

550 

- 

- 

17 

0 

.072 

0 

.094 

0 

.103 

0 

.125 

0. 

.140 

0. 

.168 

0. 

.677 

0. 

.595 

0. 

.21 

0. 

33 

The  subscript  ex  refers  to  experimental  data 
the  subscript  tn  refers  to  theoretical  model  prediction 
The  temperature  sensitivity  was  computed  by  evaluating  the 
burning  rate  at  initial  temperatures  of  248,  298,  and-348°K 
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HMX/HTPB  PROPELLANTS 

1000  PSI  -  10  PERCENT  ERROR  BAND 


0.00  0.05  0.10  0.15  0.20 

EXPERIMENTAL  BURNING  RATE  -  IN/SEC 


Figure  26 

Theoretical  Versus  Experimental  Burning  Rate 
for  HMX/HTPB  Propellants:  Type  I 
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HMX/HTPB  PROPELLANTS 

10  PERCENT  ERROR  BAND 


0  0.2  0.4  0.6  0.8  1 

EXPERIMENTAL  EXPONENT 


Figure  27 

Theoretical  Versus  Experimental  Exponent 
for  HMX/HTPB  Propellants:  Type  I 
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HMX/HTPB  PROPELLANTS 

1000  PSI  -  20  PERCENT  ERROR  BAND 


EXP.  TEMPERATURE  SENSITIVITY  -  %/K 


Figure  28 

Theoretical  Versus  Experimental  Temperature 
Sensitivities  for  HMX/HTPB  Propellants:  Type  1 
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4.2.2  AP  Composite  Propellants:  Type  II 

The  second  class  of  propellants  to  be  examined,  Type  II,  contain  only  AP 
crystals  dispersed  in  a  HTPB  binder  medium.  These  oropellants  were  also 
formulated  by  Hercules,  Inc. [30]  For  the  series  of  21  propellants,  the  solids 
loading  is  constant  at  87.37  percent  and  the  particle  size  ranges  from  0.7  to 
450  microns . 

Table  6  presents  the  formulation  data  for  these  AP/HTPB  propellants. 

Table  7  presents  the  experimental  data  and  the  theoretical  prediction  for  the 
burning  rate  at  500,  1000  and  2000  psi,  the  pressure  exponent  and  an  initial 
temperature  of  298°K  and  the  temperature  sensitivity  at  1000  psia.  Figures 
29,  30  and  31  plot  the  theoretical  HYPEM  values  versus  the  experimental  values 
for  the  burning  rates  at  1000  psi,  the  pressure  exponents,  and  the  temperature 
sensitivities,  respectively,  shown  in  Table  7.  On  each  plot  there  is  a  number 
by  each  data  point  which  corresponds  to  the  propellant  identification  number 
in  Table  6  and  7.  As  was  previously  seen,  the  burning  rate  and  pressure 
exponent  plots  have  10  percent  error  bands  while  the  temperature  sensitivity 
plots  have  20  percent  error  bands.  The  solid  line  represents  perfect 
agreement  between  experiment  and  model  prediction. 


Table  6 

AP/HTPB  Propellant  Formulations  -  12.63%  HTPB:  Type  II 


D 

No  a 

%AP 

400^ 

448m 

1.22 

200m 

195m 

1.63 

90m 
71. 2m 
1.37 

50m 

44.2m 

1.45 

20m 

22.6m 

1.68 

6m 

5.23m 

1.88 

2m 

1.89m 

1.31 

•  7m 
0.69m 
2.72 

1 

- 

- 

31.58 

- 

13.68 

- 

42.11 

2 

- 

- 

- 

- 

55.79 

- 

- 

31.58 

3 

- 

31.58 

- 

- 

24.21 

- 

- 

31.5 

4 

42.11 

- 

- 

- 

13.68 

• 

- 

31.5 

5 

- 

- 

31.58 

- 

13.68 

31.58 

10.53 

- 

6 

- 

- 

- 

31.58 

24.21 

- 

31.58 

- 

7 

- 

31.58 

- 

- 

24.21 

- 

31.58 

- 

8 

42.11 

- 

- 

- 

13.68 

• 

31.58 

- 

9 

- 

- 

42.11 

- 

13.68 

31.58 

- 

- 

10 

- 

31.58 

- 

- 

24.21 

31.58 

- 

- 

11 

42.11 

- 

- 

- 

13.68 

31.58 

- 

- 

12 

- 

31.58 

- 

31.58 

24.21 

- 

- 

- 

13 

- 

- 

31.58 

- 

55.79 

- 

- 

- 

14 

- 

- 

42.11 

- 

45.26 

- 

- 

- 

15 

- 

31.58 

- 

- 

55.79 

- 

- 

- 

16 

42.11 

- 

- 

- 

45.26 

- 

- 

- 

17 

31.58 

31.58 

- 

10.52 

13.68 

- 

- 

- 

18 

31.58 

- 

- 

42.11 

13.68 

- 

- 

- 

19 

- 

42.11 

- 

31.58 

13.68 

- 

- 

- 

20 

- 

31.58 

- 

42.11 

13.68 

- 

- 

- 

21 

42.11 

- 

- 

31.58 

13.68 

- 

- 

- 
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As  can  be  seen  in  Figures  29  and  30  the  HYPEM  model  predicts  the  burning 
rates  and  exponents  for  the  majority  of  these  propellants  within  ±10  percent. 
Figure  31  indicates  that  the  HYPEM  model's  burning  rate  temperature 
sensitivity  prediction  within  ±20  percent. 


Table  7 

Predicted  and  Experimental  Burning  Rates,  Exponents  and 
Temperature  Sensitivities  for  AP/HTPB  Propellants:  Type  II 


Burning  Rate  -  in/sec 

Exponent 

Temperature 

Sensitivity 

500psia  lOOOpsia  2000psia 

%/°K 

No 

rex  rth  rex  rth  rex  rth 

nex  nth 

°p , ex  ap , th 

1 

0. 

.603 

0. 

.718 

1. 

,160 

1, 

.238 

2. 

,210 

2. 

,089 

0. 

937 

0. 

.771 

- 

- 

2 

0. 

.881 

0. 

.774 

1. 

.450 

1. 

.252 

2. 

,290 

2. 

,027 

0. 

689 

0. 

.695 

- 

- 

3 

0. 

.632 

0. 

.627 

1. 

.170 

1, 

.041 

1. 

.910 

1. 

,714 

0. 

798 

0. 

.725 

0 

26 

0 

21 

4 

0. 

.471 

0. 

.535 

0. 

.870 

0. 

.908 

1. 

,700 

1. 

.530 

0. 

926 

0. 

.758 

- 

- 

5 

0. 

,737 

0. 

.721 

1. 

.160 

1, 

.168 

1. 

,740 

1. 

.770 

0. 

620 

0. 

.648 

0 

20 

0 

21 

6 

0. 

,680 

0. 

.704 

1. 

,100 

1, 

.150 

1. 

,770 

1. 

.832 

0. 

690 

0. 

.689 

- 

- 

7 

0. 

.631 

0. 

.624 

1. 

.090 

1. 

.025 

1. 

,820 

1. 

.644 

0. 

764 

0. 

.699 

- 

- 

8 

0. 

,500 

0. 

.531 

0. 

.901 

0, 

.893 

1. 

.600 

1. 

.461 

0. 

839 

0. 

.730 

- 

- 

9 

0. 

.676 

0. 

.618 

1. 

.030 

0, 

.967 

1. 

.590 

1. 

.387 

0. 

,617 

0. 

.583 

0 

23 

0 

21 

10 

0. 

637 

0. 

.601 

0. 

,978 

0. 

.933 

1. 

.490 

1. 

.342 

0. 

,613 

0. 

.580 

0 

26 

0 

20 

11 

0. 

.449 

0. 

.508 

0. 

.706 

0 

.799 

1. 

.170 

1. 

.158 

0. 

,691 

0. 

.594 

- 

- 

12 

0. 

.407 

0. 

.439 

0. 

.561 

0, 

.601 

0 

.761 

0, 

.797 

0. 

,451 

0. 

.430 

- 

- 

13 

0. 

.601 

0. 

.554 

0, 

.834 

0 

.767 

1, 

.160 

1. 

.035 

0. 

,474 

0. 

.451 

0. 

20 

0 

23 

14 

0. 

.521 

0. 

.522 

0. 

.718 

0 

.723 

0, 

.955 

0. 

.968 

0. 

,437 

0. 

.446 

- 

- 

15 

0. 

.536 

0. 

.505 

0 

.785 

0 

.688 

1, 

.120 

0. 

.923 

0. 

,532 

0, 

.435 

0. 

17 

0 

23 

16 

0. 

.368 

0 

.412 

0, 

.539 

0 

.555 

0 

.856 

0, 

.739 

0, 

.609 

0, 

.422 

0. 

32 

0 

25 

17 

0. 

.240 

0. 

.303 

0. 

.330 

0 

.401 

0, 

.436 

0. 

.508 

0. 

,431 

0. 

.371 

0. 

22 

0 

25 

18 

0. 

,375 

0. 

.384 

0. 

.524 

0. 

.526 

0, 

,708 

0. 

.695 

0. 

,458 

0. 

.428 

- 

- 

19 

0. 

.332 

0. 

.390 

0. 

.469 

0 

.531 

0, 

.630 

0. 

,692 

0. 

.462 

0. 

.414 

- 

- 

20 

0. 

.393 

0. 

.417 

0. 

,536 

0 

.573 

0, 

.732 

0. 

,755 

0. 

,365 

0. 

.428 

- 

- 

21 

0. 

.304 

0 

.346 

0. 

.445 

0 

.467 

0, 

.652 

0, 

.613 

0. 

,550 

0. 

.412 

0. 

25 

0 

29 
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Figure  29 

Theoretical  Versus  Experimental  Burning  Rate 
for  AP/HTPB  Propellants:  Type  II 
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AP/HTPB  PROPELLANTS 

10  PERCENT  ERROR  BAND 


EXPERIMENTAL  EXPONENT 


Figure  30 

Theoretical  Versus  Experimental  Exponent 
for  AP/HTPB  Propellants:  Type  II 
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AP/HTPB  PROPELLANTS 

1000  PSI  -  20  PERCENT  ERROR  BAND 


0.0  0.1  0.2  0.3  0.*  0.5 

EXP.  TEMPERATURE  SENSITIVITY  -  %/K 


Figure  31 

Theoretical  Versus  Experimental  Temperature 
Sensitivities  for  AP/HTPB  Propellants:  Type  II 
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4.2.3  Composite  Propellants  Containing  More  HMX  than  AP:  Type  III 

The  third  class  of  propellants  to  be  examined,  Type  III,  contain  both  HMX 
and  AP  crystals  dispersed  in  a  HTPB  binder  medium.  In  these  propellants  the 
percentage  of  HMX  is  between  60  and  78  percent  and  the  percentage  of  AP  is 
between  5  and  10  percent.  These  propellants  were  formulated  by  Hercules, 

Inc. [118]  For  the  series  of  13  propellants,  the  solids  loading  (amount  of 
oxidizer)  ranges  from  70  percent  to  78  percent.  The  particle  size  ranges  from 
4  to  200  microns  for  the  HMX  crystals  and  between  5  and  200  microns  for  the  AP 
crystals . 

Table  8  presents  the  formulation  data  for  these  HMX/AP/HTPB  propellants. 
Table  9  presents  the  experimental  data  and  the  theoretical  prediction  for  the 
burning  rate  at  400,  1000  and  1600  psia,  the  pressure  exponent  at  an  initial 
temperature  of  298"K  and  the  temperature  sensitivity  at  1000  psia.  Figures 
32,  33  and  34  plot  the  theoretical  HYPEM  values  versus  the  experimental  values 
for  the  burning  rates  at  1000  psi,  the  pressure  exponents,  and  the  temperature 
sensitivities,  respectively,  shown  in  Table  9.  As  shown  previously,  the  data 
points  numbered  and  error  bands  are  included. 

As  can  be  seen  in  Figure  32,  the  HYPEM  model  predicts  the  burning  rates 
for  the  majority  of  these  propellants  within  ±10  percent.  Figure  33  indicates 
that  the  model  predicts  the  pressure  exponent  between  +20  and  -10  percent. 
Figure  34  indicates  that  the  HYPEM  model's  burning  rate  temperature 
sensitivity  prediction  consistently  over  estimates  the  experimental  prediction 
by  more  than  20  percent. 


Table  8 

HMX/AP/HTPB  Propellant  Formulations  Containing  More  HMX  than  AP:  Type  III 


No 

D 

%HTPB  a 

%HMX 

%AP 

200/i 

147/i 

1.86 

58/i 

42.4/i 

2.68 

4/i 

4 . 46/i 
1.99 

200/i* 

200/i 

1.60 

50/i 

42.2/1 

1.96 

5  M 
4.49/1 
1.65 

01 

22.00 

_ 

78.00 

_ 

• 

_ 

_ 

02 

22.00 

- 

73.00 

- 

- 

5.00 

- 

03 

22.00 

- 

63.00 

- 

- 

10.00 

- 

04 

22.00 

- 

68.00 

- 

- 

- 

10.00 

05 

22.00 

- 

- 

68.00 

10.00 

- 

06 

22.00 

- 

34.00 

34.00 

- 

- 

10.00 

07 

22.00 

- 

- 

68.00 

- 

- 

10.00 

08 

22.00 

68.00 

- 

• 

10.00 

- 

- 

09 

30.0) 

- 

65.00 

- 

- 

5.00 

- 

10 

30.0) 

- 

60.00 

- 

- 

10.00 

- 

11 

30.0) 

- 

- 

60.00 

- 

- 

10.00 

12 

30.00 

- 

- 

65.00 

- 

- 

5.00 

13 

30.00 

- 

60.00 

- 

- 

- 

10.00 

Data  was  not  available  to  determine  the  distribution  data  of  the 
200/1  particle  size  HMX  crystals.  The  value  shown  is  assumed. 
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Table  9 

Predicted  and  Experimental  Burning  Rates,  Exponents  and 
Temperature  Sensitivities  for  HMX/AP/HTPB  Propellants 
Containing  More  HMX  than  AP:  Type  III 


Burning  Rate  -  in/sec 

Exponent 

Temperature 

Sensitivity 

400psia  lOOOpsia  1600psia 

%/°K 

No 

rex  rth  rex  rth  rex  rth 

nex  nth 

CTp > ex  ap,th 

1 

0 

061 

0 

074 

0 

112 

0 

116 

0 

117 

0 

144 

0 

496 

0 

486 

2 

0 

065 

0 

073 

0 

099 

0 

116 

0 

128 

0 

146 

0 

485 

0 

497 

0 

00 

r-4 

0 

34 

3 

0 

070 

0 

085 

0 

110 

0 

137 

0 

141 

0 

172 

0 

504 

0 

510 

- 

4 

0 

093 

0 

088 

0 

142 

0 

155 

0 

184 

0 

211 

0 

488 

0 

629 

0 

20 

0 

41 

5 

0 

099 

0 

121 

0 

167 

0 

203 

0 

214 

0 

266 

0 

558 

0 

565 

- 

- 

6 

0 

114 

0 

106 

0 

183 

0 

188 

0 

232 

0 

258 

0 

513 

0 

638 

0 

14 

0 

37 

7 

0 

137 

0 

124 

0 

223 

0 

221 

0 

276 

0 

305 

0 

509 

0 

644 

0 

34 

0 

35 

8 

0 

052 

0 

055 

0 

084 

0 

090 

0 

109 

0 

117 

0 

532 

0 

549 

- 

- 

9 

0 

043 

0 

050 

0 

070 

0 

081 

0 

092 

0 

106 

0 

546 

0 

538 

- 

10 

0 

049 

0 

059 

0 

077 

0 

097 

0 

098 

0 

126 

0 

499 

0 

547 

0 

25 

0 

41 

11 

0 

117 

0 

093 

0 

176 

0 

163 

0 

216 

0 

229 

0 

443 

0 

641 

0 

37 

0 

44 

12 

0 

090 

0 

087 

0 

137 

0 

144 

0 

174 

0 

199 

0 

473 

0 

593 

0 

31 

0 

39 

13 

0 

075 

0 

062 

0 

113 

0 

114 

0 

134 

0 

161 

0 

423 

0 

689 

0 

22 

0 

48 
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HMX/AP/HTPB  PROPELLANTS 

1000  PSI  -  10  PERCENT  ERROR  BAND 


Figure  32 

Theoretical  Versus  Experimental  Burning  Rate  for  HMX/AP/HTPB 
Propellants  Containing  More  HMX  than  AP:  Type  III 
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HMX/AP/HTPB  PROPELLANTS 

10  PERCENT  ERROR  BAND 


EXPERIMENTAL  EXPONENT 


Figure  33 

Theoretical  Versus  Experimental  Exponent  for  HMX/AP/HTPB 
Propellants  Containing  More  HMX  than  AP:  Type  III 
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HMX/AP/HTPB  PROPELLANTS 

1000  PSI  -  20  PERCENT  ERROR  BAND 


0.0  0.1  0.2  0.3  0.4  0.5 

EXP.  TEMPERATURE  SENSITIVITY  -  %/K 


Figure  34 

Theoretical  Versus  Experimental  Temperature  Sensitivities  for 
HMX/AP/HTPB  Propellants  Containing  More  HMX  than  AP:  Type  III 
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4.2.4  Composite  Propellants  Containing  More  AP  than  HMX:  Type  IV 

The  fourth  class  of  propellants,  Type  IV,  to  be  examined  contain  both  AP 
and  HMX  crystals  dispersed  in  a  HTPB  binder  medium.  In  these  propellants  the 
percentage  of  AP  is  between  72.37  and  77.37  percent  and  the  percentage  of  HMX 
is  between  10  and  15  percent.  These  propellants  were  formulated  by  Hercules, 
Inc. [29,30]  For  the  series  of  10  propellants,  the  solids  loading  was  constant 
at  87.37  percent.  The  particle  size  ranges  from  2  to  400  microns  for  the  AP 
crystals  and  between  4  and  20  microns  for  the  HMX  crystals. 

Table  10  presents  the  formulation  data  for  these  AP/HMX/HTPB  propellants. 
Table  11  presents  the  experimental  data  and  the  theoretical  prediction  for  the 
burning  rate  at  500,  1000  and  3000  psia,  the  pressure  exponent  at  an  initial 
temperature  of  298°K  and  the  temperature  sensitivity  at  1000  psia.  Figures 
35,  36  and  37  plot  the  theoretical  HYPEM  values  versus  the  experimental  values 
for  the  burning  rates  at  1000  psi,  the  pressure  exponents,  and  the  temperature 
sensitivities,  respectively,  shown  in  Table  11.  As  shown  previously,  the  data 
points  numbered  and  error  bands  are  included. 

As  can  be  seen  in  Figure  35  the  HYPEM  model  predicts  the  burning  rates  for 
the  majority  of  these  propellants  within  +15  to  -10  percent.  Figure  36 
indicates  that  the  model  over  predicts  the  pressure  exponent  by  20  percent. 
Figure  37  indicates  that  the  HYPEM  model's  burning  rate  temperature 
sensitivity  prediction  is  within  ±20  percent. 


Table  10 

AP/HMX/HTPB  Propellant  Formulations  Containing  More 
AP  than  HMX  -  12.63%  HTPB:  Type  IV 


%HMX 

%AP 

2  Op 

4p 

400p 

5  Op 

2  Op 

2n 

D 

22. 3p 

4.55  p 

436p 

48. 6p 

19. 3p 

1 . 96p 

No 

a 

2.89 

2.02 

1.16 

1.43 

1.80 

1.38 

1 

_ 

10.00 

42.11 

- 

13.68 

21.58 

2 

- 

10.00 

42.11 

21.58 

13.68 

- 

3 

10.00 

- 

42.11 

- 

13.68 

- 

4 

10.00 

- 

42.11 

21.58 

13.68 

- 

5 

10.00 

- 

- 

42.11 

13.68 

21.58 

6 

- 

15.00 

42.11 

- 

13.68 

16.58 

7 

- 

15.00 

42.11 

16.58 

13.68 

- 

8 

15.00 

- 

42.11 

- 

13.68 

16.58 

9 

15.00 

- 

42.11 

16.58 

13.68 

- 

10 

- 

15.00 

- 

42.11 

13.68 

16.58 
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Table  11 

Predicted  and  Experimental  Burning  Rates,  Exponents  and 
Temperature  Sensitivities  for  AP/HMX/HTPB  Propellants 
Containing  More  AP  than  HMX:  Type  IV 


Burning  Rate  -  in/sec 

500psia  lOOOpsia  3000psia 

No  rex  rth  rex  rth  rex  rth 

Exponent 

nex  nth 

Temperature 

Sensitivity 

%/°K 

°P , ex  °p , th 

1  0.483 

0.405 

0.630 

0.657 

1.120 

1.316 

0.417 

0.656 

_ 

_  _ 

2  0.302 

0.296 

0.374 

0.408 

0.627 

0.619 

0.368 

0.409 

0.23 

0.26 

3  0.397 

0.410 

0.523 

0.666 

1.210 

1.321 

0.536 

0.651 

4  0.271 

0.301 

0.348 

0.417 

0.599 

0.624 

0.403 

0.403 

0.26 

0.27 

5  0.558 

0.547 

0.759 

0.878 

1.340 

1.659 

0.528 

0.616 

0.20 

0.24 

6  0.544 

0.382 

0.704 

0.605 

1.260 

1.173 

0.445 

0.624 

7  0.291 

0.303 

0.359 

0.423 

0.595 

0.658 

0.343 

0.430 

0.25 

0.25 

8  0.342 

0.377 

0.450 

0.594 

0.943 

1.144 

0.472 

0.617 

0.33 

0.23 

9  0.259 

0.298 

0.329 

0.413 

0.582 

0.629 

0.395 

0.414 

0.22 

0.26 

10  0.609 

0.516 

0.775 

0.814 

0.963* 

1.238  0.332 

0.631 

0.15 

0.24 

The  Pressure  is  2000  psia  for  these  burning  rates 
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AP/HMX/HTPB  PROPELLANTS 

1000  PSI  -  10  PERCENT  ERROR  BAND 


0  0.2  0.4  0.6  0.8  1 

EXPERIMENTAL  BURNING  RATE  -  IN/SEC 


Figure  35 

Theoretical  Versus  Experimental  Burning  Rate  for  AP/HMX/HTPB 
Propellants  Containing  More  AP  than  HMX:  Type  IV 
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AP/HMX/HTPB  PROPELLANTS 

10  PERCENT  ERROR  BAND 


EXPERIMENTAL  EXPONENT 


Figure  36 

Theoretical  Versus  Experimental  Exponent  for  AP/HMX/HTPB 
Propellants  Containing  More  AP  than  HMX:  Type  IV 
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AP/HMX/HTPB  PROPELLANTS 

lOOO  PSI  -  20  PERCENT  ERROR  BAND 


0.0  O.l  0.2  0.3  0.4  0.5 

EXP.  TEMPERATURE  SENSITIVITY  -  %/K 


Figure  37 

Theoretical  Versus  Experimental  Temperature  Sensitivities  for 
HMX/AP/HTPB  Propellants  Containing  More  AP  than  HMX:  Type  IV 
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4 . 3  Summary  of  Combustion  Model  Results 

The  HYPEM  composite  combustion  model  computed  the  burning  rate  and 
pressure  exponent  for  61  different  composite  propellants.  In  addition,  model 
computed  burning  rate  temperature  sensitivity  for  30  individual  propellants. 
These  propellants  contain  HMX  and/or  AP  oxidizer  crystals  of  many  different 
particle  sizes  and  combinations  in  a  HTPB  inert  binder. 

Several  important  observations  can  be  made  about  the  behavior  of  the  model 
when  examining  Figures  26  through  37.  This  can  be  accomplished  by  comparing 
the  experimental  data  with  the  HYPEM  model's  prediction  and  examining  the 
propellant  formulations  presented  in  Tables  4,  6,  8  and  10.  The  following 
observations  are  categorized  by  propellant  type. 

1.  Type  I  propellants  which  contain  only  HMX  particles  (Tables  4  and  5, 
Figures  26,  27  and  28):  Burning  rate  and  pressure  exponent  prediction  is 
quite  good  except  for  propellants  containing  only  coarse  or  only  fine  HMX 
crystals.  Some  of  the  wide  distribution  propellant  properties  were  also 
poorly  predicted.  Wide  distribution  propellants  are  those  which  contain  very 
fine  and  very  coarse  particles.  The  predicted  temperature  sensitivity  is 
consistently  higher  than  the  measured  values  for  all  propellants. 

2.  Type  II  propellants  which  contain  only  AP  particles  (Tables  6  and  7, 
Figures  29,  30  and  31):  Burning  rate,  pressure  exponent  and  temperature 
sensitivity  prediction  were  all  reasonably  good.  Some  of  the  "wide 
distribution"  propellants  resulted  in  poor  but  acceptable  model  predictions 
that  were  usually  better  than  20  percent. 

3.  Type  III  propellants  which  contain  more  HMX  crystals  than  AP  crystals 
(Tables  8  and  9,  Figures  32,  33  and  34):  Burning  rate  prediction  was 
acceptable  but  the  pressure  exponent  prediction  was  higher  than  the 
experimental  values  for  propellants  with  fine  AP.  In  addition,  the  greater 
the  difference  between  the  HMX  and  fine  AP  particle  size  the  greater  the 
discrepancy.  As  with  Type  I  propellants  containing  only  HMX  particles,  the 
temperature  sensitivity  of  Type  III  propellants  is  over  predicted  by  the 
model . 

4.  Type  IV  propellants  which  contain  more  AP  crystals  than  HMX  crystals 
(Tables  10  and  11,  Figures  35,  36  and  37):  Burning  rate  prediction  is  usually 
within  20  percent  and  exhibits  no  noticeable  trends.  The  pressure  exponent  is 
over  predicted  for  those  propellants  with  fine  AP  or  "wide  distributions"  of 
AP.  The  HMX  particle  size  does  not  appear  to  be  a  factor.  Temperature 
sensitivity  prediction  is  usually  better  than  20  percent. 

In  general  the  HYPEM  model  does  a  satisfactory  job  in  predicting  burning 
rates  and  pressure  exponents.  As  with  many  other  models  of  this  type,  it  does 
have  a  difficult  time  with  propellants  containing  a  wide  distribution  of  HMX 
and/or  AP  particle  sizes.  In  addition,  the  HYPEM  model  over  predicts  the 
burning  rate  temperature  sensitivity  of  propellants  that  contain  all  or  mostly 
HMX. 


The  HYPEM  model's  difficulty  in  handling  "wide  distribution"  propellants 
is  due  to  the  fundamental  assumption  that  a  propellant  can  be  broken  down  into 
separate,  roonodisperse  pseudo  propellants.  A  consequence  of  this  assumption 
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is  that  different  particle  size  pseudo  propellants  burn  independently  of  each 
other.  Unaccounted  for  interdependencies  between  large  and  small  particles  is 
one  reason  the  HYPEM  model  poorly  predicts  combustion  characteristics  of  "wide 
distribution"  propellants.  Possible  interdependencies  include  flame,  heat 
transfer,  and  mass  interactions.  Some  of  these  effects  have  been  minimized 
and  the  model's  performance  improved  by  the  introduction  of  the  mass 
interaction  parameter  and  flame  temperature  corrections  discussed  previously 
in  Section  3. 1.6.1  and  3. 1.6. 2. 
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CHAPTER  V 
PARAMETRIC  STUDIES 


It  has  been  shown  that  the  HYPEM  composite  combustion  model  predicts  the 
burning  characteristics  of  a  wide  range  of  propellants  containing  HMX  and/or 
AP  crystals.  The  model  will  now  be  used  to  predict  the  effects  of  particle 
size,  chamber  pressure,  oxidizer  solids  loading  and  propellant  formulation 
changes  on  composite  propellant  burning  characteristics.  The  HYPEM  model's 
burning  rate,  pressure  exponent,  surface  temperature  and  temperature 
sensitivity  prediction  is  useful  tailoring  a  propellant  for  a  specific 
application.  In  addition,  the  systematic  variation  of  propellant  parameters 
provides  insight  to  the  controlling  mechanisms  of  propellant  combustion. 
Finally,  these  studies  will  address  the  limiting  conditions  of  the  combustion 
model . 


5 . 1  Controlling  Mechanisms 

Understanding  the  controlling  mechanisms  in  a  combustion  model  is 
important  in  order  to  know  the  limitations  of  the  model.  Figure  38  plots  the 
percent  heat  feedback  of  the  HMX  monopropellant  flame  and  the  primary  flame  to 
the  propellant  surface  versus  HMX  particle  size.  The  heat  transfer  from  the 
primary  flame  is  divided  into  that  portion  which  acts  on  the  HMX  crystal 
(dashed  line)  and  that  portion  which  acts  on  the  binder  (dotted  line).  For 
all  particle  sizes  the  primary  flame  energy  going  to  the  binder  is  dominant. 
This  is  due  to  the  large  amount  of  energy  required  to  heat  up  the  binder 
because  of  its  large  surface  decomposition  energy,  500  cal/g.  A  large 
percentage  of  the  binder  acts  as  an  ablative  material  that  does  not  react  with 
the  HMX  but  does  act  as  a  very  good  heat  sink.  For  small  HMX  crystals  the  HMX 
combustion  is  controlled  purely  by  the  primary  flame.  As  the  particle  size  is 
increased  the  HMX  monopropellant  flame  becomes  more  dominant  (solid  line) , 
eventually  dominating  the  HMX  crystal  combustion. 

A  second  way  of  examining  the  controlling  mechanisms  in  the  HMX  portion  of 
the  HYPEM  model  is  to  examine  the  relative  contribution  of  the  four  terms  in 
the  denominator  of  the  time  averaged  burning  rate  expression  for  the  HMX 
pseudo  propellants.  The  equation  is  written  below: 

(D'  +  2A) 

r - (166, 

'  D'  2A 

-  +  cign  +  frb -  +  a-frb)tb 

1  rX  rb  J 
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FLAME  HEAT  FEEDBACK  VS  HMX  PARTICLE  SIZE 

80%HMX/20%HTPB  1000PSI  298K  SIGMA=1 


1  10  100  1000 

HMX  PARTICLE  SIZE  -  MICRONS 


Figure  38 

Flame  Heat  Feedback  to  the  Propellant  Surface  Versus 
Oxidizer  Particle  Size  for  HMX  Propellants 
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In  Equation  (166)  D'  is  the  statistical  intersection  diameter  of  the  HMX 
crystal  and  A  is  the  distance  from  the  HMX  crystal  edge  to  the  edge  of  the 
assumed  binder  annulus.  The  terms  in  the  denominator  represent  the  oxidizer 
burning  rate,  rx,  oxidizer  ignition  delay  time,  t^  ,  binder  burn  through  time 
for  a  given  HMX  crystal,  r^,  and  the  binder  delay  time,  cb  •  f rb  is  the 
fraction  of  binder  that  reacts  with  the  HMX. 

The  contributions  in  percent,  of  the  four  physical  mechanisms  that  control 
the  burning  rate  are  illustrated  in  Figure  39.  The  contribution  of  each  of 
the  four  mechanisms  has  been  calculated  on  a  percentage  basis  and  plotted 
versus  the  HMX  particle  diameter.  The  HMX  crystal  burning  rate  term  (solid 
line)  is  the  greatest  for  fine  particles  but  the  other  mechanisms  contribute 
significantly.  At  10  microns,  the  ignition  delay  time  (dashed  line)  becomes 
the  largest  contributor  and  then  increases  until  it  is  completely  dominant  for 
100  micron  and  larger  particles.  For  the  larger  particle  sizes,  the  oxidizer 
rate  drops  to  less  than  10  percent  contribution.  The  binder  rate  (dotted 
line)  and  the  binder  burn  through  time  (dash-dotted  line)  contributions  remain 
roughly  constant  at  10  and  15  percent  respectively. 

Figure  40  plots  the  percent  heat  feedback  of  the  AP  monopropellant  flame, 
the  primary  flame,  and  the  final  flame  to  the  propellant  surface  versus  AP 
particle  size.  This  plot  shows  a  different  behavior  from  the  similar  plot, 
Figure  38,  for  the  HMX  pseudo  propellants.  From  Figure  40,  it  is  evident  that 
the  primary  flame  (dashed  line)  dominates  the  combustion  for  small  AP 
particles.  The  primary  flame  is  the  reaction  between  the  pyrolysis  products 
of  the  binder  surface  decomposition  and  the  AP  surface  decomposition  products. 
As  the  particle  size  increases  the  AP  monopropellant  flame  (solid  line)  gains 
in  importance.  At  the  same  time  the  final  flame  (dotted  line),  between  AP 
monopropellant  flame  products  and  binder  pyrolysis  products,  contributes  to 
the  surface  energy.  As  the  AP  particles  become  very  large,  the  AP 
monopropellant  flame  dominates  the  surface  energy  feedback.  The  obvious 
difference  between  Figures  38  and  40  is  the  way  the  primary  flame  is  treated. 
This  difference  arises  from  the  fact  that  in  the  HMX  portion  of  the  HYPEM 
model,  all  the  binder  does  not  react  with  the  HMX.  The  two  plots  are  similar 
in  that  the  oxidizer  monopropellant  flame  importance  increases  as  particle 
size  is  increased. 
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Figure  39 

Percent  Contribution  of  Terms  in  the  Denominator  of  the 
HMX  Rate  Equation 
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FLAME  HEAT  FEEDBACK  VS:  AP  PARTICLE  SIZE 

80%AP/20%HTPB  1000PSI  298K  SIGMA=1 


AP  PARTICLE  SIZE  -  MICRONS 


Figure  40 

Flame  Heat  Feedback  to  the  Propellant  Surface  Versus 
Oxidizer  Particle  Size  for  AP  Propellants 
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5 . 2  Particle  Size 

One  of  the  primary  factors  that  a  propellant  formulator  can  use  to  tailor 
a  propellant's  performance  is  the  size  of  the  oxidizer  particles  in  the 
propellant.  The  following  four  graphs  plot  the  burning  rate,  pressure 
exponent,  propellant  surface  temperature  and  temperature  sensitivity  versus 
particle  size  for  both  HMX  and  AP/HTPB  propellants.  The  solid  line  in  each  of 
these  curves  represents  propellants  that  are  88  percent  HMX  and  12  percent 
HTPB.  The  dashed  line  represents  propellants  that  are  88  percent  AP  and  12 
percent  HTPB.  For  these  figures,  and  the  ones  that  follow,  the  mode  width 
parameter,  crQX,  for  each  propellant  evaluated  to  create  the  curves  was  2.0. 


5.2.1  Burning  Rate 

Figure  41  plots  the  burning  rate  versus  particle  size.  In  Figure  41  the 
chamber  pressure  is  1000  psia  and  the  propellant's  initial  temperature  is 
298°K.  Figure  41  shows  that  the  burning  rate  for  HMX/HTPB  propellants  is  not 
a  strong  function  of  particle  size  except  for  small  HMX  particles.  Similar 
results  have  been  shown  experimentally. [30]  Referring  back  to  Figure  38  for 
small  HMX  crystals,  the  primary  flame  dominates  the  HMX  combustion.  For  these 
smaller  particle  sizes  the  relatively  cool  primary  flame  between  HMX  and  HTPB 
forces  a  relatively  slow  burning  rate.  Referring  back  to  Figure  39  for  large 
HMX  crystals,  the  HMX  crystal  ignition  delay  time  dominates  the  combustion  and 
also  forces  a  slow  burning  rate. 

Figure  41  shows  the  characteristic  "S"  shaped  burning  rate  curve  that  is 
experimentally  observed  for  AP/HTPB  propellants .[ 30]  Referring  back  to  Figure 
40  for  small  AP  crystals,  the  primary  flame  dominates  the  AP  combustion.  For 
these  smaller  particle  sizes  the  very  hot,  close  to  the  surface  primary  flame 
between  AP  and  HTPB  forces  a  relatively  fast  burning  rate.  Referring  back  to 
Figure  40  for  large  AP  crystals,  the  AP  monopropellant  dominates  the 
combustion  and  forces  a  slower  burning  rate.  As  is  evident  in  Figure  41,  a 
wide  variation  in  burning  rates  can  be  achieved  by  varying  the  AP  particle 
size . 


5.2.2  Pressure  Exponent 

Figure  42  plots  the  pressure  exponent  versus  particle  size.  The  pressure 
exponent  was  determined  by  evaluating  the  burning  rate  at  500,  1000  and  2000 
psi  by  using  a  least  squares  fit  of  the  rate  equation: 

r  -  cPn  ;18) 

where  n  is  the  pressure  exponent.  Figure  42  shows  a  gradual  decrease  in 
pressure  exponent  from  0.6  to  0.5  as  HMX  particle  size  is  increased.  The 
value  of  0.6  for  small  HMX  crystals  is  due  the  reaction  order,  1.15,  of  th 
kinetically  controlled  primary  flame  (see  Figure  38).  The  pressure  exponent 
is  approximately  1/2  the  value  of  the  reaction  order  for  kinetically 
controlled  flames.  The  decrease  in  pressure  exponent  as  particle  size 
increases  is  due  to  increased  diffusional  effects.  The  slight  increase  in 
exponent  for  very  large  HMX  crystals  is  due  to  HMX  monopropellant  kinetic 
flame  effects.  This  flame  has  a  reaction  order  of  1.39.  The  reason  the 


130 


BURNING  RATE  -  IN/SEC 


NWC  TP  6992 


BURNING  RATE  VS:  OXIDIZER  PARTICLE  SIZE 

88%0XIDIZER/12%HTPB  1000PSI  298K  SIGMA=2 


Figure  41 

Burning  Rate  Versus  Oxidizer  Particle  Size  for  both  HMX/HTPB 
and  AP/HTPB  Propellants 
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PRESSURE  EXPONENT  VS:  PARTICLE  SIZE 

88%OXIDIZER/12%HTPB  298K  SIGMA=2 


Figure  42 

Pressure  Exponent  Versus  Oxidizer  Particle  Size 
for  both  HMX/HTPB  and  AP/HTPB  Propellants 


132 


NWC  TP  6992 


exponent  only  increases  slightly  is  because  of  large  HMX  ignition  delay  times 
which  control  the  combustion  for  very  large  particles  (see  Figure  39) . 

Examining  Figure  42  for  the  AP  pressure  exponent  shows  that,  for  small  AP 
crystals  in  which  the  primary  flame  is  kinetically  controlled,  the  exponent  is 
high  due  the  high  reaction  order  associated  with  this  flame,  2.0.  As  the  AP 
particle  size  is  increased  the  final  diffusion  flame  and  the  diffusional 
effects  of  the  primary  flame  cause  a  decrease  in  the  pressure  exponent  (see 
Figure  40).  The  exponent  begins  to  rise  again  for  large  AP  particles  as  the 
kinetic  AP  monopropellant  flame  begins  to  dominate  the  combustion.  This  flame 
has  a  reaction  order  of  1.64. 


5.2.3  Propellant  Surface  Temperature 

The  next  particle  size  plot,  Figure  43,  shows  surface  temperature 
variations  with  differing  particle  sizes.  In  Figure  43  the  chamber  pressure 
is  1000  psia  and  the  propellant's  initial  temperature  is  298°K.  The  surface 
temperature  plotted  is  a  mass  weighted  average  of  the  individual  pseudo 
propellant  surface  temperatures  computed  by  the  combustion  model.  As 
expected,  these  trends  follow  the  same  pattern  as  the  burning  rate  variations 
shown  in  Figure  41.  The  HMX  propellants  show  little  variation  in  their 
surface  temperatures,  while  the  surface  temperature  varies  considerably  for 
the  AP  propellants.  In  addition,  the  AP  surface  temperatures  are  200  to  400°K 
higher  than  the  HMX  propellant  surface  temperature.  These  results  are 
consistent  with  experimental  observations . [63-65 ] 


5.2.4  Temperature  Sensitivity 

Figure  44  is  a  plot  of  the  temperature  sensitivity  versus  oxidizer 
particle  size.  The  temperature  sensitivity  w-r  compute^  by  computing  the 
burning  rate  at  initial  propellant  temperatures  of  248,  298  and  348°K.  The 
combustion  pressure  for  these  propellants  is  1000  psia.  Figure  44  indicates 
that  HMX  propellants  will  have  a  higher  temperature  sensitivity  than  will  AP 
propellants  and  that  both  types  will  show  an  increased  sensitivity  as  particle 
size  is  increased.  For  very  large  particles  the  AP  propellant  temperature 
sensitivity  surpasses  that  of  HMX.  It  has  been  shown  experimentally  that 
coarse  AP  particles  can  dramatically  increase  the  temperature  sensitivity,  and 
coarse  HMX  particles  cause  a  slight  increase  in  temperature  sensitivity .[ 155 ] 
In  the  section  comparing  experimental  data  with  the  HYPEM  model's  prediction 
(section  4.2),  the  HYPEM  over  predicted  the  temperature  sensitivity  when  the 
propellant  contained  all  or  mostly  HMX  as  a  oxidizer.  Based  on  those  results, 
the  temperature  sensitivity  values  of  the  HMX  propellants  in  Figure  44  are 
probably  20  to  40  percent  high. 
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SURFACE  TEMPERATURE  VS  PARTICLE  SIZE 

88%0XIDIZER/12%HTPB  1000PSI  298K  SIGMA=2 


Figure  43 

Propellant  Surface  Temperature  Versus  Oxidizer  Particle  Size 
for  both  HMX/HTPB  and  AP/HTPB  Propellants 
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TEMPERATURE  SENSITIVITY  VS  PARTICLE  SIZE 

88%OXIDIZER/l2%HTPB  1000PSI  SIGMA=2 


Figure  44 

Temperature  Sensitivity  Versus  Oxidizer  Particle  Size 
for  both  HMX/HTPB  and  AP/HTPB  Propellants 
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5 . 3  Pressure 


To  examine  a  propellant's  behavior  as  particle  size  changes  requires  the 
use  of  monomodal  particle  size  propellants,  ie .  propellants  containing  only 
one  particle  size  distribution.  In  reality  however,  in  order  to  meet  the  high 
solids  loading  of  today's  modern  propellants,  several  distribution  of  particle 
sizes  must  be  used.  A  realistic  propellant  might  consist  of  fine  20/im 
particles  and  coarse  200/im  particles.  Figures  45  and  46  show  the  burning  rate 
and  surface  temperature  versus  pressure  for  propellants  that  are  88  percent 
solids  loaded  containing  a  50/50  mix  of  fine  20/im  particle  and  coarse  200/im 
particles.  In  Figure  45  and  46  the  initial  propellant  temperature  is  298°K 
and  the  particle  distributions  have  a  sigma  of  2.0.  Four  types  of  propellants 
are  shown.  The  solid  line  represents  propellants  containing  only  HMX 
crystals,  while  the  dashed  line  represents  propellants  containing  just  AP 
crystals.  The  remaining  two  lines  depict  propellants  containing  fine  HMX 
mixed  with  coarse  AP  (dotted  line)  and  fine  AP  with  coarse  HMX  crystals 
(dotted-dash  line).  Since  Figure  45  is  plotted  on  a  log-log  scale  the  slope 
of  the  curve  is  the  pressure  exponent. 

All  propellants  in  Figures  45  and  46  show  an  increase  in  burning  rate  and 
surface  temperature  as  pressure  increases.  Increasing  the  pressure  causes  an 
increase  in  the  flame  heat  transfer  to  the  surface  since  the  various  flames 
are  closer  to  the  surface.  The  propellants  with  only  HMX  have  the  lowest 
burning  rates  and  surface  temperatures.  Conversely,  the  AP  propellants  have 
the  highest  burning  rates  and  surface  temperatures.  The  two  combination 
propellants  lie  in  between.  The  two  propellants  with  fine  HMX  (bottom  two 
lines)  and  the  two  propellants  with  fine  AP  (top  two  lines)  both  behave  in  a 
similar  fashion. 

Referring  back  to  Figures  38  and  40  which  plot  the  percent  heat  feedback 
from  the  various  flames  versus  particle  size,  it  was  concluded  that  the 
primary  flame  controls  the  combustion  for  both  HMX  and  AP  fine  particles.  For 
the  propellants  in  Figure  45  with  fine  HMX  crystals,  the  cool  primary  flame 
suppresses  the  rate  below  that  of  the  AP  fine  crystal  propellants  which  have  a 
hot  primary  flame.  As  the  pressure  increases  the  flame  standoff  distance  of 
the  primary  flame  becomes  even  smaller.  Because  of  the  very  hot  AP  primary 
flame,  propellants  with  fine  AP  will  increase  in  burning  rate  faster  than 
propellants  containing  fine  HMX  as  is  shown  in  Figure  45. 

Figure  47  shows  the  effect  of  pressure  on  temperature  sensitivity  for  the 
same  bimodal  propellants  discussed  In  Figures  45  and  46.  Each  propellant 
contains  a  50/50  mix  of  200/20/m  particles  with  12  percent  HTPB.  The  values 
in  the  Figure  47  at  1000  psi  agree  well  with  observed  temperature  sensitivity 
in  the  literature .[ 155-159 ]  Figure  47  shows  the  least  amount  of  temperature 
sensitivity  fluctuations  with  the  HMX/HTPB  propellants  (solid  line). 
Propellants  containing  all  or  some  AP  show  more  temperature  sensitivity 
variations  as  pressure  is  increased.  These  trends  are  partially  supported  by 
the  literature  Experimental  studies  indicates  the  effect  of  pressure  on 
temperature  sensitiviv  has  minimal  effect  on  HMX/HTPB  propellants  ard 
generally  decreases  or  results  in  little  change  in  the  temperature  sensitivity 
of  AP/HTPB  propellants. 
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BURNING  RATE  VS:  PRESSURE 
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Figure  45 

Burning  Rate  Versus  Pressure  for  Bimodal  Propellants 
Containing  HMX  and/or  AP  Particles  (12%  HTPB) 
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SURFACE  TEMPERATURE  VS:  PRESSURE 

50/50  20um/200um  PARTICLES  298K  SIGMA=2 


10  100  1000  10000 

PRESSURE  -  PSI 


Figure  46 

Surface  Temperature  Versus  Pressure  for  Bimodal  Propellants 
Containing  HMX  and/or  AP  Particles  (12%  HTPB) 
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TEMPERATURE  SENSITIVITY  VS:  PRESSURE 
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Figure  47 

Temperature  Sensitivity  Versus  Pressure  for  Bimodal 
Propellants  Containing  HMX  and/or  AP  Particles  (12%  HTPB) 
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5.4  Oxidizer  Solids  Loading 

The  next  series  of  graphs  show  the  combustion  characteristics  of  composite 
propellants  as  a  function  of  oxidizer  mass  fraction  or  solids  loading.  The 
solids  loading  ranges  from  60  percent  (fuel  rich)  to  98  percent  (fuel  lean)  in 
each  figure.  Each  graph  contains  three  lines  which  represent  monomodal 
propellants  containing  fine  10^m,  medium  100/im  and  coarse  400/«n  oxidizer 
particles.  The  particle  distributions  have  a  sigma  of  2.0. 

Figures  48  and  49  show  burning  rate  versus  solids  loading  for  HMX/HTPB 
propellants  and  AP/HTPB  propellants,  respectively.  In  both  graphs  the 
pressure  is  1000  psi  and  the  initial  propellant  temperature  is  298 °K.  Figure 
48  for  HMX/HTPB  propellants  indicates  that  the  higher  the  solids  loading  the 
higher  the  burning  rate  or,  conversely,  the  addition  of  any  HTPB  to  the 
propellant  will  lower  the  burning  rate.  Behavior  of  this  type  is  expected 
since  HMX  burns  stoichiometrically  as  a  monopropellant,  and  the  addition  of 
fuel  will  dilute  and  weaken  the  reaction. [30]  As  expected  from  previous 
theoretical  and  experimental  results,  the  differing  particle  sizes  have  little 
influence  on  the  burning  rate. 

The  AP/HTPB  curves  in  Figure  49  show  a  peak  at  approximately  90  percent 
solids  loading  for  the  10/im  AP  particle  curve  and  at  95  percent  solids  loading 
for  the  100/m  curve.  A  peak  is  expected  since  AP  requires  fuel  to  burn  in 
stoichiometric  proportions.  The  stoichiometric  solids  loading  for  an  AP/HTPB 
mixture  is  close  to  88  percent.  The  reason  the  burning  rate  peak  is  not  seen 
for  the  coarser  particles  is  that,  for  these  larger  particles,  more  of  the 
energy  is  released  in  the  final  flame  while  the  primary  flame  becomes  less 
important.  The  final  flame  has  a  large  flame  standoff  distance  so  less  energy 
can  be  conducted  to  the  surface  of  the  propellant. 

Figures  50  and  51  depict  the  pressure  exponent  variation  with  differing 
solids  loading  for  HMX/HTPB  and  AP/HTPB  propellants.  The  trends  presented 
agree  with  the  previous  theoretical  results.  The  pressure  exponent  of 
HMX/HTPB  propellants  is  only  affected  slightly  by  changes  In  propellant 
characteristics.  The  AP/HTPB  propellants  in  Figure  51  show  much  broader 
variation  in  the  pressure  exponent.  All  three  particle  size  propellants  at 
very  high  solids  loading  approach  their  monopropellant  pressure  exponents, 
since  the  combustion  processes  are  dominated  by  the  kinetic  monopropellant 
oxidizer  flames. 

The  final  two  graphs  in  this  series.  Figure  52  and  53,  plot  the 
temperature  sensitivity  variations  with  differing  oxidizer  solids  loadings. 

In  both  figures  the  higher  the  solids  loading  the  lower  the  temperature 
sensitivity.  The  figures  indicate  fuel  rich  propellants  will  have  a  high 
burning  rate  sensitivity  to  initial  propellant  temperature.  When  the 
conditions  are  closer  to  stoichiometric  conditions,  solids  loading  around  90 
percent,  the  propellants  containing  the  most  coarse  particles  have  the  highest 
temperature  sensitivity.  This  agrees  with  the  results  already  presented 
showing  temperature  sensitivity  variation  with  different  particle  sizes 
(Figure  42) . 
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Burning  Rate  Figure  48 

“  —  —A  for  Mononodal 


1 


BURNING  RATE  -  IN/SEC 


NWC  TP  6992 


BURNING  RATE  VS:  SOLIDS  LOADING 

AP/HTPB  1000PSI  298K  SIGMA=2 


SOLIDS  LOADING  -  %AP 


Figure  49 

Burning  Rate  Versus  Oxidizer  Solids  Loading  for  Monomodal 
AP/HTPB  Propellants 
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PRESSURE  EXPONENT  VS:  SOLIDS  LOADING 
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Figure  50 

Pressure  Exponent  Versus  Oxidizer  Solids  Loading  for 
Monomodal  HMX/HTPB  Propellants 
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PRESSURE  EXPONENT  VS:  SOLIDS  LOADING 
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Figure  51 

Pressure  Exponent  Versus  Oxidizer  Solids  Loading  for 
Monomodal  AP/HTPB  Propellants 
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TEMP.  SENSITIVITY  VS:  SOLIDS  LOADING 
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Figure  52 

Temperature  Sensitivity  Versus  Oxidizer  Solids  Loading  for 
Mono'modal  HMX/HTPB  Propellants 
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TEMP.  SENSITIVITY  VS:  SOLIDS  LOADING 
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Figure  53 

Temperature  Sensitivity  Versus  Oxidizer  Solids  Loading  for 
Monomodal  AP/HTPB  Propellants 
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5 . 5  Formulation  Changes 

The  final  series  of  plots  presents  the  results  of  calculations  in  which 
the  propellant  formulation  was  changed.  The  following  three  figures  show  the 
effect  of  coarse  to  fine  particle  ratio  on  burning  rate,  pressure  exponent  and 
temperature  sensitivity  for  bimodal  propellants  consisting  of  200/im  and  20/jm 
oxidizer  particles.  In  each  plot  the  bimodal  propellant  property  is  plotted 
versus  percent  200pm  particles.  The  far  left  edge  of  the  plot,  0%  coarse, 
implies  the  propellant  consists  of  all  20pm  fine  particles.  The  far  right 
edge  of  the  plot,  100%  coarse,  implies  the  propellant  consists  of  all  200pm 
coarse  particles.  In  addition,  on  each  plot,  four  types  of  propellants  are 

presented.  The  solid  line  represents  propellants  containing  only  HMX 

particles,  the  dashed  line  represents  propellants  containing  only  AP 

particles,  the  dotted  line  represents  propellants  containing  coarse  HMX  and 

fine  AP  and  the  dash-dotted  line  represents  propellants 

containing  coarse  AP  and  fine  HMX.  The  mode  width  parameter,  OqX,  for  each 
oxidizer  distribution  is  2.0. 

Figure  54  shows  the  burning  rate  behavior  versus  the  percent  of  coarse 
particles.  The  propellants  containing  fine  AP  show  the  greatest  variation  in 
burning  rate  with  the  fine  AP/coarse  HMX  propellants  yielding  the  greatest 
variation  (dotted  line).  Conversely,  propellants  containing  fine  HMX  crystals 
show  the  least  variation  in  burning  rate  with  the  fine  HMX/coarse  HMX 
propellants  yielding  the  least  variation  (solid  line). 

Figure  55  shows  the  pressure  exponent  variation  for  the  same  propellants. 
The  propellants  containing  coarse  HMX  have  the  highest  exponent  while  those 
with  coarse  AP  have  the  lowest  exponent.  The  amount  of  fine  AP  crystals 
combined  with  coarse  HMX  in  a  propellant  causes  the  widest  variation  in 
burning  rate.  This  type  of  propellant  also  yields  the  highest  exponent  as  is 
seen  in  Figure  55.  Unfortunately  these  two  facts  run  counter  to  each  other  in 
the  design  of  well  behaved  propellants,  ie.  fine  AP/coarse  HMX  mixtures  yield 
varied  burning  rates  but  also  the  highest  pressure  exponent. 

Figure  56  plots  the  temperature  sensitivity  versus  percent  coarse 
oxidizer.  In  Figure  56,  the  temperature  sensitivity  is  the  highest  for 
propellants  with  the  most  coarse  HMX  crystals.  Propellants  containing  mostly 
fine  AP  particles  will  have  the  least  burning  rate  sensitivity  to  initial 
propellant  temperature.  These  observations  indicate  that  the  very  hot  primary 
flame  associated  with  small  AP  crystals  is  relatively  insensitive  to  initial 
propellant  temperature.  Furthermore,  the  raonopropellant  kinetic  flame 
associated  with  larger  particles  is  more  sensitive  to  initial  propellant 
temperature . 
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Figure  54 

Burning  Rate  Versus  Percent  Coarse  Oxidizer  Crystals  for 
Bimodal  Propellants  Containing  HMX  and/or  AP  with  HTPB 
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Figure  55 

Pressure  Exponent  Versus  Percent  Coarse  Oxidizer  Crystals  for 
Bimodal  Propellants  Containing  HMX  and/or  AP  with  HTPB 
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Figure  56 

Temperature  Sensitivity  Versus  Percent  Coarse  Oxidizer  Crystals 
for  Bimodal  Propellants  Containing  HMX  and/or  AP  with  HTPB 
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CHAPTER  VI 
CONCLUSIONS 


In  this  thesis  the  High  Energy  Petite  Ensemble  Model  (HYPEM)  has  been 
presented.  The  model  is  a  steady  state,  statistically  based  composite 
propellant  combustion  model  which  is  able  to  predict  the  combustion 
characteristics  of  propellants  consisting  of  HMX  and/or  AP  mixed  together  in  a 
HTPB  fuel  binder  medium.  The  overall  structure  of  the  combustion  model  is 
based  on  the  Petite  Ensemble  Model  (PEM)  which  mathematically  separates  the 
propellant  into  unique  pseudo  propellants.  These  monodisperse  pseudo 
propellants  depend  upon  the  statistical  makeup  of  the  individual  oxidizer 
particle  distributions  in  the  propellant.  The  pseudo  propellant  burning  rates 
are  then  statistically  combined  together  to  determine  the  total  burning  rate 
of  the  propellant. 

The  HYPEM  combustion  model  treats  the  two  oxidizer  types,  HMX  and  AP,  in 
separate,  unique  fashions.  The  burning  rate  of  the  HMX  pseudo  propellants  is 
computed  by  averaging  the  HMX  combustion,  HMX  ignition  delay,  binder  burn 
through  and  binder  delay  time  as  the  propellant  surface  regresses.  Time 
averaging  is  required  due  to  the  unsteady  nature  of  HMX  crystals  burning  in  a 
composite  propellants.  The  burning  rate  of  the  AP  pseudo  propellants  is 
computed  by  space  averaging  over  the  entire  burning  surface  of  the  pseudo 
propellant.  Both  the  HMX  and  AP  pseudo  propellants  use  surface  energy 
balances  to  compute  the  surface  temperature  of  the  individual  oxidizer  type 
and  its  surrounding  binder. 

The  HYPEM  composite  propellant  combustion  model  was  used  to  predict  the 
burning  rate  and  the  pressure  exponent  for  61  individual  propellants 
containing  one  or  both  types  of  oxidizers.  These  propellants  contained  HMX 
particles  from  4jim  to  400pm  and  AP  particles  from  0.7pm  to  450pm.  The  burning 
rate  and  pressure  exponent  results  were  most  accurate  for  propellants 
containing  only  one  type  of  oxidizer,  however,  the  results  for  the  combination 
propellants  were  usually  within  ±20  percent.  The  HYPEM  model  was  also  usee,  to 
com >ute  the  temperature  sersitivity  of  30  individual  propellants. 
Experimentally  measured  values  of  temperature  sensitivity  compared  favorably 
with  the  predictions  by  the  combustion  model  for  propellants  containing  mostly 
or  ill  AP.  The  HYPEM  model  consistently  overpredicted  the  temperature 
sensitivity  for  propellants  containing  all  or  mostly  HMX  by  20  to  40  percert. 

Various  parametric  studies  were  performed  with  the  HYIEM  combustion  model. 
The  experimentally  observed  small  burning  rate  and  pressure  exponent 
dependence  on  HMX  oxidizer  particle  size  was  correctly  predicted  by  the  model. 
In  addition,  the  experimentally  observed  large  burning  rate  dependence  on  AP 
particle  size  was  shown  by  the  model.  The  experimental  observation  that 
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larger  HMX  and  AP  crystals  will  increase  the  temperature  sensitivity  was  also 
verified  by  the  model.  Parametric  studies  showing  the  effect  of  pressure  on 
burning  rate,  surface  temperature  and  temperature  sensitivity  were  also 
presented. 

Parametric  studies  were  performed  examining  the  effect  of  propellant 
solids  loading  and  th6  ratio  of  coarse  to  fine  oxidizer  crystals  on  bimodal 
propellant  burning  rate,  pressure  exponent  and  temperature  sensitivity.  The 
bimodal  propellant  results  indicated  a  strong  dependence  on  temperature 
sensitivity  for  fine  AP/coarse  HMX/HTPB  propellants  but  little  dependence  with 
HMX/HTPB  propellants.  Considerable  burning  rate  variations  could  be  achieved 
by  varying  the  ratio  of  coarse  HMX  with  fine  AP  oxidizer  crystals.  However, 
these  same  propellants  exhibited  the  highest  pressure  exponent.  Studies  like 
these,  and  similar  studies  where  the  formulation  is  systematically  changed, 
can  be  used  to  understand  propellant  behavior  and  can  aid  in  the  development 
of  new  propellants. 
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APPENDIX 


Sample  Computer  Printout 

A  sample  computer  printout  is  presented  which  represents  a  test  run  of  the 
HYPEM  computer  code.  In  this  test  case  the  computer  propellant  consists  of  44 
percent  HMX,  44  percent  AP  and  12  percent  HTPB.  The  HMX  crystals  consist  of  a 
single  distribution  of  particles  with  a  mean  particle  size  of  200/im  with  a 
sigma  of  2.2.  The  AP  crystals  consist  of  a  single  distribution  of  particles 
with  a  mean  particle  size  of  50/im  with  a  sigma  of  1.8.  The  propellant  burning 
rate  is  evaluated  at  500,  1000  and  2000  psia  and  at  initial  propellant 
temperatures  of  248.15,  298.15  and  348.15°K.  The  first  page  of  the  output 
lists  all  HYPEM  model  inputs  and  initial  conditions.  After  calculations  are 
completed  for  a  particular  initial  propellant  temperature  at  all  three 
pressures,  a  summary  is  printed  which  lists  the  burning  rate  at  the  three 
pressures  along  with  the  computed  pressure  exponent.  After  calculations  are 
completed  for  all  three  initial  propellant  temperatures,  each  at  all  three 
pressures,  a  summary  is  printed  which  lists  the  burning  rates,  pressure 
exponents  at  all  three  initial  propellant  temperatures  and  the  temperature 
sensitivity  at  each  of  the  three  pressures.  Other  important  output  includes 
the  various  surface  temperatures  and  flame  heights.  For  the  HMX  pseudo 
propellants  the  percent  contribution  from  the  various  terms  in  the  time 
averaged  burn  rate  equation  are  output  and  for  the  AP  pseudo  propellants  the 
percent  heat  transfer  to  the  surface  from  the  AP  flame,  primary  flame  and 
final  flame  are  output. 
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'PSEUDO  PROPELLANT  INPUT  PARAMETERS  for  at l  2  MOOE(S)  for  all  PRESSURES  and  TEMPERATURES 
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THE  BURN  RATE  for  this  MODE  is  0.2763  cm/sec  or  0.1088  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  727.11 
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THE  BURN  RATE  for  this  HOOE  is  0.3840  cm/sec  or  0.1512  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  737.04  K 
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THE  BURN  RATE  for  this  HOOE  is  0.5455  cm/sec  or  0.2148  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  745.56  K 
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'  MODE  1  AP***  23  PARTICLES,  the  PRESSURE  IS  500.00  psia  and  the  INITIAL  PROPELLANT  TEMPERATURE  is  298.15  DEGREES  KELVIN 

D  R  TSOX  TSB  TS  8F  XAPK  XPFK  XPFD  XFFD  HTAP  HTPF  HTFF  TF  TSQDAP  TSQDB  ITER 

: CRONS  cm/s  K  K  K  . MICRONS .  X  HEAT  TRANSFER  K  K  K 
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THE  BURN  RATE  fr  this  MOOE  is  0.3335  cm/sec  or  0.1313  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  734.37  K 
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THE  BURN  RATE  for  this  MODE  is  0.4635  cm/sec  or  0.1825  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  744.27  K 
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THE  BURN  RATE  for  this  MODE  is  0.6571  cm/sec  or  0.2587  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  752.22  (C 


FINAL  RESULTS  FOR  AN  INITIAL  PROPELLANT  TEMPERATURE  OF  298.15  K 
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'*  "JOE  1  AP**#  25  PARTICLES,  the  PRESSURE  IS  500.00  psia  and  the  INITIAL  PROPELLANT  TEMPERATURE  Is  348.15  DEGREES  KELVIN 
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MICRONS  cm/S  K  K  K  . MICRONS .  X  HEAT  TRANSFER  K  K  K 
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THE  BURN  RATE  for  this  MOOE  is  0.4069  cm/sec  or  0.1602  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  T40.91  K 
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THE  8URN  RATE  for  this  MOOE  is  0.5656  cm/sec  or  0.2227  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  750.75  K 


MOOC  1  AP***  23  PAnI 1CLES,  the  PRESSURE  IS  2000.00  psia  and  the  INITIAL  PROPELLANT  TEMPERATURE  is  348.15  DEGREES  KELVIN 
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THE  BURN  RATE  for  this  MOOE  is  0.7994  ;ir/sec  or  0.3147  in/sec  and  the  AVERAGE  SURFACE  TEMPERATURE  is  758.11  K 
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